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Summary 

A simple model is set up of a star enclosed by an expanding envelope. The 
star is emitting a continuous spectrum and the envelope scatters radiation in 
the neighbourhood of a resonance frequency. The profiles which result 
from these assumptions are of the emission-line type. In varying degrees 
they show violet absorption borders, central reversals and occultation of the 
red edge, which are characteristic features of lines in P Cygni spectra. 

The profiles are obtained by numerical solution of a pair of integral 
equations for the emissivity in the envelope and the flux of radiation from 
the system. ‘The critical parameter for the appearance of central reversal 
is the ratio of expansion velocity to turbulent velocity in the envelope. 
(v/v* ~ 3.) Noticeable occultation of the emission component sets in when 
the effective radii of envelope and star fall below the ratio R,/Ry ~ 1°25. 

For profiles in which the net energy in the line exceeds the underlying 
continuous spectrum, the Zanstra recombination cycle is introduced. The 
computed contours show, in addition to an increase in intensity of the 
emission components, an enhancement of the associated central reversal. 
In some cases this results in obscuring the violet absorption border. 

The theoretically derived profiles encompass a wide range of lines, not 
only of P Cygni type but of some types of Be emission lines as well. 





1. Introduction 

The existence of an extensive atmosphere subject to large-scale motions is 
now widely accepted in interpreting the emission-line contours of early-type stars. 
Envelopes of this type have been postulated as early as 1904—by Halm in a paper 
read to the Royal Society of Edinburgh—and are now used in discussions of the 
spectra of planetary nebulae*, novae after outburst (2, 3), Wolf-Rayet stars 
(4), and the stars of P Cygni class (§). Since these macroscopic motions are of 
a scale larger than can reasonably be ascribed to thermal agitation or turbulence, 
one infers that they take the form of a radial expansion. For planetary nebulae, 
novae and W-R stars this model satisfactorily explains most of the observed 


features. 
* e.g. O. C, Wilson (1). 
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In the case of P Cygni lines, however, the profiles are often highly complex. 
Previous analyses of their features by straightforward application of a radially 
expanding atmospheric model, while successful for the simpler profiles, have 
for the more detailed ones seemed to require the introduction of a number of 
ad hoc assumptions concerning the acceleration of this motion at various heights 
in the envelope. In particular the often present ‘‘ central reversals’’ have been 
taken as ground for inferring a sharp deceleration of flow high up in the envelope. 
Interpretations of this type* throw the burden of explanation on the mechanical 
aspects of the flow, without examining the necessity of the additional postulates. 

The point of view adopted in this paper is that the contours of the emission 
lines can be explained in terms of simple, somewhat idealized systems, if a close 
rein is kept on the origin and destination of the radiation under discussion. 
General equations are set up and the profiles computed in a number of particular 
cases to show that the characteristic features of the complex P Cygni profiles 
are consequences of the coupling of the Doppler displacements (resulting from the 
radial expansion) with the geometry assumed in the model. Special assumptions 
about the velocity gradient are not required. It is however necessary to assign 
a significant role to resonance scattering, thereby making more precise the 
mechanisms of absorption and emission which operate to form a line contour. 


2. Atmospheric models and their mathematical idealization 

Failing any evidence to the contrary we make the obvious assumptions through- 
out the following that the stars with which we are dealing are each 

(1) spherically symmetric systems in their physical and dynamical structure ; 

(2) describable by variables which are time-independent (that is, they are in 
a “‘ steady state’’). 

In view of the above assumptions the only mechanisms that one can now 
allow for producing an emission line are those which concentrate energy flux 
into specific frequencies by depletion—since radiation transfer is subject to 
energy conservation—of other frequencies of the energy spectrum. For an 
absorption line one has the reverse. ‘There are two such frequency-changing 
processes which are recognized as being operative in the atmospheres of early- 
type stars: 

(1) Matter, which absorbs and re-emits light in the same frequency but in 
different directions as seen from a frame of reference in which the matter is at 
rest, changes the frequency (as well as the direction) of a ray when the scattering 
is viewed by an observer with respect to whom the matter is in motion (8). This 
is the manner in which the Doppler effect enters our calculations. 

(2) Absorption and re-emission by a process of photoelectric ionization 
and recombination convert radiation of high frequency into radiation of lower 
frequencies. ‘This is the Zanstra mechanism (9g), according to which (in the 
microscopic view) there is a strong probability of a quantum of radiant energy 
absorbed beyond the series limit reappearing after absorption in the form of 
quanta in the lines of the series—due to the possibility of electron capture at any 
excited level, followed by a sequence of cascade jumps down to the ground state. 
(One must of course admit the possibility of a similar conversion of quanta from 
lines of higher excitation into lines of lower excitation.) 

* vide, e.g., (6) and (7). 
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These two processes are not mutually exclusive and there is no method as 
yet known for deciding when one or the other is not operative in the formation of 
a particular line. When, however, the equivalent width of the emission com- 
ponents of a complex line is greater than the equivalent width of the adjacent 
absorption components one may safely conclude that process (1) cannot provide 
an adequate explanation and that therefore process (2) must be at work. ‘The 
mere existence of a small emission component alongside a deeper absorption 
component does not however in itself imply recombination, as has usually been 
assumed. Likewise, when the absorption component is well displaced to the 
violet or is greatly broadened one may infer that process (1) cannot be ignored. 

For the representative case our simplified model takes the form of a perfectly 
absorbing body (bounded roughly by the reversing layer) which is radiating a 
continuous spectrum. We call this the star proper, or more simply the star. 
Enclosing the star we conceive of an expanding envelope, with arbitrary velocity 
gradient, and in which the ionization and excitation decrease radially outward 
from the reversing layer. ‘The outermost regions of this envelope will be the 
ones chiefly active in resonance scattering; whereas the Zanstra mechanism, 
when it is operative, will be at work in the regions closer to the star. ‘The 
distinction between these regions of operation is only qualitative, and at best 
poorly defined. ‘The recombining region, when it exists at all, may be expected 
to overlap and blend into the scattering region. Again, the field of radiation 
from the star may be so intense that for a considerable distance outward from the 
inner boundary of the envelope the temperature of the electron-ion assemblage 
is so high that very little recombination takes place at all. 

Clearly the simpler situation (1 above) is the one in which, for the line con- 
cerned, no new radiation of that frequency enters the stream beyond the surface 
of the star proper; that is to say, when there is no region of recombination. 
This case will be treated first and, as will be seen, the region of scattering will be 
idealized still further. The residual generality will still be sufficient to provide 
for most of the distinguishing characteristics in the typical P Cygni profile. 
The more complicated system, in which there is a recombining region lying 
between the surface of the star and the region of scattering, will be taken up at the 
end of the paper, and will be shown to yield P Cygni profiles in which the emission 
components are of great intensity relative to the absorption. ‘The case of a purely 
recombining envelope, not enclosed by any region cool enough to permit 
resonance scattering, is probably more characteristic of the hotter Wolf—Rayet 
stars and will be given only passing reference here. The profiles for this latter 
case have already been derived, using comparatively simple methods (2, 10, 11, 12). 

The tracing of the alterations in a given radiation field, in the course of its 
passage from the inner boundary of the envelope out te a distant observer, is 
properly a problem of solving an equation of transfer (8) subject to the appropriate 
boundary conditions, velocity distribution and structural symmetry. However, 
the difficulties encountered in treating even the case of plane-parallel stratification, 
in the lowest approximation of directional dependence and velocity gradient 
(13, 14), have been considerable. For the more involved case of a spherically 
symmetric envelope they would be well-nigh insurmountable. In order to 
avoid having to treat a transfer problem, and yet preserve the significant effect 
of scattering by a moving medium, we shall suppose that the envelope may be 

10* 
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treated as geometrically very thin although remaining optically dense. Thus let 
x, be the scattering coefficient per unit mass for frequency v, p(r) the density 
at distance r from the centre, R’ the interior radius of the envelope and R” the 
exterior radius. Then the optical thickness in the radial direction is 


a 
— «,p(r) dr. (1) 
/ R’ 


We assume that this quantity remains finite as we compress the envelope into a 
thin shell, R’>R”=R. Consequently, for a ray making an angle cos"! » with 
the normal to the R’ surface, the optical thickness 
oa rdr 
r) ds ds = ; 2h 
Jp eas (Rr?) 


becomes in the limit 





7, / (2) 
as R’-> R”. ‘'Thisis simply the expression for the optical thickness in any direction 
of a laminar medium. In like fashion, if B,(r) be the emissivity (Ergiebigkeit) 
at a point in the envelope, the intensity of the emergent ray, which is given by* 

“R’ 


-R° 
| Burdenp| — |" eye(r) de} nyo) ds, 
/re=R Jer J 
becomes, in the approximation R’— R’, 


“RY “R’ 
BYR) | exp { = | K,p(r’) as «,p(r) ds 

J peR’ =r J 

= B(R){I— exp [—7(v)/1]}. (3) 
By a consistent use of this limiting process we can take account of directional 
dependence in a manner more adequate to our purposes than would be provided 
by the traditional approximation of an outward and an inward beam. The 
approximation of a thin shell need not be given a literal interpretation: R can be 
considered as an effective ‘‘ mean radius” for even an extended envelope. 

Apart from its control over directional variation of the specific intensity, the 
general equation of radiative transfer incorporates another important relationship 
which we must be careful to retain even though we dispense with the transfer 
equation itself; namely, the energy conservation law. Its application will 
provide us with the necessary equations for calculating our contours. 


3. Profiles based on the scattering mechanism 


3.1. The integral equations of the problem.—Consider a perfectly absorbing 
body of radius Ry (see Fig. 1 a) radiating in all frequencies with specific intensity F. 
Let it be erclosed by a ‘‘thin shell’’ (in the sense above) of radius R, in which 
matter is moving radially outwards with velocity v; and let the optical depth 
for a ray normal to the shell be 7, as in (1), v being the frequency as seen by an 
observer moving along with the matter in the shell. If J(v,) be the intensity 
of a ray striking the shell at an angle cos“! 1 to the normal, with apparent frequency 
v, then by use of (2) Jexp[—7(v)/] is the fraction of the ray transmitted without 
suffering any scattering (Fig. 15); hence 


I(v, w){I —exp[—7(v)/u]} (4) 





* cf., e.g., S. Chandrasekhar, Radiative Transfer, p. 9, equation (50), Oxford, 1450. 
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is the amount of the ray which is scattered. The total flux of radiation thus lost 
in scattering by all rays which cross from the inner to the outer face of the shell 
is therefore 


| Hem) —expl—1(v)/q]}ye deo = 20 J Ho.n){x exp [-2(0) aD) dp (5) 


per unit area, where the + sign indicates that the integration is to be carried over 
the half sphere subtended by unit surface element on the outer face of the shell. 


cos > le" 


Let B(v) be the emissivity of matter in the shell. (Again v is the apparent 
frequency.) ‘Then, according to (3), the intensity of a ray emerging with direction 
cos! yw is* 

s fe 1s 
B(v){1—exp [ —7(v)/]}, (6) 
and in consequence the energy flux out of the shell is 


| Be) -expl=7(>)/[H 1} [el deo. (7) 


The integration here is to be carried over the whole unit sphere of solid angle 
{ dw = 47, since the emission as a result of scattering will emerge from both the 
inner and outer faces of the shell. Formula (7) can be integrated out to give 
1 
2n | B(v){1 —exp[—7(v)/ |e] ]} |u| du =27B(v)[1 — 2E,(7,)]. (8) 
J-1 
The exponential integral of order n is defined as 


. , dt 
E(x) I, exp(—tx) =. (9) 





* This expression is the limb “‘ darkening ”’ law for a thin shell. Note that in point of fact 
it makes the limb of the disk »=o brighter than the centre »=1 in the ratio 1 : {1—exp[—r(v)]}. 
The brightness of the limb is of course one of the striking features in photographs of planetary 
nebulae where the condition of a thin shell is well approximated. 
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If there are no sources or sinks of radiation within the shell, then the energy 
lost from the star’s radiation field in unit area of the shell by scattering (5) must 
equal* the amount emitted by the shell (8); thus: 


anB(v)[X— 2Ey(7,)] =27| I(v,u){1—exp [—(v) /u] dy. (10) 


Now the radiation striking the inner surface of the shell can, in this system, 
come only from other portions of the shell or from the central star itself. A ray 
judged as having frequency v and angle cos“! to the normal by an observer 
moving with the scattering material (that is to say, radially outwards with 
velocity v) would be judged by an observer at rest at the centre of the star as 
making an angle cos! with the normal to the shell, and as having frequency v’. 
Since v’ —v is due to the Doppler shift, 
vy =v[I+pv/c], (11) 

in which c is the velocity of light. Similarly an observer stationed at the point 
on the other half of the shell where the continuation of the same ray strikes the 
shell again would judge the frequency as being v”, if he participates in the local 
motion, wheret 

v” =v'[I1+po/e] =r[1 + 2p0/c]. (12) 
See Fig. 1c. 

Rays making an angle with the normal at the point of incidence whose cosine 
falls between 1 and [1—(R,/R,)?]"? =, originate at the surface of the central 
star (see Fig. 1d). The latter has been assumed to be radiating with intensity F 
in all frequencies and directions (Fig. Ie), so 

U(v,p)=F, oy <p. (13) 
For all other directions of incidence the rays originate on the shell, the intensity 
being given by (6) and the frequency by (12), so 

U(v, w) = B(v"){1 — exp [—7(v")/]} 

= B(v[1 + 2uv/c]){1 —exp [ —(1/u)r(o[ 1+ 2u0/c])]}, O<p <p. (14) 
If (13), (14) be substituted into (10) there results an integral equation 


1 
BU)[x—2Ey(r,)]= | F{x—exp[—r(»)/p] ed 


+ [ B(v[1 + 2nv/c]){1 — exp [ — (1/)7(v[1 + 2u0/c]})]} 


x {1—exp[—7(v)/p]}udu 
for the determinatien of the emissivity B(v). 

Once B(v) is determined, the combined flux leaving the system of star plus 
envelope in any frequency is known. Apart from an amplitude factor this 
combined flux is the spectrum contour. ‘The flux is the sum of contributions 
from: 

(i) the star as seen through the forward face of the envelope, 

(ii) the back half of the shell as seen through the forward half, and 

(iii) the forward half, or what is usually called the disk. 

* In equation (10) we have equated the energy in each frequency separately. That is to say 
we assume coherence. Since the frequency referred to as v is always that judged by the scatterer, 


this represents no departure from standard procedure. However, see references (15, 16). 
+ Second-order terms in (v/c) are of course ignored throughout. 
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After escaping from the outer surface of the envelope the radiation suffers no 
redistribution in frequency or direction and the flux can therefore be evaluated 
at any surface enclosing the system of star and envelope. It is most convenient 
to take for this purpose the outer surface of the shell and judge frequencies in a 
frame of reference at rest with respect to the star; from there outward the inverse- 
square dilution law holds. Quantitatively (i) is given by 


-1 
2m | Fexp {—(1/u)r(v[1 —pv/c]) }u dp, (16) 


My 
where, in the integrand, F is the intensity of the ray as it leaves the surface of 
the star and exp {(—1/)r(v[I1—jev/c])} is the fraction of this which traverses 
the shell without scattering. The optical path length (1/)7(v[1—pv/c]) is that 
which the scatterer would apply to a ray which originated at the stellar surface 
with frequency v. Similarly (11) is given by 

“wy 


an | B(v[1 + pv/c}){ 1 —exp[ —(1/4)7(v[ 1 + pv/c])]} 
“0 
x exp { —(1/p)7(v[1 —pv/c]) }u dp. (17) 
The first two factors in the integrand are together the intensity (6) of the ray 
leaving the shell with frequency »[1+v/c] as seen there, but with frequency v 
as seen in the chosen stationary frame. ‘The factor exp {(—1/u)r(v[ 1 —p:v/c])} 
has the same significance as in the expression for (i). Lastly, (iii) is given by 


“1 
2n| B(v[1 —pv/c]){1 —exp[ —(1/n)r(v[1 — pv/c}) Ju dp, (18) 
Jo 


following the same reasoning as for (ii); but here no factor for absorption is 
required since the radiation is already outside the shell. If (16), (17) and (18) 
be summed one obtains the expression for the spectrum contour, apart from a 
multiplicative constant, 


rl 
F(v)=F | exp { —(1/u)r(v[ 1 — pv/c]) }u dy 


+ |" B(v[1 + pv/c]) {1 — exp [ —(1/4)r(v[ 1 + pv/ce])]} 


x exp { —(1/u)r(v[ 1 —pv,c]}) ju du 


“y 
+ | B(v[ 1 — pv/c}){ 1 —exp[ —(1/n)7r(r[ 1 —pev/e})]} dp. (19) 
/0 


3.2. The form of solution of the integral equations.—Equations (15) and (19), 
which are now the fundamental equations of the problem, take on a neater form 
if the frequency of each point in the spectrum is measured as displacement from 
the normal frequency of the resonance line in units of the Doppler shift produced 
by a velocity v equal to the velocity of expansion of the envelope. Set 


or v=v(I+yv/c), (20) 


wherein vy is the frequency of the ‘‘normal”’ line centre. ‘Then 7(v), B(v), A(v) 
become functions 7*(y), B*(y), A*(y). In the expression 7(v[I—pv/c]), for 
example, 
v[I —pv/c] = ve(I + yu/c)(1 —pv/c) = vl 1+(y—p)2e/c]; (21) 
hence 
7(v[ I —pv/c]) =7*(y—p). 
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No confusion will result if the asterisks are dropped, it being understood that 
henceforth the unit of frequency is always (20), unless the contrary be specifically 
indicated. ‘The normal line centre corresponds to y=0. Equations (15), 
(19) become respectively 


1 
BUy)[x— Es(o(y) =F | (1 —expl 29+ 21) ]}u dy 


‘My 


+ a By + 2){1 —exp[—7(y + 2pu)/p]}{1 —exp[—7(y)/H]}u dp, (22) 


1 
F(y)=F I exp[—7(y—p)/m)u du 


+ B(y + u){1 —exp[—7(y +4)/u]}}lexp —7(y—#)/e le dy 


+{ By-H){x-expl—r(9— wind. (23) 


Given the parameter y, and the function 7(y), equation (22) can be solved for 
B(y) and the result used in (23) to determine A(y). Clearly the method of 
solution will have to be numerical. 

Its appearance to the contrary notwithstanding, (22) is a very simple equation. 
For any value of y the right-hand side of the equation involves only values of 
B(y +2) in which y+2u>y, since the integrand contains the factor » while B 
is always finite. (Note the absence of the equality sign.) ‘Thus the emissivity 
at any frequency is explicitly given by the emissivity in frequencies further to 
the violet.* ‘This is a consequence of the fact that each scattering particle in 
the envelope is receding from all other particles in the system.t The deter- 
mination of B(y), for large positive y first, is therefore indicated. 

Assume for y large that r(y)—>0. (The optical depth falls off to zero outside 
the resonance line.) ‘Then on the left-hand side of (22) the factor [1 — 2E,(7(¥))] 
is of order 27. ‘The first integral on the right-hand side is of order F(1—,.,)r 
whereas the second integral is of order r?.__ For y large, therefore, 

By) > }F(1— 4): (24) 
Thus in calculating B as a function of y one need only start the integration of (22) 
sufficiently far in the violet for the last integral on the right to be negligible in 
view of the accuracy desired. (24) can be used in the integral to estimate its 
magnitude. ‘Then the right-hand side of (22) can be evaluated at successive 
values y, of y as one moves towards the red end of the line. ‘The spacing of the 
points y, at which B(y,) is determined can be chosen to maintain the ease and 
accuracy of the numerical integrations involved in the succeeding evaluations. In 
equation (23) the last two integrals are each of order 7 whereas the first one on the 


right-hand side is of order F[’ edu. For y large, therefore, 

Fy) > 4F(1—p,°) (25) 
or 

F (y)> (1+ py) B(y). (26) 





* This necessity of proceeding in one direction is a recurrent feature in problems of a moving 
atmosphere. For example, in reference (13), the integration must be started at the violet end. 
If it is begun at the red end the solution remains indeterminate. 

+ If the envelope were contracting the red frequencies would determine the violet ones. The 
velocity v in (20) would then be negative. 
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In view of the arbitrariness of the multiplicative factor mentioned in the para- 
graph following equation (15) above, it will be convenient to arrange that 

$F(1—p,*) =1, (27) 
so that according to (25) the profile of the line rejoins the continuous spectrum 
with the value unity. 


It is necessary now to assign the function r(y). According to (1), when 
R’ > R’=R, 
+R" 
T.>k, | p(r) dr; 
- 
hence 
7(y) =(constant) «(y). (28) 
When a resonance line is broadened by thermal agitation or turbulence, the 
absorption (read: scattering) coefficient is given by (17), 
v—vy c\2 
x, =(constant). exp; —{ —— -4] >, (29) 
Yo wv 
where the most probable velocity 
v* = \/(v" + U7"), 
of an atom in the Maxwell distribution, is the resultant of the thermal and turbulent 
velocities. In the frequency units (20) here employed, (29) becomes 
«(y) =(constant) exp ( — f?y?) (30) 
and consequently (28) gives 
7(y) =7(0) exp (— B®y*), (31) 
where the new parameter 
B=v/v* (32) 
is the ratio of the directed velocity (radially outward) to the most probable random 
velocity (relative to the local rest frame) of the atoms comprising the material. 
For increasing velocities of expansion the turbulent velocity must be expected to 
increase too, so that the parameter is unlikely to take on extreme values.* 

3-3. The profiles for a range of values v/v*.—The parameter 4, which appears 
in the equations (22), (23) measures the ratio of the radius Ry of the nuclear star 
to the effective radius R, of the envelope (Fig. 1d): 

Hy = V[1—(Ro/R,)?]. (33) 
In this section we let ».,->1, which gives the case of a point star, postponing to 
the next section consideration of the changes produced by p, “1. 

When p.,—>1, subject to condition (27), the first integrals on the right of (22), 

(23) become respectively 


F{ (2-expl—1(y)/H]yedu> (1 —expl—r(y)]}F | de> (1—expl—n(y)} 
(34) 


7] 
F | exp[—7(y—p)/u]udu—>exp[—r(y—1)]F | pedu—>exp[—r(y—1)]. 


~§ 


(35) 





* The profiles obtained below indicate values of 8 of order 3, when the velocity of expansion 
is ~ 200 km/sec. Not all the width in v* is necessarily due to turbulent velocity in our method 
of thin shells. ‘The assumed velocity of expansion v of the envelope is only an average: any 
gradient in v would be mixed up with the turbulent width. 
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Consequently the equations to be solved become, for y, = 1, 
B(y)[1 — 2E,(7(y))] = {1 —exp[ -—7(y)]} 


1 
+ | Bly + 2y){1—exp[—r(y + 2y)/a]}( exp — 709) uD a di (36) 
“1 
F (y) =exp[—7(y—1)] +) Boru) —exp[—7(y+p)/u]}exp[—7(y—p)/p] udu 


<. 
+) Boy) —exp[—7(y—#)/m]} udp. (37) 
The limiting values (24), (25) are given by 


By)>3, F(y)>1, (38) 


for y large and positive. ‘The equations (36), (37), (38) have been numerically 
integrated, after setting 7(0) =1 in (31), for three values of v/v*: B =1, 1/10, 10. 
The results are shown in Tables I, II, III, and A(y) is plotted for these three 
cases in Fig. 2. ‘The scales of the drawings have been chosen to keep the turbulent 
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width the same throughout (constant v*); whereas the unit of frequency used 
in the tables is such that for each profile the expansion velocity gives unit Doppler 
displacement (v=1). Equations (36), (37) may be written symbolically as 

B(y) = B,(y) + BAY), 

F (y)=F fy) +F (vy) +F2(y). 
Here B,(y) is the contribution to the emissivity which arises from initial or 
primary scattering of the star’s emission, and B,(y) is that from scattering of the 


diffuse or secondary radiation (that is, rays coming from other parts of the envelope). 
F (y), F;(y), F2(y) are the contributions to the resultant profile which were listed 
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respectively as (i), (ii), (iii) in the paragraph (q.v.) following equation (15). 
These functions have also been tabulated in order to show the relative importance 
at each frequency of the various components. 
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From Fig. 2 it will be noticed that the profiles are characterized by an absorption 
component to the violet of the line centre. Since we arranged for conservation 
of energy absorbed and re-emitted, and since we gave the star no cross-section 
for capture of radiation scattered by the shell (point star), it was assured that 
there would therefore also be emission components of total equivalent width 
equal to the equivalent width of the absorption component. (‘lhis provides a 
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convenient check of the accuracy of the profiles.) What was perhaps unexpected 
was that for values of 8 not much greater than unity the re-emission components— 
one from the forward face of the envelope, F,(y) above, the other from the rear 
of the envelope F ,(y)—would begin to separate so as to give an apparent “ reversal” 
which falls almost exactly at the line centre. That is to say, as a result of the 
presence of the scattering shell a portion of the flux in the neighbourhood of a 
resonance line reaches the observer after having been redirected by the shell. 
The energy built up in this scattering envelope has a maximum, which for the 





ee ee ee Se ee 


























i 1 l j 
4 2 fe) 2 a 


E Q U E N 





Fic. 2.—Dependence of the profiles on the velocity ratio B=v/v*. 
vu ts the velocity of expansion and v* is the “‘ most probable random’? velocity. 


forward half of the shell falls at a frequency to the violet and for the rear half of 
the shell falls at a frequency to the red of the normal frequency, as judged by a 
stationary observer. 

3.4. Effect of occultation on the calculated profiles—To investigate the effect 
on the profiles of a finite radius of the “ perfectly absorbing” star—hence the 
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ability to occult part of the rear half of the envelope—equations (22), (23), (27) 
were integrated for three additional values of 4, namely 0-8, 0-6 and 0-0, with 
B=~+/10. The results are shown in Tables IV, V, VI, which should be compared 


with the corresponding one for », =1, B= 1/10 in Table II. 
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R, is the radius of the star; R, is the mean radius of the envelope. 














Fic. 3.—Effect of occultation. 
For #4, =0 equations (22), (23), (27) simplify greatly: (22) becomes 
B(y)[1 —2E,(7(y))] = Fl {r—exp[—7(y)/p] }u du 
= 4F(1-—2E,(r7(y))] 
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and (27) gives 4F=1; hence B(y)=1. When these values for F, B are 
substituted into (23) the contour follows: 


Fy) =F expl—ay—w)inlndu+ | B(y—n){x-expl—a(y—H)/uD ede 


=2} expl—7(y—m)/ulu du + | {r-exp [-r(y—p)/e]}e de 


rl 
=$+ | expl —7(y—p)/H]e dp. 


Since the integral in the last line above is always < }, the net effect of the two terms 
is an absorption line with no accompanying emission component. 

The four curves for A(y) are shown in Fig. 3. The diminishing effect on 
the emission components, it will be seen, does not become noticeable until the 
occultation is well advanced.* 


4. Profiles based on the Zanstra mechanism 


4.1. Emission from the recombining shell.—According to the Zanstra~Menzel- 
Beals theory of emission in early-type stars, the energy appearing in lines of the 
Lyman and Balmer series is supplied from radiative frequencies beyond the 
head of the series. In the region immediately adjacent to the inner boundary 
of the envelope the radiation field is so intense that almost all the atoms are 
ionized.t 

In this assemblage of ions and electrons, recombinations can take place. 
The probability p of the electron being captured directly into the ground state 
is, according to calculations by Cillié (20), only of order 4. In the other (1—/p) 
captures the electron appears on an excited level and if undisturbed rapidly 
cascades down to the ground state, emitting at each jump a line quantum. On 
our assumption of steady state (of the gross physical variables) these recombina- 
tions cannot be permitted to upset the degree of ionization, so simultaneously 
the photoelectric ionization of an equal number of neutral atoms must be 
occurring. ‘Thus, for each quantum hv, of radiation newly appearing in the line, 
a quantum hy, of energy greater than the ionization potential of the atom con- 
cerned has been lost to the radiation field in the ultra-violet. 

If ~S, be the total flux of radiation emitted from the surface of the star in 
frequencies beyond the series limit, and 7S,exp(—7,) be the flux in these same 
frequencies which escapes through the envelope (thereby defining 7,), the energy 
lost by conversion in the envelope is represented by 


7S[1—exp(—7,)]. (39) 
Since roughly one quantum hy, appears in the line for each quantum hr, lost in 
the ultra-violet continuum, the energy converted by the recombining region 
contributes 
7S,( vo v)I I— exp ( = 7,)] (40) 
to the flux in the immediate neighbourhood of vy. The precise distribution of 
this flux among the frequencies adjacent to vg will be determined by various 
* For 4, =0°6 the radius of the nuclear star is already 0-8 of the radius of the shell. The last 
profile of the sequence, the one with 4,;=o in which occultation is complete, is the line produced 
in a plane atmosphere. This is the same profile arrived at by O. C. Wilson (18), and again by 


A. B. Underhill (19), whose discussions are restricted to moving plane atmospheres. 
+ Upwards of all but one in 10°, in representative cases. 
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factors which ultimately determine the shape of each component (steepness of 
sides, flatness of top, etc.). We shall assume that apart from a uniform expansion 
of the envelope this distribution is mainly determined by the random motion of 
the atoms doing the emitting, so that the spectral variation of this additional 
flux can be written as 


, 
Ty 


7F'(v)=S, > rz, 1—exp(—n) (41) 


where 7,’ is the optical thickness of the recombining region in frequency v. 
That is to say, at any instant, the number of atoms capable of absorbing radiation 
of frequency v is a direct measure of the number that have just arrived in that 
state by emission of radiation of that frequency. ‘The flux in any frequency is 
also expressible in terms of the specific intensity emerging from unit area on the 
shell : 


7F'(v)={I'(v,n) |u| do, (42) 
where the integration is taken over both faces of the shell; that is, over the solid 
angle fdw=47. ‘To obtain the variation of /'(v,) with direction we again make 
the approximation that the recombining region is a thin shell of radius R,. ‘The 
dependence of /’(v, 1) on » is then given by a formula analogous to (3). Making 
use of (41), (42) we find that its explicit expression is 


Wy ai Se ¥o Te I—exp[—7(v)/[ul], 
a 2 »v, fr,’ dv 1—2E,(7,') ae (43) 





For the early-type stars no direct estimate of the flux outside the visible 
region can be made. The temperatures of these stars is obtained (5, 9) by 
assuming that all the ultra-violet radiation is converted to line radiation. In 
the notation adopted here, this implies that 


T=. (44) 
Further, in supposing that the actual physical state of the envelope can be 
approximated by two thin shells of different radii, one of which is characterized 
by the Zanstra emission process and the other by a process of scattering by a 
moving medium, we have inferentially supposed that the first is also characterized 
by a very high degree of ionization*, while the second is composed of a sufficiently 
high proportion of neutral atoms to maintain resonance scattering. In view of 
this separation of function we may, without stretching the approximation, 
assume that the recombining shell is so far ionized as to be transparent to radiation 
in the frequency of the line. Hence 


T, 0. (45) 


Nevertheless 7,’/{7,’ dv can still be finite. (As, for example, with a coefficient 
of the form (31) when 7(0)->0.) Lett 


in dv i en (49) 





* That this is compatible with equation (44) is admittedly open to argument. 

t In this fashion we account for the possible non-coherence of the emitted radiation, as discussed 
by L. Spitzer (15). When one uses, as we shall, a coefficient of the form (31) to determine (46) the 
emission is distributed in frequency in the fashion postulated for planetary nebulae by H. Zanstra 
(16), Section 3. 
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With use of (44), (45), (46) the expression for the emergent intensity is considerably 
simplified. (43) becomes* 


; S.vea(v 
I'(v, p) = a . 


4.2. The equations for the two-shell model.—lf there is a recombining shell 
(in the sense of the preceding section) of radius R, lying between the star and 
the scattering shell, the right-hand side of equation (10) will involve not only 
the incident radiation given by (13) and (14) but also a term which measures 
the flux of energy scattered by unit area of the outer shell from rays which 
originated at the inner one. 

The rays which originate at the forward face of the recombining shell R,, 
and reach the scattering point, lie within a pencil whose extremes are the ray 
which is perpendicular to and the one which just grazes the R, shell. The 
cosines of the angles these rays make with the normal to the scattering shell R, 
fall in the range I >>. ‘Those rays from the R, shell which by-pass the 
occulting star R,, and traverse the forward face of the R, shell to reach the 


(47) 


2 cos" eg 


(2) (b) 


Fic. 4. 


scattering point, lie within a pencil bounded by the ray which just grazes the 
star and the one which just grazes the R, shell. ‘The corresponding cosines 
fall in the range p, > Spy. (See Fig. 4a.) 
My = [1 —(Ro/R,)*]*?, Be =[1—(Ry/R,)?}*”, 
where Ry, R,, R, are respectively the radii of the star, the recombining shell, 
and the scattering shell. A ray, which strikes the R, shell at the angle cos"! , 
leaves the R, shell at the angle cos !,.’, where 
(1 —p*) =(1—p’?)(1 — 1"). (48) 
Matter in the R, shell may have a radial velocity, say v’, of which one must 
take account. A ray which appears to the scatterer at R, to have frequency v 
and direction » seemed to the emitter at R, to have frequency v’ and direction 
pw’, where (Fig. 4) 
v[I+pvr/cl, (49) 
v[I +pv/c], (50) 





* Note that the singularity at «=o, introduced by letting 7,’ -» 0 would make the radiation 
density $f ldp non-convergent if calculated on the above formula; that is to say, if the limit were 
taken first and the integration performed later. However, in our application we need only the 
flux fz p dw, which remains well behaved, since »/ in (47) is finite for all p. 

Il 
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according as the ray is, respectively, from the approaching or the receding disk. 
If (49), (50) be written in the Doppler units (20), then a ray which appears at 
the R, shell to have frequency y will at the R, shell have had frequency y’, where 


V4 =Yte— ype’ =ytm(y')— yp’, 
m1) = [w2(I — 4g") + 4" )¥* (51) 
y =yt ety’ =y+m(u')+yv’, 
according as the ray is from the advancing or the receding disk. Here 


y=u'/v (52) 
is the ratio of the two expansion velocities. 
Proceeding as in (4), (10), (14), we find that the flux in frequency y which, 
having originated at R,, is scattered upon reaching R, is (omitting the 27) 


|. I'(y,', 4’ {1 —exp[—7(y)/m)}u du+ { : I'(y_’,w'){1—exp[—r(y)/u]}udp. (53) 


If the variable of integration be changed according to formula (48), then (53) 
becomes 


myu’)= 1 
(r—mat)| [7 Oem) —expl—ay)im(u’ De de 


mj") = 1, 
+ s “ae I'(y_’,p’){1— exp | —a(y)im(a I} de’ (54) 


The intensity of the recombination emission has been written down in (47). 
In the present variables 
Ya a! » Yo OY 
I'(y',w')=38, — mae (55) 
‘Thus (54) is given by 
2 S.vo9 f . ’ , , 
(1 — pe io ) o(y,’){1 —exp | —r(y)/m(u’)]} du 
r “ 0 
rl(uy* — ay") 1—p,*)}'" 


+ | : o(y’){1—exp | —x(y)/m(u’)1} du’ (56) 


If y’ is substituted from (51), and if it be agreed that the positive square root 
be always the one taken, then the above expression for the amount of new radiation 
scattered can be written as the single integral 

“em “1 
(1 —1,)| o(y + m(n') — yp’){1—exp[ —7(y)/m(u’)]} du’. (57) 
4Y. ~ —[y*—42)/—n,*)}" 

‘This is the term to be added to the right-hand side of (15) in order to obtain the 
integral equation of the emissivity of the scattering siiell in the two-shell model. 

Proceeding in the same manner one finds, for the term to be added to the 

right-hand side of (19) to represent the flux of the emission from the recombining 
envelope which penetrates the outer shell without scattering, the expression 
2 S.vo 


(I—H2 o(y—yp')exp{—ztly—m(u’)]/m(u’)} du’. — (58) 


rl 
4¥¢ | — (Coa? — na" — 9)? 
Comparing (57) and (58) one can understand the significance of the component 
factors, of which one may have lost sight during the changes required to simplify 
them. ‘The second factor in the integrand of (58) is the fraction of the original 
ray which penetrates the outer shell when it arrives there at an angle whose 
cosine is (cf. (48)) « =m(yz’) and with a frequency which is y when judged in the 
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rest-frame of the star (as was agreed upon in the paragraph preceding equation (16)), 
but is y—m(yu') when judged by the scatterer. On the other hand, the first 
factor gives the strength of the emission along that ray in the same frequency y, 
but corrected now for the fact that the recombining shell, at which the ray 
originates, has velocity y in units of the radial velocity of the outer shell and for 
the fact that the cosine of the angle made at this inner shell is »’. In contrast, 
the first factor in the integrand of (57) has its frequency determined by the fact that 
the frame of reference for judging frequencies as y for use in equation (15) is 
taken to be that of the material at the scattering point (cf. the paragraph preceding 
equation (I1)) and must therefore be adjusted by a term m/(j’) to reduce it to the 
rest-frame of the star and by a term yp’ to bring it to the rest-frame of the point 
and direction of emission. ‘The second factor of the integrand of (57) is then 
the fraction of the ray momentarily absorbed by the scattering material when the 
frequency it sees is y and the cosine of the angle made with the normal is x = m(j’). 

The addition of terms (57), (58) to our equations introduces a number of new 
parameters. ‘lhe variation of these parameters so greatly increases the variety 
of profiles obtainable with our methods that an exhaustive investigation of all 
resultant types would run to exceedingly great length. We necessarily contine 
ourselves, therefore, to deriving the profiles for a few illustrative cases. 

Since the effect of occultation on any profile is readily predictable, and has 
already been quantitatively investigated for the case of a single scattering shell, 
we shall restrict the calculations to point stars. ‘Then the equations for the two- 
shell model are given by adding the terms (57) and (58), for «, =1, to equations 
(36) and (37) respectively : 


S.vq (' 
B(y)| 1 — 2E,(7(y))| = (1 — 22") = | a(v + mu) — yp) 1 —exp[ —r(y)/ m(ye) |} dye 
Ve - l 


1 
{1 —exp[—r(y)]}4 | By + 24){1—exp| —7(y + 2u)/u]} 


« {1—exp[ —7(y)/p) ue dp, 


| 


” a aM 
F(yj=(1- re) | a(y— yy.) exp | — rly — m2) |/m(2)} du 
AV l 


@ ¢ 


“1 
+exp|—7(y—1)]+ | Biy — pw) {1 —exp[—r(y—p)/p) pu du 


rl 
+ | Bly + p)iI— exp| (y+) elfexpl—7(y—p)/p]y du. (60) 


4.3. Profiles for the limiting case of small recombining shell.—When the re- 
combining shell is of radius small compared to the effective radius of the scattering 
shell (in this point star model) then (Pig. 4@) po =|1—(R,/R,)?|!2~1. Lf one 
lets jz. >I in (58), (60), but in such a manner that S,(1 —j,*)— S,’ remains finite, 
it is easily verified that these equations become identical with (36), (37) except 
for the first term on each right-hand side, where the factor 

So. 


| 
av . 
I+ — | o(y—ye)dp -1+X(y), (61) 
4% J-1 — ; 
representing the flux in frequency y of the star plus recombining envelope, 


replaces the unit flux (27) from the star alone.* 


* When the combined flux (61) does not vary too rapidly with y the resultant profile F(y) 
will be reasonably well given by superimposing the corresponding contour of Section 3.3, F*(y> 
say, on to the function 1+ 2(y); thus F(y)=F *(y) .[1+2(y)]. 

11* 
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The integration of the equations has been carried out for several choices of 
expansion velocity and turbulence in the recombining shell. The profile shown 
in Fig. 5 together with the values given in Table VII is the computed result for 
the case 

y=v/v=2,  p'=v'/'*=3, B= v/v®=3-16; (62) 
that is, when the recombining shell is expanding twice as fast as the scattering 
shell (v’ =2v), and when the turbulent velocity in each shell is approximately a 
third of the associated expansion velocity (v=3°16v*, v’ =3v'*). The optical 
depth in the scattering envelope is again taken to be of the form of (31) 


7(y) =} exp(—*y?*), (63) 


and a similar choice is made for the recombining envelope, in accord with (46). 
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Fic. 5.—The profile of Section 4.3. 
The star radiates the continuous spectrum 1. & is the flux due te recombination. 
As a result of scattering, a distant observer sees the line contour # . 


The flat-topped, symmetrical contour (fine line in the figure, 1+ 2 in the 
calculation) underlying this complex profile of P Cygni type is the shape the line 
would take if there were no scattering envelope. For then the emission from the 
recombining envelope (and from the continuous spectrum) would be the only 
visible feature. Flat tops indicate emission from an expanding medium of 
spherical shape. ‘The effect of the scattering layer is to produce an absorption 
component, displaced to the violet in proportion to the velocity in the line of 
sight of the scattering atoms. Since the expansion velocity of the scattering 
material was taken, for this illustrated case, to be but half that of the recombining 
matter, the absorption attacks the emission profile where the latter is still quite 
intense. For radiative balance the energy thus absorbed must be re-emitted. 
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The energy reappears in the peaks on either side of the line centre, producing an 
apparent “central reversal’’. 

4.4. Profile when the shells coincide.—When the recombining and scattering 
shells have the same effective radius, ..==[1—(R,/R,)*]"*. If the expansion 
velocity of matter in each shell is also the same, then y=1. This is the particular 
case we should have had, had we not made the distinction in Section 4.1 between 
recombining and scattering shells. ‘The equations resulting from (59), (60) 
have been integrated with 

Seo o(y) = 107(y) = 10 exp(— 1oy?). (64) 


Al 
2v, 


‘The results are shown in Table VIII and the profile in Fig. 6a. It will be noted 
that the rather large recombination flux adopted has almost entirely filled in the 
violet absorption border. However, the most striking feature of the profile is 
again the fact that the emissivity built up in the outer shell by the scattering 
process separates into two emission components leaving a deep “ reversal’’ 
near the line centre. ‘This may be seen from Fig. 66 where the contributions 
to the profile of starlight F,’, envelope emission F,", and diffused light from 
the two faces of the scattering shell 4, +A, are separately shown. 
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‘The absorption line, produced in the A,’ curve at a displacement from the 
line centre corresponding to the velocity of the scattering shell, is the usual 
reversal of parallel light in passing through an absorbing medium. ‘The attenu- 
ation of the symmetrical flat-topped contour = into the highly asymmetric Fy" 
is, however, quite striking, and shows to what extent the scattering-recombining 
shell is transparent to radiation from its opposite face. Frequencies to the red 
of y =0-5 escape almost completely. The diffuse emission as a result of scattering 
is, as we have now learned to expect, concentrated into two peaks with a minimum 
near the line centre. 
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lic. 6a. The pre file of Section 1.4. 
The star radiates the continuous spectrum 1, = is the flux due to recombinatu 


As a result of scattering, a distant observer secs the line contour 7. 




















Fic. 6 b.—The components of profile 6 a. 
The star, seen through the scattering envelope, provides the flux F 4’. The flux due to recombination 
gives the component F 9”, when seen through the scattering envelope. Re-radiation by the envelope 
adds tne contour (F ,+ Fs). 


The equivalent widths of the A and 1+ profiles are—within the limits 
of computational accuracy—equal. 

4.5. Profile for an intermediate case.—-Vinally, in order to verify that no 
untoward effects intervene for choices of R, intermediate to the limiting cases 
R,=R, and R,=R, considered in Sections 4.3, 4.4, above, the integral 
equations ((59), (60)) were solved for \2R,=R,, v v=}: 


I S.y, 
Hy =1, F(I—"y)=1, pe= o(y) =? (1 — ps) = 107(y) = 10 exp (— 10y*). (65) 
\ 2y, ; 
The values found for the various components are shown in Table [X and the 
profile in Fig. 7. The profile is of a familiar type. 
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Note that in setting o(y)«x7(y), we have taken v*=v'* and hence have 
adopted as turbulence ratio for the recombining shell 


, 


" 
7 — —_ 
P= a= 


a value which, when y = }, gives half the value in the scattering shell. Emission 
after scattering from the rear and forward faces is not sufficiently displaced in 
frequency. As a consequence the profile fails to show central reversal. 


v 
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It is interesting to contrast Tables VIII, IX. But for the presence of a 
scattering envelope the contours of Figs. 6 and 7 would have been similar: an 
emission profile 1 + & symmetric about the line centre. In the model described 
by Table IX, however, not as much of the energy finds its way into emissivity 
B,”, B, of the scattering shell, because the latter is farther removed in position and 
velocity from the recombining shell than was the case for Table VIII. 

5. Conclusions 

5.1. The contours and the physical model.—In appraising the extent to which 
our assumptions about the principal physical features of a P Cygni type envelope 
are verified by the profiles, we can do no more at this stage than compare the 
calculated contours with the P Cygni type profiles considered as a class.* 
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Fic. 7.—The profile of Section 4.5. 


The star radiates the continuous spectrum 1. & is the flux due to recombination. 


As a result of scattering, a distant observer sees the line contour #. 


Relative to this more limited objective, the agreement between the calculated 
profiles and the observed ones is good. ‘To illustrate this we reproduce a number 
of figures (Nos. 8-11) kindly provided by C. S. Beals from his publication “ The 
Spectra of the PCygni Stars” (§). ‘The numbers appearing on the left of 
Figs. 8 and 9 give the type of profile according to Dr Beals’ classification. For 
our purposes these are then representative observed profiles. Except for types 
II, VI, whose broad wings are indicative of Stark effect (which we have not 
attempted to take into account), each of these observed profiles has its mate among 


* A quantitative analysis of the observed profiles of a particular star in terms of the parameters 
used in our calculations would require an investigation all its own. However, such an investigation 
appears entirely feasible. 
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the calculated profiles of Figs. 2, 3, 5, 6, 7 or among the profiles that may readily 
be interpolated from them. 

Special attention should perhaps be drawn to one or two specific cases; 
for example, the remarkable parallelism of the Balmer series in HD 190603, 
Fig. 10, to the sequence of profiles showing the effects of occultation in Fig. 3. 
Also noteworthy is the similarity of the calculated profile in Fig. 6 with the 
observed hydrogen emission in Be spectra. Fig. 11 shows the Ha line in 
HD 218393. Despite the complexity of the latter profile it is readily seen that 
all its features can be accounted for by interpolating from the profile of Fig. 6 
towards the profile of Fig. 7. 























Fic. 8.—Typical profiles of P Cygni-hke stars. (After Beals.) 


‘Thus, as a class, the P Cygni stars come within the range of our model. ‘This 
agreement with the observed profiles permits us to recast the assumptions about 


the model into assertions about these early-type stars. 
(1) They are enclosed by an expanding atmosphere, whose outer regions are 
chiefly engaged in scattering. This is borne out by the orientation of absorption 
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and emission components in the spectral profiles: absorption minimum on the 
violet, emission maximum over the central-red frequencies. (When the latter 
is double, then the included minimum falls at the natural frequency of the line.) 

(2) Since occultation is not a noticeable feature of the profiles the envelopes 
extend to at least 25 per cent in excess of the radius of the star proper. 

(3) In addition, because in many of the profiles the equivalent width of 
emission exceeds that of absorption, the inner regions of the envelope must be 
emitting line radiation as a result of photoelectric ionization and recombination, 


-—---- 





























Fic. 9.—Other profiles of P Cygm stars. (After Beals.) 


According to the analysis advanced in this paper, the oft-occurring “ central 
reversal’’ is not primarily the result of absorption by matter with zero velocity 
in the line of sight. It is rather a consequence* of the fact that the average 
velocity of expansion is greater than the resultant of all non-uniform velocities 
in the envelope. ‘These latter include the thermal velocities, turbulent velocity 
and any gradient or differential velocity present in the expansion. “Central 
reversals’’ arise in the computed profiles when the velocity of expansion exceeds 
the random velocities by a factor of order \/10. ‘This value provides a useful 
estimate of the velocity ratio in the envelopes of representative early-type stars. 
It is therefore noteworthy that the value thus found from the line contours agrees 
quite well with the value 2,3 found for the same ratio by Chandrasekhar (14), 
who obtained his estimate from consideration of the mechanical stability under 
radiation pressure of an expanding planetary nebula. 


* The reader is referred again to Fig. 6 6. 
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If the expansion in P Cygni-type envelopes is primarily produced by Lyman-« 
pressure on the hydrogen constituent, Chandrasekhar’s calculation applies 
equally well here.* Consequently, the prevalence among PCygni stars of 
profiles which are best fitted by a velocity of expansion at least three times the non- 
uniform velocities is not just a fortuitous fact. If the velocity of expansion were 
to fall below this value, with respect to the mean non-uniform velocity, pressures 
in the envelope would build up until the velocity increase became sufficient to 
restore equilibrium. 


aS 








H.D. 190603 


























KM/SEC, -200 


FIG. 10, 


Although this interpretation of the velocity ratio is satisfactory, the scale 
of velocities is still to be explained. ‘The observed displacements of the absorption 
components toward the violet imply expansion velocities of several hundred 
kilometres per second. Then, according to the equilibrium criterion, the 
corresponding non-uniform velocities must be of order one hundred kilometres 
per second. Indeed, the widths of the individual components themselves suggest 
a value of this order. But at temperatures around 25000 deg. not more than 


* The ‘‘ thin shell approximation ’’ that we have used replaces, in effect, the extended envelope 
by an equivalent planetary atmosphere. 
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twenty kilometres per second of this can be thermal velocity. Of course one 
would not have expected the expansion to be simply 1/10 times the thermal 
velocity, since an ejection velocity, produced solely by radiation pressure in the 
envelope, itself bespeaks a velocity gradient of comparable magnitude. And 
any velocity gradient itself contributes to the non-uniform velocities which 
widen the frequency range of the scattering coefficient. Still, it is difficult to 
see quantitatively why the expansion could not stabilize itself at smaller values 
than those observed. One cannot, it would seem, avoid having recourse to 
“turbulence” to provide the necessary breadth to the scattering coefficient. 
(An explanation in terms of the unexplained is not, however, very satisfactory.) 





















































Fic. 11.—The Ha line in HD 218393. (After Beals.) 


A point not emphasized in the previous paragraphs is the extent to which 
the choice of scattering coefficient determines the shape of each component, as 
distinct from the juxtaposition of emission and absorption components, in a 
complex line. Reference to the calculations, or a quick glance at the computed 
profiles, shows that each component and especially the wings show strongly 
the Gaussian form assumed in (29). Comparison with the observed profiles 
is in this respect also satisfactory. If reliance could be placed on the shapes of 
the components in the observed profiles being significantly free from instrumental 
distortion, it would be possible to conclude that the departures from uniform 
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motion (expansion) in the envelope are principally of the Gaussian or random 
type. ‘To the extent that the non-uniform motions are a result of turbulence, 
or of non-coherence in the scattering (15, 16), this Gaussian form is under- 
standable. But if the principal factor in broadening the scattering coefficient 
is the velocity gradient, the shape of the components should not be Gaussian. 

The conclusion, that PCygni atmospheres have an overall velocity of 
expansion approximately three times as great as their mean non-uniform velocities 
(whatever form these take), quantitatively supports Beals’ radial ejection hypo- 
thesis. The above investigation has shown however that scattering (or re- 
emission after absorption) plays an important role in the outer regions of these 
atmospheres, where it acts to dam back the radiation from the inner regions until 
multiple scattering has altered the frequency of the radiation enough to permit 
escape. ‘This makes it unnecessary to postulate any significant decelerations of 
the motions in this region in order to explain the central reversals, although such 
decelerations may be present. 

5.2. The computations.—The correspondence of the computed profiles to 
the observed ones attests to the validity of the approximations necessarily made in 
order to derive them. Chief among these was the averaging process involved 
in treating the atmosphere or envelope as one or more thin shells. Mathematically 
this meant that the equation of transfer was given short shrift in order to concentrate 
on the boundary conditions of the problem. ‘The importance of boundary 
conditions in problems of this sort was to a certain extent foreshadowed by 
Milne’s classic treatment of the radiation problem in planetary nebulae (21).* 

In introducing the method of approximation in Section 2, the only justification 
given for the process was its convenience. It was desired to study the conse- 
quences of a radiating system in which the several parts are moving with respect 
to one another and in which the occultation or obscuration of one part of the 
system by another may be significant. For this purpose it was necessary to have 
an adequate representation of the dependence of the radiation field on direction, 
particularly as required to relate the rays reaching one component of the system 
to the rays leaving another part. ‘The traverse of the radiation field through 
the envelope is not of immediate concern, when one’s interest is solely in the 
flux leaving the system and reaching the observer. All that one requires of the 
approximation is that it reproduce as well as possible the sum of radiation entering 
and leaving the various parts of the system in any one frequency; and this, as 
has been seen, can be accomplished without explicitly using the equation of 
transfer. 

From a purely pragmatic view, then, the procedure adopted is the correct 
one for calculating profiles produced by an atmosphere in motion. ‘The functional 
difficulties associated with using the integro-differential equation of radiation 
transfer are thereby avoided. Secondly, the effect of the boundary conditions 
on particular solutions can be interpreted easily in physical terms. Lastly, the 
two integral equations, for emissivity and flux, which form the basis of these 
calculations yield very readily to numerical computation. 

There is some point now in enquiring why the “thin shell’’ provides a good 
approximation to an extended envelope for the calculation of line contours. If it 
were hypothetically possible to follow the trajectory of an individual light quantum 
as it leaves the photospheric surface, one could classify the quanta which reach 


* See also (14). 
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the observer's recording apparatus in each frequency in the following fashion: 
(i) those whose trajectory is a single straight line (no scattering), (ii) those whose 
trajectory consists of two straight segments (single scattering), (iii) those whose 
trajectory consists of three straight segments (double scattering), etc. If one 
were to compare the numbers in corresponding classes when the envelope is 
regarded as a “ thin shell’’ and when it is not, the numbers in (i) would be equal— 
by our definition of optical depth. In (ii) they can differ only as a result of 
velocity and density changing rapidly with depth in the extended envelope, 
and allowance was made for this in the “turbulent width’’ of the absorption 
coefficient. Or, in terms of the two segments of the trajectory which defines 
this class, there is no effective change introduced by variation in length of the 
segment which joins the photosphere to the scattering point. Not until class 
(iti) is reached does the ‘thin shell” description begin to depart from that of 
the extended envelope. ‘Ihe former gives an adequate accounting only of those 
quanta which traverse a considerable distance between scattering points and 
hence suffer a sizeable Doppler shift. Reference to Tables I, Il, LI shows 
however that singly scattered radiation B, represents the greater part of the total 
scattered radiation density B. Small changes in the B, columns will not notice- 
ably change the contours. 


It would be impossible to conclude this paper* without recording the author’s 
indebtedness to Dr C, S. Beals, not only for suggesting the problems treated and 
making possible the progress of the work, but for his constant guidance through 
the observational material on P Cygni stars and its interpretation. 


Dominion Observatory, 
Ottawa: 
1951 November 17. 
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Summary 

The general problem of the motion of a star system in the smoothed 
gravitational field which the system itself produces has been considered 
previously for the special case in which the stars are of equal mass. This 
restriction is now shown to be unnecessary. By a suitable definition of the 
velocity-distribution function the three equations which are involved in this 
problem can be written without any reference to the distribution of stars of 
different mass. The complexity of the problem is not affected by the 
increased generality. 

For a star system with complete spherical symmetry the problem is reduced 
to solving a single non-linear partial differential equation. A number of 
special solutions have been obtained, including solutions which correspond 
to isothermal and polytropic gravitating gas-spheres. For three different 
models the variation of the projected density with radius has been found, 
and this is compared with observations of the variation of luminosity with 
radius in globular clusters. ‘The best agreement is obtained in the case of 
a model which has never previously been considered, but which is essentially 
a modification of Plummer’s model in which terms involving the angular 
momentum have been introduced. 





1. General theory of self-gravitating star systems with stars of different mass.- 
The two-fold relationship between the velocity-distribution function of an 
isolated star system and the gravitational potential of the system itself has been 
pointed out in a previous paper on this subject (1). ‘The generalized Boltzmann 
equation expresses the fact that changes of the velocity-distribution are controlled 
by the gravitational forces. The Poisson equation relates the gravitational 
potential and the mass-density which produces it. Finally the mass-density 
is obtained from the velocity-distribution function by integration over the velocity 
components. In the previous paper it had been supposed that the system 
consisted of stars of equal mass. ‘This assumption is quite unnecessary. 
Obviously the results are unaffected if we assume, not that all the stars have the 
same mass, but that the average mass of a star is the same in all regions of the 
system. It now appears that, by a change in our definitions, we can formulate 
the problem in such a way as to include stars of very different masses, even when 
the distribution of masses is a function of position. Most fortunately, too, the 
change leaves the mathematical treatment quite unaffected. 

As before, we adopt any convenient non-accelerated Cartesian frame of space- 
coordinates, x, y, z, and corresponding velocity components u, v, w. Then the 
mass-velocity-distribution function, f(t, m, x, y, z,u, Vv, w), is such that the number 
of stars at time ¢ which have masses between m and m+ dm, which are in a volume 
dx dy dz, surrounding the point (x,y, z), and for which the velocity-components 
lie in the ranges, u to u+du, v to v+dv, w to w+ du, is f.dmdx dy dz du dv dw. 

Again we assume that the gravitational forces can be derived from a potential 
function (2(t, x,y,z). As a consequence, the Boltzmann equation now becomes 


of. of. af. af aa, ana ama 
ot + ax T° dy os. Oxdu* dydc + Jz dw (1) 
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From the definition of /, it is clear that p(t, x,y,z), the mass-density, is given 
by the equation 


p(t, x, y, 2) = JJ [fomte mM, X,Y, %,u, Vv, W) dm du dv du, (2) 
where the integrations are over the entire ranges of m,u,vandw. Finally, from 
Poisson’s equation we also have 

V2Q = — 4n7Go(t, x, y, 2). (3) 
The effect of non-uniform mass is removed if we define a function 


F(t, x,y, 2, U,V, W), 
such that 


F(t, x, y, 2, U, V, W) = | mf dm. 


If we multiply (1) by m, and integrate over all m, we obtain 


OF OF OF OF OQOF OQOF odRoF 
= tux tux +w— + ——4+——4F¢ 2 
ot Ox dy dz oxou dydv dz dw 


Similarly (2) can be written 


p=|]| F du dv dw. 


Equations (3), (4) and (5) contain the essence of the self-gravitating problem, 
but do not refer to individual masses in any way. ‘The mathematics of the 
previous paper is quite unaffected, except that the constant m is replaced by unity, 
and the function f which appears should be replaced by the function F’, defined 
above. 

It is not altogether surprising that the problem can be formulated without 
reference to the distribution of mass, for we are only concerned with the gravi- 
tational field of a smoothed-out distribution of matter. In such a field the 
acceleration of a star is independent of its mass. It is only when the effects of 
encounters are considered that stars of different masses experience different 
accelerations. It should be emphasized, of course, that in comparing theory 
and observation in this kind of investigation we usually compare the theoretical 
mass-density with either the number-density of observed stars, or the observed 
luminosity from a given element of volume. 

2. General theory for a spherical star system.—The treatment of star systems 
with spherical symmetry introduces one more variable into the fundamental 
equation than in the problem of a system stratified in plane parallel layers. In 
the latter case there were only three variables, namely the time, one space 
coordinate and one velocity coordinate; in the spherical problem, however, there 
are two velocity variables, namely the velocity components along and perpendicular 
to the radius vector. It is convenient to take spherical polar coordinates, r, #, ¢, 
with the origin at the centre of the system. The corresponding components 
of velocity are II, radially outwards, and © and ®, the components transverse 
to the radius vector. Since complete spherical symmetry is assumed, the 
components © and ® do not occur separately, but only in the combination 
©?+?, It is therefore simpler to put 7=0?+? (so that 7 is the square of 
the transverse speed), and to use 7 as the second velocity variable. ‘The square 
of the angular momentum per unit mass is r*7. 
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The integrated velocity-distribution function, F, now depends on the four 
variables, t,7, II and T only. The function is defined in such a way that at time 
t the sum of the masses of all stars in an element of volume dr at a distance r from 
the centre of the system which have radial velocity components between II and 
Ii+dII, and which have transverse velocity components whose squares lie 
between JT and 7'+dT, is F(t,r, II,7)drdildT. Thus F must satisfy an 
equation of continuity, namely 

OF OF diloF dToF 
ot Or dO dt oT 
The equations of motion are 


=0, 


dil T GM(t,r) 
S750 
dT 2uT 
ao 
In the first of these G is the gravitational constant and M(t,r) is the total mass 
at time ¢ inside the sphere of radiusr. ‘The equation of continuity now becomes 
oF oF (5 “) OF 2iiToF 


Fn += + ~ _ all _ ~.. =O. (6) 


The mass-density, p(t,7), is obtained from F by integration, so that we have 
etyr)=| | F(tr, I, T)all a. 


J 7=-0/ N=- 
Lastly, Poisson’s equation can be more simply written in the form 
0M/or = 4nr°*p(t,r). (8) 
By a suitable change of variables it is possible to replace (6), (7) and (8) by a 
single partial differential equation. In effect we take the Fourier transform 
of F with respect to the variable II, and the Laplace transform with respect to 7. 
It is therefore convenient to define a function M(t,r,s,p) by the equation 
M(t,r,s,p)=47| [ [| — x®F(t,x, IL, T)exp(ills—pT) dll dT dx. 
J2=0/ T=0/ M=-@ 
From the definition it will be seen that M(t, r, 0, 0) is the same as M(t, r), the mass 
at time ¢ mae the sphere vl aps r. It is also clear that 
ry o nat. 9,0:9)= { [" -F(t,r, 11, T)exp(ills—pT)dldT. (9) 
a: T=0/ Il=—a 
We now multiply (6) by exp(iils—p7), and integrate with respect to II and T. 
The following equations, all of which follow from (g), are required: 
1 0M 
4nr* drat’ 
0*M 
2nr dros =4nr*er*ds’ 
is OM 
exp(ills—p7T)dIldT = ~ Zeit Br’ 
_ is @M 
~ 4m drop’ 
i aM | ip &M 
4nr* ards * 4nr* Orosop* 


[ [ exp (ills—pT) alt , a 


{| [ 1 exp (ills—pT) dll aT = 


IE: 
|| 7 Sijexp (ills pT) dll aT 


{| TIl oF exp (ills— pT) dll dT = —, 
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With the aid of these equations, (6) is transformed, and after some simplification 
is written 

2M aM s2M_ GsaM 2p aM 

1 arat rds - 7 drop + oy sy M(t, r,O, 0) 7  drasdp =0. (10) 
This single equation is a necessary consequence of Poisson’s equation and 
Boltzmann’s equation. Conversely, any solution of (10) will lead, by application 
of the inverse Fourier and Laplace transformations, to functions F, p and Q 
which satisfy equations (6), (7), (8). ‘These will correspond to a possible physical 
model, if one condition is satisfied, namely that the function F(t, 7, I, 7) is non- 
negative. 

It should be emphasized that this treatment completely disregards the possible 
effects of close encounters. At a later stage we shall have to consider how impor- 
tant encounters are in the physical systems to which the theory will be applied, 
but as there is no immediate likelihood of obtaining an exact solution which 
allows for encounters, we shall examine equation (10) for solutions, even though 
in each the gravitational field is really that of a smoothed-out distribution of 
matter. 

3. Some solutions of the steady-state problem: the isothermal model.—The 
dynamical theory of spherical star systems has been studied quite extensively. 
Much of the early work is summarized by Martens in the introduction to his 
monograph on the subject (2). More recent contributions have been made by 
ten Bruggencate (3), Shiveshwarkar (4), Finlay-Freundlich (5, 6), and Kurth 
(7, 8,9). ‘Two fundamental results which reappear in the present paper are due 
to Plummer(1o) and to Eddington(11). Forthe most part these papers produced 
solutions which were too general, either because they did not require that the 
gravitational potential was due to the system itself, or because they were based 
on the virial theorem, which is less restrictive than the combination of the 
Boltzmann and Poisson equations. 

The aim of this present investigation is to find solutions of (10) which are as 
varied in type as possible. ‘There seems to be no clear way of finding the general 
solution, and we therefore have to obtain particular solutions by the method 
of separation of the variables. If this is necessary, then the use of transforms 
yields no advantage, and we can combine (6), (7) and (8) by more direct methods. 
The experience gained in studying the stratified star system suggests that we 
attempt to find a function F(t,r, II, 7) which satisfies the equations (6), (7) 
and (8), and which has the fairly general form 

F= saat { M M(t, nj .a(t,r).b(t,r). P{Ila(t, r), Tb(t, r)}. (11) 
In this, a and d are functions of ¢ and r, yet to be determined, and P is a function 
of p( =a) and q(= Tb), also to be determined. Without loss of generality we 


may assume 2 poe 
| P(p, q) dp dq =. 
: q=0/ p= — « 
It therefore follows that 


p(t,r)=| 
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Thus the expression (11) gives a form for F which satisfies both (7) and (8). 
The expression for F is substituted in (6), which becomes an equation in ¢,r, p 
and gq. By repeated differentiation we can separate the functions M,a,b and P, 
and so find each of them in turn, either exactly or else as the solution of a 
differential equation. ‘The method is quite straightforward, but the labour of 
examining all the possibilities arising is extremely tedious. ‘The conclusion of 
the matter is that there are only four solutions of type (11), and in none of them 
is the time explicitly involved. Consequently they all correspond to models 
in a steady state, and they are therefore of a special form in which F is solely a 
function of the two integrals of motion, namely the total energy per unit mass, 
3(11?+ 7) —Q(r), and the square of the angular momentum per unit mass, r?7. 
This result follows from a more general theorem due to Jeans (12), and has 
also been established for steady spherical systems by Shiveshwarkar (4). 
The simplest of the four solutions is one for which 


a® = b =1/2h? (a constant), 
P=7" exp { —(II? + T)/2h*}. 
‘The mass-density is related to the gravitational potential by the formula 
p =poexp {(Q —Q)/h*}, 


where py and {, are the values of the density and the potential at the origin. 
Now we also have 


dQ/dr = —GM/r, 


d (~ ” . 


from which it follows that 


dr 


= —4nGr°*p. 


dr dr 


The potential function 2 must therefore satisfy the differential equation 


1d dQ 
aa ( z) = — 47Gp,exp {(Q—Q,)/h?}. 
‘The actual form of F in this case is 
F =(87)-"*p,h- exp {(2Q —2Q, — Il? — T)/2h?}. 


Equation (12) is the equation satisfied by the gravitational potential of an isothermal 
gas sphere. It will be seen that the velocity-distribution is Maxwellian. ‘This 
is likely to be the form ultimately taken up by a spherical system after a verv 
large number of encounters have occurred. ‘There is essentially only one solution 
of (12) which will satisfy the boundary conditions:Q =Q,) and dQ/dr=o at r=o. 
(The second of these conditions merely expresses the fact that the gravitational 
attraction must vanish at the centre of the system. Models in which this attraction 
does not vanish at the centre have been considered by von Zeipel (13), but it 
seems unlikely that such models can correspond to a real situation.) In (12) 
it is convenient to write 


Q=Q,—h*u, 
Thus the equation becomes 
I a =) _ 47Gpoc* et, 


xtdx\” dx) 
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The constant ¢ is at our disposal, and can be chosen so that 
c? =h?/47Gpp. 


This makes the equation completely dimensionless, and we have 

1d du a 

de (= z) rod 
with the initial conditions that u and du/dx vanish at x=0. ‘This equation, 
originally studied by Emden (14), has more recently been considered for its 
astronomical implications by Chandrasekhar (15). It leads to a model of 
infinite radius and infinite total mass, which clearly will not correspond to a 
physical system. It may, however, represent the innermost part of a physical 
system, where the density is high. The general equation (10) applies equally 
well to a model where the form of velocity-distribution is fundamentally 
different in different regions—such as an isothermal core with an envelope 
satisfying quite a different law—but by restricting the form of solution to (11) 
we are prevented from finding such a combination. 

4. Eddington’s model.—The second solution of type (11) was originally 

obtained by Eddington (11), and is one in which a is constant, and 6 = a*(1 + ar?), 
where « is a constant. The function P is of the form 


P= Aq'exp {—3(p*+q)}. 
It is clearly necessary that the index yp is zero, for if » were negative the function 
F would be infinite for T’=0, so that at any point of the system there would be 
an infinite mass with a purely radial motion. Similarly, if ~ were positive, there 
would be no particles with purely radial motion. 
This solution leads to a relation between the density and the gravitational 
potential, namely 


—_ Po ; ee 2 
p= exp {(2-0,)/h*}. 


Consequently the function Q satisfies the equation 

I dQ 7G, 

ni (ft) =~ Peo 0-29 
This replaces (12) and has, of course, the same boundary conditions. The 
change of variables previously used will not remove the constant «, so that a one- 
parameter family of solutions can be obtained. From the form of Eddington’s 
distribution function it can be shown that there is no upper limit to the kinetic 
energy a star may have, and so there can be no “velocity of escape”: the model 
cannot have a finite boundary. ‘Thus it may be deduced that if the constant « 
is negative, the system is unlikely to correspond to a physical system. For if 
a is negative, (I +«r*) will vanish for a finite value of r, and the density is infinite 
there. 

If the constant « is zero, we recover the isothermal model previously considered. 

If « is positive, the solution corresponds to a model which has both an infinite 
radius and an infinite total mass. However, it differs from the isothermal model 
in two respects. ‘The velocity-dispersion in the previous case was the same in 
all directions and at all points. In this modification, due to Eddington, the 
velocity-dispersion in the radial direction is the same at all points, but the dis- 
persion of the transverse velocities decreases with increasing radius. Thus 
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there is star-streaming towards the centre at all eccentric points. Secondly, 
the two solutions differ in the way in which the density varies with radius for large 
values of r. In the isothermal model p(r) is of order r~* for large r, whereas 
in the modified form p is of order r-*(Inr)“". Thus the density tends to zero 
rather more rapidly, but not sufficiently rapidly to give a finite total mass. 

5. The polytropic model.—A third type of solution of the form (11) is one in 
which 6 =a?, and a is a function of r. The function P is now given by 


P=q‘(q+p? + constant)’, (13) 
or P=q'‘(constant — q—p*)*. (14) 
The arguments used in the previous section enable us to dismiss the cases in 
which Ao as being inapplicable to a physical system. In the case A=0 we 
shall certainly obtain a model of infinite radius unless we impose a finite limit on 


the energy. If we use the form (14) we can impose such a limit. In this case 
F is a function of the total energy only, and we have 


F=A(2Q-T?-Ty, 112+ 7 <2Q, 
F=0, Il? + T>2Q. 


Clearly the index » must be positive, otherwise F becomes indefinitely large as 
(II?+ 7) approaches 2Q. The function Q vanishes at the boundary, and is 
positive inside it. The density is obtained directly from F, namely 


p=A[|(2Q-IP-TydlldT, (2Q-I-T)>0 


Ar! =u! aii 
= res iy rey (2Q)"+82, 
In consequence of the Poisson equation, the function 2 must be a solution of 
the equation 


Ld (ot) — — 1422S! yn, (15) 


Pdr dr (w+ 3)! 
This can be simplified if we introduce the velocity-dispersion. By direct 
integration it will be found that 
TI? = 2Q/(2n +5), T=2i0. 

Hence there is no star-streaming, and the velocity-dispersion is isotropic at_all 
points, but it will decrease as r increases. If h? is the maximum value of II’, 
attained at r=0, we have 

202(0) = (2. + 5)h?. 
It is now convenient to change variables in (15), putting 

20 = (24. +5)h%, 

r= CX. 
If we choose the constant ¢ so as to satisfy the equation 
C=(m+h)! (uly At 2p+ 5) eh, 

we find that (15) reduces to the form 


I d dif 3/2 — 
ax(*%) +y"t+ =0. 


dx 
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This is the familiar Emden equation for the equilibrium of a polytropic gas sphere 
(14). ‘The polytropic index n is in our case the index (u+#). The boundary 
conditions are also those of Emden’s problem, namely 4 =1, dib/dx =0 at x=0. 
From Emden’s results we know that the star-system has a finite radius and finite 
mass if m<5, and that the radius is infinite but the total mass finite ifm=5. Inthe 
case in which n> 5, both radius and total mass are infinite. Effectively, therefore, 
there is a one-parameter family of solutions of type (14), with » as the parameter. 
The case n = 5, 4. = § is one for which (16) has an elementary solution, namely 
p=(r+ px?) ™, 
p=po(I + $x"). 
This special case has been considered by Plummer (10), and was found to agree 
fairly well with the density distribution observed in globular clusters. ‘This 
solution will be examined again in comparing these models with each other, and 
with observations. 
The alternative form (13) is one in which no limit is placed on the total energy 
a star may possess, ‘There may be stars with extremely high total energy, but it 
is reasonable to suppose that the number of stars of a given total energy tends to 


zero as the energy tends to infinity. This is arranged by taking as the distribution 
function 


F=A(II?+ T—2Q)". 


The index » must be positive to ensure that Fo as T,II> 0. It is also 
necessary that Q is everywhere negative, so that for all II and 7 the function F 
remains finite. It is then found that 


TI? =(—2Q)(2.-5),  T=2Tl?. 
In order that the mean-square velocity should be finite, we require 2u >5. ‘Thus 
the velocity-dispersion at any point is isotropic, and since the absolute value of 
Q increases with 7, the velocity-dispersion also increases with the radius. 
The differential equation satisfied by the potential function is analogous 


to (15). By a similar change of variables it can be reduced to the differential 
equation 

1df(/ dé ‘ 

= (= x) == fle, (17) 
This is a polytropic equation with negative index. We know that the model 
is of infinite radius, and by considering (17) it can be shown that for large values 
of r, p is of order r-*+’, where v=4/(2u—1). ‘Thus the density falls off more 
slowly than in the isothermal case, and for this reason we shall see that the solution 
under discussion does not correspond to a physical model. 

6. A modification of the polytropic model.—The remaining solution of the 
special type (11) we have been considering is related to the polytropic model 
in the way that Eddington’s is to the isothermal one. We have that 6 = a*(1 + ar*), 
where « is a constant, but a is a function of ry. As in the polytropic case, the 
function P takes one or other of the two forms (13) or (14). Again it is evident 
that the index A must be zero. ‘The function F is seen to be a power of a linear 
combination of the total energy and the square of the angular momentum. 
Solutions of this type have been considered by Kurth (7), though he has not 
exhausted all the possibilities of this model. 
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As before, to obtain a finite system it is necessary (though not sufficient) to 
take P to be of the form (14), so that F is given by the pair of equations 
F= A(2Q—-Il?- T-ar*T), IT? + T(t +ar*) <2Q, 
F=0, Il? + 7(1 + ar?) > 2Q. 


The index » must be positive. The potential function 2 vanishes at the boundary, 
and is positive inside it. The density is obtained by direct integration, and 
comes out to be 





= Anil? p! (2Q2)4+3/2 
sites (w+)! I+ar? * 


It therefore follows that Q satisfies the differential equation 
1d 2 dy ne 4: 4AnGp! (2Q)08 
r* dr dr (u+3)! I+ar? 
This can be simplified a little if we introduce po, the central density, and h?, the 
maximum value of II*, occurring atr=o. It is convenient to put r=cx, where 
ac?= +1, and also to put Q=ky. After substituting these new variables in 
(18), we obtain 





(18) 





a ae a = By " 9 >* 
x? dx dx h?*+3(2u +5) +32 | 1 + x? 


1d 24) a te yrs 


The constant k can be chosen so as to make the quantity in square brackets 
numerically equal to unity, so that we finally have 


. “= , as u+3/2 
( ; *) wo Eom (19) 


_—_— 
dx 


The boundary conditions for (18) were Q=Q,=A?(u+%) and dQ/dr=o at 
r=o0. After the transformations made, the boundary conditions of (19) are 


p=py=hu+4)/k and dj/dx=o at x=0. 


x? dx T+x2° 


‘Thus there are two parameters in the equation, the initial value 4) and the index pu. 
‘There is the further complication that the constant « may be positive or negative. 
If x is positive, the denominator on the right-hand side of (19) is (1+); if « 
is negative, the denominator is (I—x*). It should also be pointed out that the 
right-hand side of (19) is proportional to the density. Consequently, if a is 
negative, y must have a zero for some x, O0<x<1, otherwise the density will 
become infinite at x=1. ‘lhe only part of the solution which is of interest is 
that part which lies between x=0 and the first zero of ¥(x). The function % 
may have a mathematical existence beyond this point, but the density distribution 
to which this would correspond cannot influence the gravitational forces inside 
the critical radius. Further, since the density vanishes at this radius, there is 
no transfer of mass across the bounding sphere. For practical purposes, then, 
y decreases from yy to zero, and then remains zero for all larger values of x. 
The velocity-dispersion for this model is found to be 


Te =22 (2445), T= 211?/(1 +ar?). 
Again it is evident that the sign of « makes a real ditlerence. If « is positive, 
there is a radial star-streaming, as in Eddington’s model; if « is negative, the 
dispersion is smaller in the radial direction than in any other. 
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Equation (19) has been solved numerically for »=4 and for »=§ (corre- 
sponding to polytropic indices nm =2 and n= 4), and with the constant « taken to 
be negative. It was found that the starting value ys, had to be sufficiently large 
(i.e. « sufficiently small), otherwise % does not vanish before x=1, and so the 
density becomes infinite. ‘The permitted values of « are in fact so small that the 
variation of p with radius is scarcely distinguishable from that found in the 
corresponding Emden polytropic model. The small difference which is found 
tends to make the density of the modified polytrope fall less rapidly near the 
centre and more rapidly in the outer regions. ‘Thus the model is more like a 
polytrope of slightly lower index. 

The special case » = $, n=5 has not been solved for negative «. We know 
that if « is zero the density tends to zero only as r tends to infinity. It therefore 
seems very improbable that the introduction of any negative value of « will cause 
y to vanish before x=1. No formal proof of this conjecture is offered, but 
experience with lower indices showed that y will vanish before x=1 only if « 
is sufficiently small, and we know that for n= 5, if « is actually zero, y tends to zero 
only as x tends to infinity. 

For very similar reasons it has not seemed profitable to investigate any models 
with 1 >§ and « negative. There is little likelihood that they will avoid infinite 
density at x =I. 

Equation (19) has also been solved numerically for » = 4, $ and § (polytropic 
indices n=2, 4 and 5) and with various positive values of «. The differences 
from the Emden polytrope naturally increase with increasing «, but the nature 
of the solution does not change fundamentally. The models with »=4 and § 
have finite radii; those with »=$ donot. The equation has not been examined 
for values of 4 exceeding §. 

There remains the modified polytrope with negative index; this is the case 
in which the distribution function is 


F=A{Il?+ 7(1+ar?)—2Q}". 


In this form there is no upper limit;to the energy, so that » must be positive, and 
the model is of infinite radius. It is therefore essential that the constant « 
should be non-negative. 

The change of variables which gave rise to (19) now shows that the function 
Q is a constant multiple of 4, where y satisfies the differential equation 


Id/.d\_ ye 
xidx\" dx) 1+2x%° 


This equation has not been examined in detail, but it is found that for large 
values of x the function behaves like (Inx)"“-", It follows that the density 
will vary as x~*(In x)~*, where 


v=(u—§)/(u —4). 


Consequently the density tends to zero rather less rapidly than in Eddington’s 
modification of the isothermal model. 

This model, like the polytropic model with the negative index (of which 
it is a modification), also has the peculiarity that the value of I]? tends to infinity 
with increasing radius. 
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7. A comparison of the four theoretical models.—It is clearly important to 
compare in more detail the various models so far discussed. This can be done 
in a variety of ways for, in addition to the parameters ,« and « (when they occur), 
each model contains two further parameters. These could be taken to be any 
two of the following quantities :—the central density, the radius, if finite, the 
total mass, if finite, and the velocity-dispersion at the centre. Since, however, 
the radius and the total mass are infinite in some cases, they are best excluded, 
and the most convenient parameters are found to be the central density, pp, and 
an effective radius, a, at which the density falls to pyp/10. With these two para- 
meters we can plot p/p, as a function of r/a for any model, and all the curves 
obtained will pass through the two points (r/a = 0, p/py = 1) and (r/a =I, p/py =0°1). 

The first comparison made is between the polytropic model for which n=4 
and modifications of it when terms representing angular momentum are included. 
In this case the constant « can be positive or negative. ‘The variation of density 
with radius is not given here, since the differences between the three models 
are very small in this respect. The difference between them, such as it is, shows 
that the polytropic model (denoted by P) has a slightly more extensive region of 
high density near the centre than the model modified with a positive parameter « 
(the model denoted by M.P.+), and a less extensive region than M.P.—. ‘The 
difference can also be expressed by saying that M.P.+ shows a density distribution 
which is rather like an Emden polytrope which has an index a little greater than 4, 
whereas M.P.— gives a distribution like a simple polytrope of index a little less 
than 4. 

In Table I some further comparisons are given for these three models. ‘The 
dimensionless quantity «a? indicates the extent of the modification which dis- 
tinguishes (18) from (15). The degree of central condensation is indicated 
by the ratio of the mass within the sphere of radius a, denoted by M(a),'to the 
total mass. 

TABLE I 
A comparison of three models of polytropic index 4. P is the 
Emden polytrope, M.P.— and M.P.+- are modifications 

Model M.P.— P M.P.-4 
aa* —0°0245 ° +O°O415 
Radius/a 6°12 6°54 7°41 
Total mass /p,a* 1°882 1*883 1°868 
M(a)/'Total mass 0°574 0°577 
h?/Gp,a* 0°455 0°478 

(=II*/Gp,a*), > 

(Il*/Gp,a*),n¢ 0°265 0269 
(47/Gp,a*), = 0°272 0269 


The table also gives the values of A?, the velocity-dispersion at the centre, and 
in the last two lines the dispersion in the radial and transverse directions evaluated 
atr=a. In model P the velocity-distribution is spherical, but in M.P.— there 
is transverse star-streaming, and in M.P.+ the streaming is radial. 

The differences between these three models appear from the table to be 
extremely small, but the differences in velocity-distribution are considerable for 
large y—though at such distances the density will be extremely low. For the 
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Emden polytrope (P) the velocity-ellipsoid is in fact everywhere a sphere, but for 
M.P.— it becomes more oblate as we move farther from the centre of the system, 
and ultimately at the boundary the axes are in the ratio 1:1:0-29. On the 
other hand, the velocity-ellipsoid for M.P.+ becomes increasingly prolate, and 
the limiting ratio of the axes is 1:1: 1-81. 

It is to be remembered that all these three models are particular cases. ‘lhe 
properties of P are determined once the index n= 4 has been chosen, but in the 
two modified solutions the constant « is still at our disposal. ‘The negative 
value of « chosen for M.P.— is quite near the extreme limit when the denominator, 
1 +ar*, in (18) vanishes at the zero of Q(r). Thus the model M.P.— is almost 
as different from P as any choice of « can make it. ‘The case of M.P.+ is very 
different, since there does not appear to be any restriction on the size of « in this 
case. The results given in Table I for this model refer to a partieular solution 
which was determined after a value of /, particularly convenient for the numerical 
integration of (19) had been adopted. ‘The choice of a smaller value of 4 would 
have led to a larger «, and would presumably have increased the differences 
between P and M.P.+. 


TABLE II 
Comparisons between the isothermal model, Eddington’s modification 
of it, Plummer’s model, and a modification of it 


Model I MI Pl M.PIl 





Radius 
rla=o'o 
i 


Density, p(r)/po 
‘000 1'000 1*000 
823 0808 0863 
‘510 0°490 0°582 
°237 0°277 "337 
163 0°162 "184 
100 o*100 *100 
064 0°065 056 
"044 0°045 032 
O31 0°032 ‘O19 


e©o0000000% 


2 
“4 
‘6 
8 
‘o 
“4 
‘0 





M(a)/poa* 0°939 0°917 
di? Gpoa*)r =o o°*401 1°444 
(11?/Gpoa*), a 0°401 1°444 
(4T Gpoa"), <a O'401 0270 
(K.E../Gpoa*)y a o'6o1 0992 





The density distributions of these three models do not in fact seem to corre- 
spond very closely to observed density distributiens in globular clusters. In 
‘lable II is given a more detailed comparison of four models, three of which have 
previously been suggested at various times as representing the observed density 
distributions. ‘These are the isothermal model (denoted by J), the modification 
of it made by Eddington (./), the polytropic model of index n=5, suggested 
by Plummer and denoted here by Pi, and finally the modification of Plummer’s 
model by the inclusion of angular momentum terms (M.P/). ‘This last model, 
it will be recalled, is one in which the constant « must be positive. 

As before, the comparison is made between models with the same central 
density, pp, and with the same radius a at which the density falls to pg/10. All 
four models extend to infinity, and in all except P/ the total mass is also infinite. 
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The differences in density distribution are quite considerable; P/ has the most 
and M.] the least massive core, as is evident from the behaviour of p(r) and also 
from the values of M(a), the mass inside the sphere of radius a. Beyond r=a 
the density falls off most rapidly in Pl, and least rapidly in M.J. ‘This feature 
in which the density falls off more slowly in M.J than in J persists so far as the 
integration has been taken, but ultimately the position must be reversed, for it 
is known that in J the asymptotic behaviour of p(r) is like 7-*, whereas in M.J 
it is like r-?(Inr)-}. 

Table II also brings out the radial star-streaming in M.J and M.P1; it is 
quite marked even at r=a, which must be regarded as comparatively near the 
centre. In both these models the parameter « was chosen simply for convenience 
of numerical integration. It is clear that the values adopted are quite large 
enough to bring out the general effect of the extra-angular momentum terms. 

Another remarkable difference revealed in Table II is the mean kinetic 
energy per unit mass. At the centre this is $h? or $(II*),_); at radius a it is 
4(11?+7),_,. This latter quantity is given in the last row of Table II. We see 
that, both at the centre and at r =a, the mean kinetic energy per unit mass is much 
greater in the two models in which angular momentum terms have been introduced. 
This is connected with the fact that in these two models the individual star orbits 
tend to be longer and thinner than the orbits in the original isothermal and 
Plummer models. (The elongated orbits are responsible for the radial star- 
streaming, since a star tends to move, for a large part of its orbit, in directions 
inclined at only small angles to the radius vector.) ‘The high kinetic energy in 
M.I and M.P! also has some bearing on the stability of these systems. Chandra- 
sekhar has shown (16) that the effects of close encounters are much reduced if 
the mean kinetic energy is high, so that these models will be affected relatively 
slowly by close encounters in the central regions. It seems that the direct 
application of Chandrasekhar’s formula will actually exaggerate the influence 
of encounters, as the formula was derived for a star field of more or less uniform 
density. However, in M.J and M.P/ many of the stars which pass through the 
central core move in orbits which extend to regions of low density, and only a 
small proportion of the stars move permanently in a region of high density. 

8. The projected density.—For purposes of comparison with observation a 
knowledge of the theoretical space density is not of immediate value. What is 
observed is the projected density, v(r), the total mass of stars in a cylinder of 
unit cross-section with its axis at a distance r from the centre of the cluster. 
This quantity is calculated from the space-density by means of the formula 


ur)=2{- p(x) (20) 


ene (x? — 1°)!” 
For a model with a finite radius, p(x) is identically zero beyond a certain fixed 
value of x (x =c, say), and the upper limit of the integral is effectively c. Clearly 
in this case v is zero for all values of r greater than c. 

There is some difficulty in comparing theoretical models which have an 
infinite radius with the observations on globular clusters which are obviously 
finite. ‘There is, however, no reason why a globular cluster should not possess 
such a velocity-distribution for the greater part of its mass as would in theory 
belong to a cluster of infinite radius. ‘The density-distribution of a cluster 
depends on its own gravitational field, and this is provided by the material nearest 
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the centre. That part of a cluster which in theory may exist at large radial 
distances does not contribute to the gravitational field nearer the centre. It is 
rather the reservoir from which are replaced the stars of large kinetic energy which 
move out from the central regions. If the reservoir is absent, then there will 
be a slow escape of stars of high energy. 

An actual cluster is subject to the gravitational field of the galaxy to which it 
belongs, in addition to its own gravity. The external field will dominate the 
movements of all stars beyond a certain distance from the centre of the cluster, 
and therefore only the stars inside this boundary can properly be said to belong 
to the cluster. ‘This fact is not allowed for in the elementary treatment of the 
dynamics of clusters, and consequently we cannot expect to obtain theoretical 
models which agree with the observed density distribution in the extreme outer 
regions. ‘This is a great drawback to the method developed by Kurth (7). He 
assumes a law of density-distribution of the form p« r~", for sufficiently large 
values of r. Unfortunately the density in the outer regions of a cluster—where 
such a theoretical distribution is made to fit the observations—is only in part 
controlled by the field of the cluster itself. 

In the present investigation the question of the finite or infinite radius is 
important in calculating the projected density. In practice the integral (20) is 
evaluated numerically up to some convenient upper limit of x, beyond which 
p(x) has not been found. A correction can then be applied which will effectively 
extend the upper limit of the integral to infinity. In all cases the limit to which 
the integral was evaluated directly was about 3a, varying slightly as between 
different models for the convenience of integration. At this distance the ratio 
p/Po has dropped to about 0-003. The behaviour of »(r) will be shown both for 
corrected and uncorrected models. It turns out that cut-off at r=3a which we 
arbitrarily imposed is much too near the centre, but the effect produced on 
v(r) by making a cut at this distance permits us to predict the effect of terminating 
a model at a much greater radius. 

‘The comparison of the projected densities of different models cannot be drawn 
until we fix a scale for the radius. ‘This is most easily done by adjusting the 
constants of each system so that all models have the same central projected 
density vy and so that at the same radius a’ the projected density falls to 9/10. 
The graph of v(r)/v_ as a function of r/a’ is given for Pl, M.Pl and J in Figs. 1, 
2and 3. Itis not very convenient to represent all the curves on the same diagram, 
but it is not difficult to compare the different models, since in each figure the same 
twenty-four points (obtained from observations) are plotted. Since the 
projected density varies by a factor of 100, it is difficult to represent the ordinates 
on a uniform scale. It seems advisable therefore to multiply the scale of the 
ordinates by 10 for all radii r exceeding a’. In each diagram there are two curves ; 
the continuous one is the graph of v(r)/v, for the model of infinite radius. ‘The 
broken curve is the graph of the same function when the contribution both to 
v(r) and to vy of the material which lies outside a certain fixed sphere is completely 
neglected. Naturally, the effect of this cut-off depends very much on the space- 
density, and especially on its variation at great distances. In PI, for instance, 
p behaves as r~® for large values of r, and consequently the outer regions really 
contribute very little to the integral for v(r). It follows that the two curves ror 
this model are extremely close together until the radius approaches the limit of 
the finite model. On the small scale used for r<a’ the two curves are indistin- 
guishable. In M.P/ and I the density behaves for large r as r-*(Inr)~** and 
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r~*, and in these two models the effect of the cut-off is more obvious. In both 
cases the finite model will appear to have a more extensive central core, but the 
projected density will fall off more rapidly in the outer regions. 

9. The comparison of these models with observation.—The twenty-four points 
which are included in Figs. 1, 2 and 3 correspond to observations made by 
Lohmann (17) on three globular clusters, NGC 5904, NGC 6341 and NGC 7078. 
(These clusters are also known as M5, Mg2 and M15, but the symbol M is in 
such demand that this nomenclature will be avoided.) Lohmann has measured 
the luminosity emitted per unit area at a large number of points on the disks 
which these clusters present. Instead of giving his results in terms of apparent 
magnitudes, he estimates the number of stars of apparent magnitude 14™ which 
would be required in each area to produce the measured luminosity. If the 
fractional distribution of stellar masses is the same throughout a cluster, then 
the luminosity of any element of volume would be proportional to the total 
mass which the volume contains. In these circumstances Lohmann’s figures 
for the luminosity (represented as a number of stars) would correspond to 
the projected density v(r) which has been calculated for theoretical models. 
Lohmann has also made actual counts of the bright stars in these clusters, and so 
has found the distribution of these stars with radius. He finds that this does 
not correspond exactly with the distribution of luminosity, so it cannot be supposed 
that the proportion of stars of given mass is the same at all points. On the other 
hand, the differences in distribution which Lohmann finds are fairly small, and 
it seems likely that the distribution of mass differs very little indeed from the 
distribution of luminosity. 

Taking Lohmann’s figures as being proportional to v(r), we find by inter- 
polation the angular radius in each cluster at which r(7)/vp falls to 4; this will 
be the standard radius a’, now measured as an angle. ‘The observations are then 
represented as points on Figs. 1, 2 and 3, in which v(r)/vp9 is plotted against r/a’. 
Three different symbols have been used, so as to distinguish between the three 
clusters. The same points occur in all three diagrams, and serve as a guide in 
comparing the theoretical curves. ‘The remarkable feature of the observations, 
and a feature which was not pointed out by Lohmann, is that within the limits 
of observational accuracy the points obtained from all three clusters appear to 
lie on the same curve. ‘This uniformity is to some extent confirmed by Shapley’s 
figures for the total angular diameters (18). From these and Lohmann’s results 
we find that the boundaries of the clusters are at radii corresponding to 6-7a’, 
5:5@’ and 8-4a’. Apart from the observational uncertainty of measuring these 
radii, it is scarcely to be expected that they would be in much closer agreement. 
It has already been pointed out that a model which is theoretically infinite in 
extent will be modified in practice by the gravitational field of the galaxy. Thus 
the position of the finite boundary of a cluster will be determined not merely 
by the mass distribution in the cluster itself but also by the past motion of the 
cluster as a part of the galaxy. Shapley’s values of the angular diameters therefore 
seem to be quite consistent with the suggestion that the three clusters are dynami- 
cally similar. 

The comparison between Lohmann’s observations and the theoretical curves 
we have obtained brings out a number of interesting points. In the first place 
it is obvious that the finite boundary, whose arbitrary imposition leads to the 
broken curves, is in all cases too near the centre. The point is stressed by 
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Shapley’s measurements, which put the finite boundary at about three times the 
radius imposed on the theoretical models. The finite models serve only to 
indicate the general effect of cutting off an infinite model at some finite distance. 

Of the three infinite models it is evident from Fig. 1 that Plummer’s model, 
Pl, does not fit these observations. The density does not fall off sufficiently 
rapidly near the centre, but it becomes too small in the outer regions. 

The isothermal model, J, shows in Fig. 3 much better agreement near the 
centre, but the density still does not decrease sufficiently rapidly in this area. 
Beyond r =a’ the projected density falls off much too slowly, and the curve in no 
way fits the observations. ‘This part of the curve can, of course, be modified 
considerably if we assume a finite boundary. (Chandrasekhar (1§) gives values 
for v(r) for an isothermal model with a finite boundary, but his bounding radius 
is little more than half the one chosen here, and this in turn is probably only about 
a third of the value required to fit Shapley’s estimates of the bounding radii.) 
The effect of cutting off the extreme outer layers is to reduce the projected density 
for r>a’, though the difference between theory and observation near the centre 
will be slightly increased. It seems rather unlikely that a boundary of radius 7a’ 
will so modify the infinite model as to bring about a real agreement between 
theory and observation. 

The model M.P/ which is derived by introducing angular momentum terms 
into Plummer’s model, the polytrope »=5, seems to fit the observations very 
closely in Fig. 2, except for the usual discrepancy in the central area. For r>a’ 
there is a general tendency for the curve to be above the points. This would 
be corrected to some extent by the exclusion of the mass beyond r =7a’, and it 
seems that the curve would then fit the observations remarkably well. It will 
be recalled that this model is derived from equation (19) in which there are two 
parameters. ‘The index »= 4, (n=5) clearly fixes the main properties of the 
model, but the initial value of ys at x=0 is also at our disposal. (This quantity 
yy is related to «, which is the constant controlling the extent to which the angular 
momentum terms are involved.) . For convenience of integration y%, was first 
taken to be 1; then to investigate the effect of changing this constant, a second 
solution was obtained with %,=0-9. ‘This agreed even more closely with the 
observations than did the first solution, and it is this model which has been used 
in Table II and in Fig. 2. Evidently we might expect to improve on this solution 
by changing the starting constant y, but it is doubtful whether the observations 
are sufficiently accurate to justify the extra labour. 

This comparison with observation, revealing a good agreement with M.Pl 
but not so good in the case of the others, might be checked in another way. If we 
take the case of the cluster NGC 5904, for example, we can find the total mass 
inside a cylinder of radius 2-0a’ from the integral 

r2-0a’ 


2r | rv(r) dr. 


This quantity can be estimated for all three models; it corresponds, as it were, 
to the moment of the area under the curves in Figs. 1, 2 and 3 about the line r =o. 
The curves are sufficiently alike for all to give about the same value. This is 
the mass inside the cylinder examined by Lohmann, and he associates with it 
an integrated apparent magnitude of 6"-8. In the same way we can obtain 
three estimates of the mass inside a cylinder of radius 6-7a’, the estimated radius 
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of NGC 5904. Since the three curves of »(r) are very different beyond r=2:0a’, 
these values are very different. Since we know the integrated magnitude of 
the part we can now calculate three values for the integrated magnitude of the 
whole, and this can be compared with Shapley’s observed value (18). Here the 
check breaks down, for Shapley’s value is 7™, suggesting that the whole is fainter 
than the part! There must be some discrepancy in the measurements of inte- 
grated magnitude, and this obviously prevents us from making this check. It 
should be added that Lohmann’s estimates of the integrated magnitudes of the 
regions he has observed are in reasonable agreement with the values given by 
Vyssotsky and Williams (19). ‘These are slightly brighter than Lohmann’s, 
but they refer to a slightly larger region. 

It is to be hoped that this discrepancy can soon be resolved, for these obser- 
vations of integrated magnitudes provide a valuable test. ‘The three models 
give entirely different results, as can be seen from ‘Table III. In this is given, 
for each cluster, the fraction of the total mass which lies in the cylinder examined 
by Lohmann. ‘This fraction can, of course, be evaluated according to each of 
the three models. 

Taste III 
The fraction of the total mass which is observed by Lohmann, 
as calculated for each model 

Model Pl M.PIl I 

Cluster 
NGC 5904 o'9gol 0582 0°323 
NGC 6341 0'937 0°694 0°459 
NGC 7078 0950 0°634 0°346 


From these ratios it is clear that if the clusters have a density distribution like 
the M.PI/ model, then the integrated magnitude of the part observed by Lohmann 
should be in each case about half a magnitude fainter than that of the whole 
cluster. For the isothermal model the difference is about one magnitude. 

10. The velocity-dispersion in the clusters.—In several dynamical theories of 
globular clusters it is possible to obtain a formula connecting either the total mass 
or the central density with the velocity-dispersion. Unfortunately we have no 
direct knowledge of the central densities, and though there are estimates of the 
total masses and of the velocity-dispersion, they are still very uncertain. To 
obtain the mass of a cluster by a method which does not involve dynamics it 
seems clear that we must use absolute luminosities. ‘The masses of NGC 5272 
(M3) and NGC 6205 (M13) have been determined in this way by Lohmann (20). 
The absolute photovisual magnitude of each is obtained by correcting Shapley’s 
absolute photographic magnitude for colour. For both clusters the quantity 
comes out to be —8™-8. Some of the luminosity is directly attributable to the 
bright stars visible in the cluster. ‘Thus the number of stars in each of the 
brighter magnitude classes is known. Lohmann assumes that the rest of the 
luminosity is due to a distribution of faint stars comparable with the distribution 
of absolute magnitudes in the neighbourhood of the Sun, as found by Kuiper 
(21,22). In this way Lohmann deduces the distribution of absolute magnitudes 
in each cluster, and hence the distribution of stellar masses. Actually Kuiper 
gives two different distribution functions, and Lohmann uses each in turn. In 
this way he derives for NGC 5272 a total mass of 6-4 x 10*© or 7°7 x 1040 ; 
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and for NGC 6205, 9:7 x 10© or 11-8 x 1040. The difference between the two 
clusters arises from the fact that NGC 5272 has more intrinsically bright stars, 
and NGC6205 must therefore have many more faint stars to reach the same 
total absolute luminosity. 

On this evidence it is reasonable to assume that a cluster of absolute magnitude 
— 8™-8 has a total mass of 9g x 108. There may easily be an error of 30 per cent in 
either direction, but only a cluster with exceptionally few bright stars will have 
the same absolute magnitude and a mass exceeding 12x 10'©. This rough 
approximation can be applied to NGC5904, NGC6341 and NGC 7078, the 
three clusters considered in Lohmann’s earlier paper. For each the absolute 
photographic magnitude has been given by Shapley (18), and the colour index 
by Becker (23). Thus we can obtain the absolute photovisual magnitudes, It is 
assumed further that the distribution of stars of different brightness is the same 
in all clusters. It therefore follows that the mass of a cluster of absolute photo- 
visual magnitude M is N solar masses, where N is given by the formula 


N=Q x 104 x 10-0& +88), 


The masses of the three clusters determined in this way are given in the first line 
of Table lV. For NGC 6341 Hachenberg (24) has obtained a mass of 25 x 1040, 
over three times the value given here. It is difficult to see how this mass could 
have such a faint absolute magnitude. 


TABLE IV 


Above the line are constants adopted from observational data on the three clusters. 
quantities below the line are calculated from them on the assumption 
that the clusters are of M.PI type 
Cluster NGC 5904 NGC 6341 NGC 7078 
‘Total mass 7°2 X 10 7°2x 10'O 9°5 x 10°O 
Distance 10°! kpc 10°3 kpc 11°5 kpc 
Radius 36°7 pe 21°0 pe 30°3 pe 
Radius at which the projected 
density is v»/10 5°52 pe 3°84 pe 3°59 pe 
Acceleration due to gravity at 
boundary ° 10-19 cm/s? 2:28 107® cm/s?___-1°45 X 10~* cm/s? 








Central density 330 O/pc? 1070 O/pc? 1460 O/pc?® 
h 4°4 km/s 5°5 km/s 6-0 km/s 
h(o) 4°2 km/s 5°3 km/s 5°8'km/s 
h(a’) 2‘2 km/s 2°7 km/s 3°0 km/s 
Tx, time of relaxation 1°08 x 10° yr. 7°6 x 108 yr. 6-2 x 108 yr. 


The distances and the angular diameters have been given by Shapley, and so 
the linear radius of each cluster is readily found. ‘The angular radius at which 
the projected density falls to a tenth of its central value is interpolated from 
Lohmann’s results, and converted to linear measure. This determines the 
scale of the model considered. On the assumption that the clusters have the 
properties of the particular M.P/ model we have examined, the central density 
and the velocity-dispersion can now be calculated. 

To obtain the central density we make use of the total mass. ‘The mass inside 
a sphere of radius R is given by the expression R®?G~!( —dQ/dr),_p. On changing 

13 
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the variables in the way which transforms (18) to (19), we find that this expression 
can also be written as 


4m yy °( — x? dib/dx) o—Rie- 


We also have for this particular solution yj, =0-9, a’/c=2-696. Hence py, can be 
found for each of the clusters. It is very noticeable that for NGC 5904 the central 
density is only about a quarter of the value obtained for the other two clusters. 
This is seen to be a consequence of the large value of a’. ‘The density of this 
cluster falls off much more slowly than in either of the others, and consequently 
the same total mass can be attained with a much lower central density. 

The root-mean-square velocity at the centre, denoted by A, is now found from 
the formula 

h? = 0-439Gp,a’*. 


Unfortunately the quantity h cannot be observed directly. What we hope will 
ultimately be determined is the root-mean-square velocity in the line of sight for a 
cylindrical column of material passing through the centre. If this quantity be 
denoted by A(o), then we have 


ea -o 


ho) =h? | ” p(x) (x) dx / Q(0) | p(x) dx. 


Similarly the root-mean-square velocity in the line of sight may also be found for 
a cylindrical column with its axis at a distance 7 from the centre. Denoting this 
quantity by A(r), we have the relation 





, he 2 x dx ey x dx 
na(e) =H | o(x)a(a) (1-5) sre / 2(0)| ” o8) Ga Sar 


+ ax" 


In this case the velocity in the line of sight in the region of greatest density is not 
radial to the cluster but is tangential. We have already stressed that the disper- 
sion falls off much more rapidly for the tangential components than for the radial 
one. It follows that h(r) decreases quite rapidly with increasing r. The 
quantities 4(o) and h(a’) have been evaluated for the three clusters, and appear 
in Table IV. 

The only observational material on velocity-dispersion in clusters is given by 
Joy (25). He gives the radial velocities of a number of variable stars in several 
globular clusters. ‘There is no indication of the part of the clusters in which the 
stars are found, but it seems that the dispersion in the velocity measurements is 
likely to correspond more nearly with the theoretical quantity 4*(a’), the dispersion 
in a cylinder at a distance a’ from the centre, than with h?(0). For NGC5904 
Kurth estimates from Joy’s results that the root-mean-square velocity is 2-00 km/s. 
This agrees quite satisfactorily with the theoretical figure of 2-2km/s. For 
NGC 6341 no results are available. For NGC 7078 the value deduced by Kurth 
is 15°5km/s. This is far away from the theoretical value of 3-0. However, it 
should be mentioned that whereas in the case of NGC 5904 the dispersion is 
calculated from 22 variables, 14 of which could be given full weight, the dispersion 
in NGC 7078 is obtained from two variables only, neither of which could be 
given full weight. ‘The disagreement between theory and observation in this case 
need not therefore be taken very seriously. Indeed, Joy writes: ‘‘ The velocities 
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for individual stars are clearly of insufficient precision to justify any conclusions 
as to peculiar motions or rotation effects.”’ It therefore seems that the observa- 
tions on velocity-dispersion should be introduced only at the end of the calculations ; 
it would be dangerous to work backwards from these observed dispersions and 
so obtain the total masses. Such a method has been adopted by Kurth, who is 
thus led to the conclusion that the mass of NGC 7078 is a hundred times that of 
NGC 5904. If this is so, it is very difficult to account for the absolute magnitudes, 
which differ by no more than o™-3. 

Two further points need to be made. Whatever the dynamical model of 
the clusters may be, we can calculate directly from the mass and the radius the 
acceleration due to gravity at the boundary. (This is also given in ‘l'able IV.) 
This is only a fraction of the acceleration due to the galaxy in the direction normal 
toits plane. From the density variation near the plane a value of 6-75 x 10~* cm/s? 
has been found (1), and the work of Oort and van Woerkom (26) leads to a value 
8-4 x 107-* cm/s? up to 10 000 parsecs from the galactic plane. Although the three 
clusters considered are in high galactic latitudes, it seems quite certain that at 
the boundaries of the clusters the attraction of the galaxy is the dominant force. 

Lastly it is both interesting and important to calculate the time of relaxation 
of these models, the time which will elapse before close encounters cause a major 
redistribution of kinetic energy. Chandrasekhar’s formula for this time Ty 
(16) is based on the assumption that the distribution of velocities is Maxwellian, 
whereas this is not so in the model we have used. However, it seems likely that 
the formula will still give the right order of magnitude. In using the formula 
we have taken the average mass of a star to be 0-5, as was obtained by Lohmann. 
It should also be noted that Chandrasekhar’s quantity v? is the mean-square 
speed, so that its value in our notation is 3h? at the centre of acluster. The actual 
values of T',, given in the last line of ‘Table IV, are not quite large enough for us 
to infer at once that close encounters play a negligible part in the dynamics of 
clusters, but it must also be remembered that only a small part of the material 
within a sphere of radius a’/10, say, is permanently within that sphere. Chandra- 
sekhar’s formula is based on the cumulative effects of close encounters, but in 
the model considered these will not be frequent. ‘The small energy changes will 
diffuse fairly rapidly through the system, but will not be sufficiently frequent in 
the history of any one star to alter its energy really drastically. ‘This is the real 
justification for the dynamical methods here used, ignoring, as they do, the 
effects of close encounters. 

We therefore conclude that the modification of Plummer’s model by the 
introduction of angular momentum terms leads to a model which is dynamically 
self-consistent, in the sense that the gravitational field which it produces is such 
as will maintain the velocity-distribution exactly. ‘This model appears to fit 
the variation of brightness of the observed disk, as measured by Lohmann for three 
globular clusters. It also provides a relation between the total mass and the 
velocity-dispersion which is in general agreement with the estimates which we 
have for these quantities. 


The University, 
Manchester, 13: 
1951 November 30. 
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THE OXFORD SOLAR SPECTROSCOPE 
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Summary 


The prismatic spectroscope, used with the Oxford solar telescope, has a 
clear aperture of 15:2 cm and a focal length of 8-96 m. The linear dispersion 
varies from 2°8 mm A~! at A 4000 to 0-5 mm A~! at A6500. The horizontal 
apparatus function, determined from the Kr line A 4319°6, has a half-width 
of 0-048 A with a marked asymmetry to the red. The practical resolving 
power at this wave-length is 93 900, less than one-third of that predicted 
by diffraction theory. Because of the asymmetry in the apparatus function 
the centres of intensity of solar absorption lines are displaced to the red by 
amounts which vary as their equivalent widths; precautions must therefore 
be adopted in the measurement of Doppler displacements. ‘The vertical 
apparatus function (up and down the slit) has a half-width of 0-167 mm. ‘This 
tends to wash out the spectra of fine details on the solar surface, but theory 
and practical experience show that spectra of solar granules can be obtained. 
For unanalysed white light, incident on a slit with a height of 19 mm, the 
general scattered light in the focal plane is 1-6 per cent. 





The solar equipment of the University Observatory, Oxford, consists of a 
solar telescope of 32cm aperture and 19-°8m focal length (13), and a large 


prismatic spectroscope. ‘This latter instrument is described, its performance 
determined and its limitations noted in the three succeeding sections of this 
paper. 

1. Description 

This auto-collimating spectroscope has a clear aperture of 15-2. cm (6in.) and 
a focal length of 8-96 m (29°4ft.). Dispersion is provided by three prisms, two 
of 60° and one, silvered on its back face, of 30°. These prisms have a height 
of 15:2cm (6in.) and a base of 24-8cm (9°8in.). The dispersion varies from 
approximately 2-8 mm A! at A 4000 to 0-5 mm A=! at 46500. When allowance 
is made for absorption in the prisms the theoretical resolving power at A 4320 
is 303 000. 

The spectroscope is mounted horizontally in a brick tunnel in the basement 
of the Observatory. ‘Thé beam from the solar telescope is turned horizontally 
to the west by a smail flat C (Fig..1) and the solar image is focused in the plane of 
the slit. ‘The slit and plate-holder unit is mounted on the T-shaped pier, B. 
The beam, emerging from the slit, passes into the tunnel and falls on the objective 
and prisms mounted on pier A. Reflected off the silvered face of the 30°-prism 
the beam retraces its course and the focused spectrum is photographed on the 
same level as, and very nearly in the plane of the slit. 

1.1. Mechanical parts.—The slit and the plate holder are mounted in a 
massive, hollow, cast-iron box some 45cm wide by 28cm high by 33cm deep 
(17-5 x 11 x 131n.°). This heavy box rests on a casting D (Plate 2, 1) which is 
bolted through spring washers and a cork pad to the cross-bar of the T-shaped 
pier; for collimation in the horizontal plane the box may be moved north or 
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south on its base D, and finally locked in position, by four push-screws. The 
box, open on the side towards the objective and prisms, is connected by a light- 
tight flexible bellows to a corresponding opening in the east wall of the tunnel. 

The slit, E, has jaws 5 cm high, one of which may be moved by a half milli- 
metre pitch screw against a spring to give an opening of any required width. 
The height of slit exposed is determined by one of a number of slots in a diaphragm 
which slides in ways almost in contact with the front face of the slit-jaws. The 
slit is carried by a tube some 6cm in diameter and 33cm long, in which it is 
rotatable through a small angle. ‘This tube passes through the front face of the 
metal box to which it is rigidly attached by a heavy metal flange. The other 
end of this tube, well inside the metal box, is closed by a rectangular aperture which 
passes the beam, here some 5 mm wide, on its way from the slit to the objective. 
Normally this aperture, and consequently the beam from the slit, are completely 
unobstructed, but by operating a lever on the front face of the metal box a mask 
may be swung into the beam so as to shadow either one or other half of the distant 
objective. ‘This provides a simple means of applying the Hartmann test for 
focal setting. 
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Fic. 1.—-Spectroscope tunnel. (Scale 1 cm=2-3 m.) 


A gun-metal plate, F, some 33 x 33 cm? is pivoted at G and H to the cast-iron 
box. A carrier for the plate holder, 1, may be moved in ways up and down this 
gun-metal plate by means of a handle, J. Rectangular openings, on the same 
horizontal level as the slit, in both the cast-iron box and the gun-metal plate allow 
the spectrum to fall on the photographic plate; the passage between the iron box 
and the gun-metal plate is rendered light-tight by a short length of bellows. 
The gun-metal plate may be rotated on its pivots through +5° so as to set the 
photographic plate tangential to the focal surface of the objective. ‘Three gun- 
metal plate holders are provided; each takes a plate 25-4 x 10cm? (10 x 4in.°), 
inserted by removing the back, and has carrying handles and lugs for clipping it 
into the carrier, as well as the usual dark slide. 

The objective and the prisms are mounted on a heavy ribbed casting, K 
(Plate 2,2)*, of surface area I-10 x 1-17 m* (43 x 46in.*). Four massive screws 
at the corner of this casting provide collimation in the vertical plane; these screws 
are seated in four metal bosses, spigoted into a 5-cm thick cork pad which covers 
the whole top of the pier, A(Fig. 1). To the casting, K (Plate 2,2), motionless 
on the pier by virtue of its weight, is pivoted a sub-base. This latter, adjustable 
for angle in the horizontal plane, carries the slides for the prism cells, L, M, N 
and the objective cell, O. The three triangular slides, to which the prism cells 


* Reproduced from Engineering, p. 335, 1939 March 24, where further details on the mechanical 
design will be found. 





Montuty Notices or R.A.S. Vor 112, Prate 2. 


1 es 
Cs 


1. Slit and plate-holder unit. 

2. Objective and prism-carrying unit, 

3- Fouc ault test of 30°-prism and OG, 1 Daly edlen a wate. 
4. Foucault test of 3 prisms and O.G. | * 
5. Typical scattered light plate. 

6. Meg b region 1°67 with full aperture. 

7. Mg b region © 1°67 with diaphragm in position. 


H.H. Plaskett, The Oxford solar spectroscope. 











The Oxford solar spectroscope 179 


may be interchangeably screwed, are carried on a Browning—Proctor minimum 
deviation link work (8) operated by a screw and carriage (not in position when the 
photograph was taken). The objective cell, O, slightly adjustable for angle on 
its slide, may be moved for focusing through 15cm by a screw and nut; its 
setting is read off a scale. 

A metal-lined three-ply box, itself lined with cork to a depth of 5 cm, sits on 
the cork pad on the top of the pier, and completely encloses, but does not touch, 
the objective- and prism-carrying unit. ‘This box, of which the back and a 
part of the top are removable to give access to the optical parts, has a base of 
1-47 x 1-40 m? (58 x 55in.®) and a height of 0-79m (31in.). A circular hole in 
the east wall of this box allows the beam from the slit to enter the objective and 
to return to the photographic plate; a bellows connects the front of the lens 
cell with the inside of this box, thus making the front component of the lens a 
window into the enclosure. ‘The spindles of the screws for focusing and changing 
the position of minimum deviation are carried through seals in the east wall of 
the box. Sprocket wheels on these spindles engage corresponding sprocket 
wheels on two shafts which run in self-aligning ball bearings along the north 
and south walls of the tunnel. These shafts, again through sprocket wheels, 
can be rotated by hand wheels on the cross-bar of the T-shaped pier (B, Fig. 1). 
Engagement of the sprocket wheels is provided by loose, flexible chains to 
minimize any transmission of vibration. The focal setting of the objective, read 
by a small telescope just above pier B, and the wave-length of minimum deviation 
can thus be altered at will from this pier without entering the tunnel. 

Other features, for example an optical bench on the stem of the T-shaped 
pier, are included in the mechanical design. ‘The details of this design, as well 
as the actual construction, are due to C. F. Casella & Co.; the proved stability 
of the instrument and its convenience and flexibility in use are a tribute to their 
skill and the excellence of their workmanship. 

1.2. Optical parts.—The prisms and the objective (with their cells), and the 
spectroscope slit, are due to Adam Hilger & Co. The prisms, cut from a large 
block of optical glass (just as it came from the fire-clay pot) and annealed after 
cutting, are of dense flint with refractive index, dispersion and transmission almost 
identical with the Jena Glass (Schott & Gen.) of type F3. ‘To minimize the 
effects of internal inhomogeneity the faces of the prisms are locally figured. 
The achromatic objective of 15-2 cm aperture is a special design of J. W. Perry, 
Head of the Applied Optics Department at Hilger’s. 1t consists of two single 
components, separated by 22cm, to reduce the chromatic difference of spherical 
aberration. ‘The minimum focal length of 895-8cm occurs at 45100 with an 
increase for single transmission of 1-5 cm at A 8000 and of 2°8cm atA4o000. ‘The 
radius of curvature of the focal surface is estimated from the design to lie between 
two and three metres. ‘Three of the four surfaces of the lens are concave to the 
slit so that the real, white-light images of the slit can be masked off to reduce 
the general scattered light. 

" The spectroscope was first assembled as a unit under good conditions of 
temperature control (see next sub-section) in 1938 September. It was immedi- 
ately evident that the definition was poor. ‘The top and the bottom halves of 
the optical train gave spectra which did not superpose and which focused in 
planes differing by 0-7cm. Even when one-half of the train was masked off 
the definition was not good, and Hilger’s willingly agreed to take the prisms 
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back for further examination and possible re-figuring. On examination with 
a large aperture interferometer, not available when the prisms were originally 
figured, Hilger’s found that the 30°-prism was satisfactory, but that further 
work could usefully be done on the two 60°-prisms. On the completion of this 
work the prisms and the objective were bloomed and returned to Oxford in 
1946. On re-examination in the spectroscope a marked improvement in de- 
finition was found; in particular the difference between the top and bottom halves 
of the optical train had disappeared. Dr Adam used the spectroscope in this 
form to obtain interferometric spectra (1) in the autumn of 1946. 

It seemed probable, however, that some further improvement in definition 
could be secured by masking off part of the optical train. Knife-edge tests in 
the focal plane, using the Na D, line from a Philips Lamp as a source (Plate 
2, 3——30°-prism alone, Plate 2, 4—all three prisms in train) fully confirmed 
Hilger’s conclusion that the 30°-prism was satisfactory, and that the remaining 
trouble arose from abrupt changes of refractive index in the two 60°-prisms. 
Under these circumstances almost any masking seemed bound to effect some 
improvement, but the best results were obtained early in 1948 with a diametral 
rectangular aperture of the full width of the prisms (15-2cm) and of vertical 
height 3-tcm. ‘The lower edge of this mask also casts a shadow on the photo- 
graphic plate for the virtual image of the slit produced by the one convex lens 
surface. ‘The effect of this mask on the definition of the solar spectrum can be 
judged from Plate 2,6,7. The upper of these two spectra was taken with 
the full aperture, and the lower with the diaphragm in position. It can be seen, 
though unfortunately not in the reproduction, that the Mg—Fe doublet at 45167 
(separation 0-18 A) and the Fe—Fe* pair at 45169 (separation 0-14 A) are both 
more clearly resolved in the spectrum taken with the diaphragm in position. 

In spite of the reduction in aperture to nearly one-quarter, the diaphragmed 
spectroscope used in conjunction with the solar telescope remains moderately 
fast. Spectra sufficiently dense for the measurement of sight-line velocities 
are obtained with a normal slit-width of 0-04mm on Ilford Rapid Process 
Panchromatic plates in the following exposure times :— 


Region 45900 Dispersion 0-71 mmA-! Exposure I second 
4227 2°23 8 seconds 
4010 2:80 64 seconds 


It 1s the performance and limitations of the spectroscope, masked as described 
above, which are discussed in the remaining two sections of the paper. 

1.3. Temperature control.—Because of the change of refractive index with 
temperature, it is essential that the temperature of the prisms be controlled within 
narrow limits. For the glass of these prisms a change of I deg. C gives a wave- 
length displacement of 0-6A; thus if between the top and the bottom of the 
prisms there were a temperature difference of 0-023 deg. C a line widening would 
be produced comparable with the minimum distance, 0-014 A, which is theoreti- 
cally resolvable at 44320. ‘The only possibility of realizing so high a degree of 
uniformity within the prisms is to maintain them at a temperature constant within 
about these limits ( + 0-023 deg. C), and to trust that the relatively high conductivity 
of their cast-iron cells will smooth out any small residual temperature gradients. 

On the advice of the Building Research Station no heating elements are 
placed inside the cork-lined enclosure which contains the prisms, the objective 
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and their base. Instead the temperature of the air outside the enclosure, that 
is in the tunnel, is maintained constant within 0-1 deg. C by means of two heating 
circuits, each dissipating one kilowatt and each controlled by its own thermostat. 
The thermostats are moderately sensitive mercury thermometers with sealed-in, 
adjustable platinum contacts. Each carries a current of 10 micro-amp; this 
current controls the potential of a grid, and the latter in turn the anode current of 
the valve which operates the on-off switch of the power supply to the heating 
circuits. ‘The thermostats and valve relays are supplied by Baily, Grundy and 
Barrett of Cambridge, and give continuous, trouble-free service. 

The thermostats control the temperature in their immediate vicinity to better 
thano-1deg.C. ‘Toensure that the whole of the tunnel is also at this temperature, 
it is essential to keep the air vigorously stirred. On the floor of the chamber 
containing pier A (Fig. 1) and at diagonally opposite corners are two powerful 
fans which run continuously ; these keep the air outside the cork-lined enclosure 
in rapid turbulent motion. In addition three fans on the roof of the tunnel are 
switched on for a few hours each day to help maintain uniformity of temperature 
between the chamber and the tunnel proper. It is of course essential that the 
moving currents of air should not actually cross the optical path and produce 
“bad seeing”. ‘The optical path from the slit to the front component of the 
objective is therefore completely enclosed by a metal-lined three-ply tube of 
rectangular cross-section, fixed to the walls of the tunnel but flexibly connected 
by a bellows to the opening in the cork-lined enclosure. ‘This metal tube helps 
to maintain uniformity of temperature across the beam, but in addition, and 
before an exposure is made, air may be blown through the tube to eliminate 
temperature gradients. 

The prisms, the objective and their base in the cork-lined enclosure are 
thermally insulated from the air of the tunnel and from the pier on which they 
stand. By aslow process of conduction the inside walls of the cork-lined enclosure 
assume the temperature imposed on the air outside, but with the temperature 
ripple greatly reduced in amplitude. ‘To prevent stratification within the 
enclosure, the air is stirred for a few hours each day (though not when an exposure 
is in progress) by a fan blade driven by a belt from a small motor on the floor of 
the tunnel. ‘The temperature of the air inside the enclosure can be read from a 
sensitive Beckmann thermometer to an estimated 0-001 deg. C, and is found to 
remain reasonably constant. ‘Thus in the 1951 observing season from April 5 
to November 10, when the temperature of the tunnel was continuously main- 
tained at 24-0 deg. C + 0-1 deg. C, the temperature inside the enclosure was read 
135 times at more or less regular intervals. ‘The root mean square deviation of 
a single reading from the mean was + 0-027 deg. C—a constancy close to the 
stipulated requirement. 


2. Characteristics 


2.1. Collimator illumination.—The Cassegrain mirror of the solar telescope 
(13) forms a central stop, the diameter of which is 0-6 that of the main mirror. 
Each point in the solar image is thus at the vertex of a hollow cone of rays, and from 
the geometrical—optical point of view this hollow cone would be expected to 
propagate itself beyond the slit to give an annular illumination of the spectroscope 
objective. Examination of the collimator illumination for a normal slit-width 
of o-04mm, and a height of 3mm (effectively a pin-hole for the 9 m distant 
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objective), shows, however, very nearly uniform illumination with no trace of a 
central shadow. ‘To get a well-developed central shadow the slit-width has in 
fact to be increased to 0-16mm. This result is in agreement with experience 
with stellar spectroscopes used in conjunction with the large Cassegrain reflectors, 
though these latter have, of course, much smaller obscuration ratios than the 
Oxford telescope. The absence of the central shadow for a normal slit-width 
is due to diffraction. 

A simple, and correspondingly artificial, picture of what is occurring can be 
obtained by regarding the mirrors of the solar telescope as self-luminous sources, 
a pin-hole image of which is being formed by the slit on the spectroscope 
objective. If s be the slit width (0-04 mm), 4 the distance of the objective from 
the slit (taken as 8-74 m) and A be the wave-length, then according to Rayleigh 
(15) the diameter of the image on the collimator of a luminous point on the 
telescope mirror will be 

D = 2bX/s = 26-2 cm for A 6000, 


that is, nearly twice the diameter of the collimator itself. ‘The actual distribution 
of intensity in this image, given by the ordinary two-dimensional Fresnel theory, is 
A(x) =const.[{C(v,) — C(v,)}*+ {$(e2) — S(e)}?) (1) 


, a T . I ete 
where C(v) = | cos—v" dv ~- + —sin 
/0 “ “ TU 
peo 
S(v) = | sin — v* dv = 

0 2 
and the approximate expressions are accurate (14) to better than a tenth of 
I per cent forv>6. In eqn. (1) v, and v, are defined by 


ee a Sy Ri AF a ae s Ba 
= a+b” 2 — oe sat" 3 aby ” 


where x is the distance measured along the arc (in the image plane) from the 
centre of the geometrical image of the slit, a is the distance from the luminous 
object to the slit (in our case taken to be 8-83 m), and b( =8-74m), s and A have 
the meanings already assigned. 

Using eqn. (1), A(x) has been calculated for a wave-length of A 6000 for slit- 
widths 0:04 and o-16mm, and normalized. Then if 7(x) be the intensity 
distribution on the collimator as given by geometrical optics, the actual intensity 
distribution will be given by 


- a 


O(x) = | T(E)A(x—£) dé. (2) 


The result of this calculation is shewn in Fig. 2. In so far as this simple treatment 
may be regarded as valid, the figure makes it clear why the shadow visible for 
a slit-width of 0-16mm disappears for a normal slit-width of 0-04 mm. 


>” 


2.2. Focal surface.—-Focal settings of the objective as a function of position 
on the photographic plate are shown in Fig. 3. ‘The abscissae are distances, 
measured in centimetres, from the left-hand edge of the plate; for this origin 
the slit is atx=—3cm. ‘The ordinates are the focal settings, determined by the 
Hartmann method, measured in centimetres and with an arbitrary origin; the 
larger the setting the farther the objective is from the photographic plate. Curve I 
is obtained by backing the objective, in position in the spectroscope, with a plane 
mirror and determining the focal setting for the slit illuminated by the D lines 
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Fic. 2.—Collimator illumination. 
Abscissae : Distances x in centimetres from centre of collimator. 
Ordinates : Illumination in arbitrary units. 
Curve A: Illumination according to geometrical optics. 
Curve B: Illumination from Fresnel theory for slit-width 0:16 mm. 
Curve C : Illumination from Fresnel theory for slit-width 0-04 mm. 
Filled circles are points computed from eqn. (1). 
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0 +4 8 12 16 +-20 cm 
Fic. 3.—Focal settings in centimetres (ordinates) against positions in centimetres in focal plane 
(abscissae). Slit at x 3 cm. 


Curves I and II are monochromatic (A 5900) focal settings for double transmission through 
objective backed by plane mirror (1), and backed by prisms (II). Curve III gives focal settings 
as functions of wave-length. 
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of Na; the slit image is moved from point to point across the plate by slightly 
rotating the mirror and the focal setting determined for each point. Curve I 
thus shows the curvature of field for double transmission through the objective ; 
its radius of curvature is 4-4 m, a satisfactorily large amount. Curve II shows 
the focal settings for the D lines after double transmission through the objective 
and prisms; the D lines are moved across the plate by changing the position of 
minimum deviation of the prisms. The difference between the monochromatic 
curves I and II clearly indicates that the prisms are having some effect on the 
focal properties of the system. Finally curve III gives the focal settings as 
normally found in taking a spectrum; 5900 is roughly central, and the plate 
covers a range from A 6073 at the left to 45725 at the right-hand edge. ‘The only 
wave-length in minimum deviation in this spectrum is that returned back on to 
the slit itself, that is, A6125 approximately. Curves II and III differ because of 
the chromatic properties of the objective, the minimum focus of which occurs 
atA5100. ‘The three curves, it should be noted, have been brought into agreement 
at x =0 by the addition of a small arbitrary constant to the focal settings. 

The unexpected feature of these focal settings is the difference between 
monochromatic curves I (objective backed by mirror) and II (objective backed 
by prisms). Its interpretation I owe to Mr J. W. Perry. ‘The slit is some 6cm 
farther away from the objective than is the photographic plate, so that the beam, as 
it emerges after its first transmission through the objective, is slightly convergent. 
For rays reflected off a plane mirror (curve I), or passing through a system of 
prisms in minimum deviation (extreme left of I1) this convergence is unaltered, 
but for rays passing through prisms away from their position of minimum devia- 
tion (as happens at right of Il) the convergence, as Herman (4) at least suggests, 
must be increased. A trigonometrical ray-trace through the system for the 
D lines on the left of the plate (that is, 45 A away from the position of minimum 
deviation), and on the right of the plate (470 A away from position of minimum 
deviation), predicts a difference in focal setting of 1:07cm. The actual difference, 
as compared with curve I, is 0-53cm. ‘The predicted increase of convergence is 
thus of the right order of magnitude, but it appears probable that the prisms are 
having some additional effect on the focal properties of the objective. This may 
arise because, out of minimum deviation, the objective is vignetted (noticeably 
so if the collimator illumination be examined from the right-hand side of the 
plate), and this vignetting may affect the field curvature which is determined by 
the Petzval sum and the extra axial astigmatism. 

2.3. Horizontal apparatus function.—Following Redman (16) the horizontal 
apparatus function has been found from the measured profile of the Kr line, 
A4319°6. This line, examined under high resolution by a number of workers 
(7), is apparently without structure and has a half-width of the order of 0-008 A. 
Since this is less than the half-width of 0-0126A of the theoretical apparatus 
function (predicted by diffraction theory), it is reasonable to suppose that the 
profile of this line, measured with the spectroscope, will give a good determination 
of the actual apparatus function. It should, however, be emphasized that van 
Cittert (2) has shown that the apparatus function depends upon the mode of 
slit illumination (coherent, non-coherent, or some combination of the two). 
Consequently, while the apparatus function about to be found may be regarded 
as typical of what the spectroscope can do, it is not necessarily applicable to 
spectra obtained under different experimental conditions. 
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The Kr line A 4319-6 was produced in a G.E.C. discharge tube run at 300 V 
d.c. A field lens, a few centimetres in front of the slit, formed an image of the 
small region of discharge on the collimator; the magnification was such that this 
image was slightly larger than the aperture over the collimator. ‘The extra- 
focal patch from the field lens gave a very nearly uniform illumination over the 
height of the slit. The slit was set at the normal or van Cittert width for this 
wave-length region, of 0-025mm. Over the slit was placed the neutral-tint 
wedge (W. iv) used in previous investigations, and one spectrum of Kr A 4319-6 
was taken through this wedge; the wedge was then removed and through a 
neutral-tint filter an exposure of identical duration was made. ‘The first of these 
spectra was thus photometrically self-calibrated under simultaneous and identical 
conditions, a considerable photometric gain; the second gave a measure of the 
flux incident on the wedge as a function of height, as well as another determination 
of line profile. ‘Three backed Ilford Special Rapid plates were taken, with two 
or more spectra per plate, and with exposure times of 4, 5$ and 64 hours per 
spectrum. After a lapse of at least two weeks the plates were brush-developed 
in Parkhurst developer. 

The neutral-tint wedge, W. iv, was re-calibrated in position on the slit of the 
spectroscope. ‘The line 44358, given by a G.E.C. “‘Osira’’ mercury lamp, was 
used as a source with the same mode of collimator and slit illumination described 
above. This lamp can scarcely be regarded as a photometrically steady source, 
but by varying the order in which exposures were made it was possible to convert 
any possible secular change of flux with time into an accidental error. ‘To 
change the flux incident on the slit by a known amount, one or more of four 
neutral filters could be inserted before the field lens. The transmission of 
these filters had been measured at the National Physical Laboratory at six different 
wave-lengths; from measures made there it was also known that the greatest 
angular deviation produced by any filter was 22 seconds of arc, corresponding to 
a displacement on the collimator of under a millimetre. Any error produced 
by this slight displacement was made accidental by changing the orientation of 
the filters in different exposures. ‘The calibration was based on 36 wedge and 
48 filter exposures, wedge and filter exposures on any one plate being of exactly 
identical duration (from 20 to 80 seconds depending upon the plate). From 
microphotometric measures of these spectra the transmission of the wedge as a 
function of height could be found in the usual way; as in previous investigations 
of this particular wedge (10, 11) a strictly linear variation of the logarithm of the 
transmission with height was found with a wedge constant of 


o,=0°275 +0001 mm}, 


This may be compared with the values of 0-268 + 0-003 found in 1922 and 
0-278 + 0:003 found in 1936 at the same wave-length but under very different 
experimental and spectroscopic conditions. Adopting 0-275 for A 4358, the 
value at the Kr line 4319-6 may, from the previously found variation of wedge 
constant with wave-length, be taken as 0-278. 

From this calibration of the wedge and from the spectra of the Kr line taken 
under the conditions already described, it was readily possible to derive profiles 
for the line. ‘The spectra were marked at intervals of 0-4mm of height on the 
line, and microphotometer tracings were made in the wave-length direction at 
these heights with an effective slit width of 0-03 mm, a height of 0-14 mm and a 
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recording magnification of 48-2. From these cross-sections at different heights 
in the wedge spectrum a provisional characteristic curve was constructed, and 
this in turn was corrected in the usual way (12) for non-uniformity of incident 
illumination. From the true characteristic curve thus derived and from the 
microphotometer deflections as a function of linear position on the plate it was 
possible to derive a profile for each measured spectrum. All the relative intensities 
(arbitrarily made 100 for the line centre) were plotted as a function of wave- 
length, the latter found from the dispersion, and a smooth curve was drawn 
through them. From five well-exposed spectra on the three plates five such 
curves were thus available. The values of intensity at a number of specified 
wave-lengths were then read off these curves, and the arithmetic mean taken of 
the five intensities at each wave-length. The root mean square residual from this 
mean was found to average 6-3 per cent—a satisfactorily small quantity when it is 
noted how exceedingly sensitive the intensity is to wave-length. 
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Fic. 4.—Apparatus function of spectroscope. 

Intensity in arbitrary units as a function of wave-length (abscissae). 
Curve A: for masked spectroscope. 

Curve B: with no mask. 

Curve C: diffraction theory. 
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The resulting mean profile is shown in Fig. 4 as a function of wave-length 
(continuous line with observed points). It is well determined except in the wings, 
which are a pure extrapolation, guided to some extent by the wings which have 
been found for the vertical apparatus function. ‘This weakness in the wings is 
the price which has to be paid for self-calibration by exposure through the wedge. 
In accordance with our initial assumption we take this line profile of Kr A 4319-0 
to be the horizontal apparatus function of the spectroscope at this wave-length. 
In so far as the Kr line has a definite width and the photographic plate a finite 
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resolving power this apparatus function errs in the sense of exaggerating the 
width of the core. Its characteristic features are the large half-width of 
0-048 A and the presence of a very strong satellite on the long wave-length side. 
For comparison there is included in Fig. 4 the profile of the Kr line obtained from 
four spectra made in the same way with the undiaphragmed spectroscope. Since 
the half-width of this profile is 0-076 A, nearly 60 per cent greater than that with 
the spectroscope diaphragmed, it is evident how much is gained by always using 
the spectroscope in the latter form. Also for comparison is given the ideal 
diffraction pattern with a half-width of 0-0126A for an infinitely narrow slit, 
assuming the absorption in the prisms to be that given by Eppenstein (3). ‘The 
practical and theoretical resolving powers, corresponding to these patterns, are 
93900 and 303 000 respectively. 

The difference between the actual and the predicted apparatus functions is thus 
considerable. It is due to aberrations in the spectroscope, largely produced by 
the sharp changes of refractive index which, from the Foucault test, are present 
in the glass of the prisms. All solar spectroscopes, whether grating or prismatic, 
show similar, though on the whole less marked, departures from theoretical 
expectations. In this respect they contrast unfavourably with the small 
laboratory spectroscope. It is, however, perhaps worth recalling Redman’s 
remark (16) that it is only in instruments of large focal ratio, such as solar spectro- 
scopes, that the approximate half-width of the diffraction image (focal ratio times 
wave-length) exceeds the minimum sized image (diameter = 0-02 mm) given by 
the ordinary photographic plate. Hence photographic laboratory instruments 
with focal ratios smaller than 20 could depart nearly as widely from theoretical 
performance as the average solar spectroscope without the user being aware of it. 

2.4. Vertical apparatus function.—The vertical apparatus function gives the 
intensity as a function of y (that is, perpendicular to the dispersion) of a point 
source on the spectroscope slit. To find this function, a fine wire of known 
diameter was stretched across and very nearly in contact with the slit. Spectra 
of the Sun, taken at 44320, then showed the shadow of the wire, the intensity 
distribution in which could be found from wedge spectra exposed on the same 
plate with identical exposure time. ‘The spectrophotometric details of the work 
were carried out as described in the preceding sub-section, except that the plates 
were not brush-developed and that a microphotometer magnification of 6-89 
was used. Five plates were available for determining the intensity distribution 
in the shadow of the wire, and from these the mean intensity as a function of y 
was found precisely as described above. The root mean square residual for a 
single determination was 2 per cent, a satisfactorily small error. Systematic errors 
arising on the one hand from imperfect contact of the wire with the slit and 
finite resolving power of the emulsion, which tend to wash out the profile, and 
on the other hand from Eberhard ef'ect which tends to deepen it, have not been 
estimated, but from the shallow, broad shape of the profile should not be large. 

The observed intensity distribution, O(y), in the profile of the shadow is 
given by a2 

A(u) du. (3) 
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Here T(y) is the shape of the shadow in the focal plane in the absence of smearing, 
that is a rectangular profile of width a (=0-057mm) and zero intensity, and 
A(y) is the normalized vertical apparatus function. A(y), thus found from the 
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observed profile of the shadow, is plotted (filled circles) as a function of y in 
Fig. 5; the abscissae are millimetres in the focal plane, counted pesitive in the 
direction towards the bottom of the focal plane. The half-width is 0-167 mm. 
For comparison the diffraction-theory apparatus function, corresponding to the 
masked objective with rectangular aperture of vertical height 3-1cm, is shown 
in the same figure as a dotted curve. This diffraction pattern in the vertical 
direction has a half-width of 0-098 mm as contrasted with the 0-025 mm in the 
horizontal direction, but in spite of this great theoretical width the observed 
vertical apparatus function is still wider. This is a consequence, presumably, of 
non-homogeneity of refractive index in the prisms together with astigmatism 
due to the spectrum not being in exact minimum deviation. 
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Fic. 5.—Vertical apparatus function (normalized). 
Abscissae are millimetres in focal plane, positive values downwards. Inset is \n A(y). 


2.5. General scattered light.—Superposed on the spectrum there is a flux of 
general scattered light. ‘This originates from that white light, admitted by the 
slit but lost to the image-forming beam by reflection and scattering on the various 
optical surfaces, which either directly or after further reflection and scattering 
arrives in the focal plane to produce there a more or less uniform illumination. 
The ratio per unit area in the focal plane of scattered-light flux to flux in the 
spectrum thus increases approximately linearly with the height of the slit and with 
the dispersion (mm A~*) of the spectrum. The ratio also increases towards the 
short wave-lengths where the true absorption of the image-forming beam increases. 

The ratio per unit area of scattered to spectral light at A 4320 is found from 
five backed plates. A small range from one of these plates is shown in Plate 2, 5. 
The four exposures on each plate are spectra of the centre of the solar image 
made with a wide slit (0-075 mm) and with precisely identical exposure times (of 
the order of 3min). ‘The spectrum at the top, with radiation unanalysed for 
colour before entering the slit, is made with a slit height of Ig mm and exposed 
through a diaphragm in the focal plane 33 mm high. ‘The image of the illuminated 
height of the slit is heavily over-exposed; on either side of it, but limited by the 
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height of the diaphragm in the focal plane, can be seen a more or less uniform fog. 
This fog is produced by the general scattered light. The spectrum at the bottom, 
also with radiation unanalysed for colour before entering the slit, is of the same 
part of the solar image, but with the flux per unit area incident on the slit reduced 
by a factor of 0:0137. Thisis achieved by interposing three of the neutral-tint 
filters, calibrated by the National Physical Laboratory, in front of the slit. For 
this spectrum the slit height is 3-2 mm, and the height of the diaphragm in the focal 
plane 9mm. From the roughly comparable densities in this spectrum and in 
the fog on either side of the top spectrum it is evident that the scattered light for 
a slit height of 19 mm is only some 2 per cent of the light inthe spectrum. ‘The two 
remaining spectra on the plate are a wedge and non-wedge spectrum to give the 
characteristic curve, corrected for any non-uniformity of flux incident on the 
wedge. They are taken through a diaphragm in the focal plane 9mm high. 
To keep the superposed scattered light to a minimum in these two spectra the 
radiation from the Sun, before entering the slit, passes through a Wratten No. 39 
filter. 














1 i | 
0 Spectrum 0 +2 +4 +6mm 
Fic. 6.—Values of J for a typical scattered-light spectrum (Plate 2, 5). 





With magnification 6-89, microphotometer tracings are made across the plate 
in the y-direction, that is, perpendicular to the dispersion. From the measured 
deflections and the characteristic curve obtained from the wedge spectrum 
(corrected for non-uniformity and assumed free from scattered light) the intensity 
J (in units of flux per unit area in the solar spectrum) is found as a function of 
y from the fog on either side of the top spectrum made with the 19-mm high slit. 
The values of J thus found for a typical plate (that illustrated in Plate 2, 5) 
are shown in Fig. 6. If S(A) is the general scattered-light flux per unit area for 
a slit hmm high, and / the flux per unit area in the adjacent spectrum (found 
from the spectrum taken with a slit 3-2 mm high), then 


J S(t9)_ S19) S( 





“ 19) , 
“7+ S(32) 1+4809) ” “T =J(I+iJ), 


on the assumption that the general scattered light varies linearly with slit height. 
Taking J to be of the order of 0-02, then the J-curve represents the required 
S(19)/J with an accuracy of 0-3 per cent—an accuracy greater than that of the 
spectrophotometry, where a characteristic curve for 4320 has been necessarily 
used for the white-light flux J. 
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Taking the J-curve then as the required ratio of scattered-light flux to flux 
in the spectrum, its variation with y has to be explained. As can be seen from 
Fig. 6, J decreases rapidly with increasing distance from the spectrum, but at 
the top of the plate (neg. y) commences to increase again beyond y= —3 mm. 
This increase arises because we are leaving the umbra of the shadow cast by the 
diaphragm for the virtual image of the 19-mm high slit, and are entering the 
penumbra; in fact, if the eye is placed in the focal plane beyond y = — 10mm the 
virtual image can be seen escaping beyond the edge of the diaphragm. For 
positive values of y, that is, towards the bottom of the focal plane, we are always 
in the umbra and are getting the true scattered light which must be incident on 
the spectrum. 

The increase of J towards the spectrum is due to the wings of the vertical 
apparatus function. ‘That this is so can be seen from Plate 2, 5, where it will 
be observed that the Fraunhofer lines extend well beyond the limits of the 
spectrum fixed by the height of the slit. ‘The flux ratio originating in this way 
towards the bottom of the plate is 

SP = [° ean) d= |" AQw) a 

where i(y) is the flux ratio per unit area in the spectrum (i(y)=1 for y<o, and 
=0 for y>o), and A(y) is the normalized vertical apparatus function. We 
then have O(y) _ S(19) 

— es 
and if we assume that S(19)/J is given by the asymptotic value to which J tends 
for large values of y, we can find O(y)/J and therefore the wings of the vertical 
apparatus function. The results for In A(y) are shown as an insert in Fig. 5, 
where the open circles are the mean values found from the five scattered-light 
plates and the filled circles are the values found in the preceding sub-section. 
Though the two sets of values do not overlap they do at least fit at their junction ; 
since both determinations give normalized values of A(y), that is, absolute values, 
this fit is quite unforced and may be regarded as satisfactory. ‘The dotted curve 
is tangent to the maxima of the theoretical diffraction pattern. For positive y 
the wing of A(y) is nearly 30 per cent greater than that given by theory; for 
negative y, where the wing cannot be determined because of the penumbra of 
the shadow of the diaphragm, there is at least an indication that A(y) is 
approaching theory. 

The asymptotic values of J(y) for y large, that is, the required ratio of general 
scattered-light flux to flux in the spectrum, are given in Table I for each plate. 
It will be noted that the mean ratio of scattered light for the two panchromatic 
plates is 1-54 per cent, and for the two ordinary plates 1-15 per cent. ‘This is the 
anticipated result. ‘The general scatered light is made up of all wave-lengths 
which enter the slit (in our case the whole spectrum of the Sun as transmitted by 
the Earth’s atmosphere), but the photographic plate can only record that part of 
this wave-length range to which the emulsion is sensitive. Hence a greater 
scattered-light ratio will be given by the panchromatic than by the ordinary 
emulsion. In view of the great height of the slit, the dispersion of nearly 
2mmaA_! and the short wave-length of 44320 where the determination is made, 
a scattered-light ratio of 1-54 per cent may be regarded as satisfactory for an 
auto-collimating spectroscope. By a preliminary analysis with an ordinary 
colour filter of the radiation entering the slit and by using a slit only a few 
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millimetres high, it is clear that even at A 4320 the general scattered light can be 
reduced to negligible proportions. 


3. Limitations 
The technical problems of solar spectroscopy may be classed into profile 
determination, measurement of line displacements and discrimination of spectra 
from different parts of the solar image. Given the mechanical and thermal 
stability, and the freedom from scattered light of the Oxford spectroscope, the 
limits within which reliable results can be obtained in these problems are fixed 
by the horizontal and vertical apparatus functions. 
Taste | 
Flux Ratio of Scattered Light/Spectrum at A 4320 
Date : Slit a Exp. 
Height Width 


1945 m s 
4 mm mm 


S(ig)/I Emulsion 


June 19 0'075 3 03 0°0146 Ilf. R. Proc. Pan. 
June 19 0075 3 00 00163 Ilf. R. Proc. Pan. 
July 19 0°075 3 00 o'or2I lif. Sp. Rapid 
July 5 19 0°075 3 00 00120 llf. Sp. Rapid 
July 5 19 0°075 3 00 O°0104 Ilf. Sp. Rapid 


3.1. Profile determination.—The observed intensity at nominal wave-length 
A in a line profile is 


O(A)= | TA)AA-X) aX, (4) 


where 7(A) is the true shape of the profile and A(A) is the normalized horizontal 
apparatus function. ‘The value of the required 7(A) derived from this equation 
is very sensitive to errors in O(A) and A(A), and becomes the more sensitive the 
more nearly 7(A) approaches in half-width the horizontal apparatus function; 
compare van de Hulst (6). The horizontal apparatus function therefore fixes 
a lower limit to the width of the Fraunhofer line for which the profile can be found. 
Thus Redman (16), with an apparatus function of half-width approximately 
0-145 A, has confined his determination of central intensities to lines of Rowland 
number 8 and greater, and Houtgast(5), with the Potsdam spectroscope with 
a remarkably sharp apparatus function of half-width 0-029 A, to lines of 7 and 
greater. ‘Io Redman’s rule (17) that “existing instruments are inadequate to 
deal with lines fainter than about Rowland number 5”’ the Oxford spectroscope 
with an apparatus function of half-width 0-048 A is no exception. 

If the spectroscope be used as a monochromator for an interferometer, then 
it is the apparatus function of the latter which appears in the integral equation (4). 
Since the resolving power of the interferometer can be very high and the width 
of its apparatus function correspondingly small, it might be thought that no 
limit exists to the practical determination of absorption line profiles with an inter- 
ferometer combined with a spectroscope, ‘l'reanor (18), however, has shown 
that, owing to the periodicity of the apparatus function of the interferometer, 
increasing amounts of false light are superposed on the fringe system as the 
resolving power of the interferometer is increased and with it the range is decreased. 
Thus using ‘l'reanor’s exact equation (9) we find for the Oxford spectroscope 
false light of 0-24 per cent for a Fabry—Perot interferometer of range 0-36 A and 
resolving power 500 000, and false light of 1-50 per cent for a range of 0-18 A and 
resolving power of 10°. ‘I‘hese are small amounts and since they are accurately 
calculable, provided the apparatus function is known for the circumstances under 

14* 
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which the spectroscope is being used (wave-length, slit-width and nature of its 
illumination), there is reason for hope that profile determination below Redman’s 
limit may be possible. 

3.2. Line displacement.—Because of the marked asymmetry in the horizontal 
apparatus function the observed profiles of solar absorption lines will also 
necessarily be asymmetrical. ‘To show this effect “‘true”’ profiles for Fraunhofer 
lines have been computed for a purely scattering atmosphere (above a bounded 
photosphere) at a temperature of 6000deg. K and with ten times the classical 
damping constant, and from these the observed profiles have been computed by 
eqn. (4). Calculations have been carried out for four lines with equivalent 
widths of 30, 69, 125 and 327 mA, corresponding roughly to Rowland numbers 
0,,2,4andg. Fig. 7 shows the “true” and observed profiles for the first and 
second of these lines; the asymmetry of the observed profile is easily visible. 
The centres of intensity and the centroids, the latter found from numerical inte- 
gration, are given in Table II. ‘The centres of intensity of the observed lines are 
displaced to the red with respect to the “true” lines by an amount which varies 
nearly linearly as the logarithm of the equivalent width of the line. The centroids 
of the lines, on the other hand, all show the same displacement within the errors 
of computation; this displacement, as is mathematically obvious from eqn. (4), 
for a symmetrical, true profile, is simply the displacement of the centroid of the 
apparatus function from its centre of intensity. 


Tas_e IT 
Displacements of Observed Lines 
W Centre Int, Centroid 
0030 A 0'003 A o-018+o0'001 A 
0069 "010 0019+ 0°002 
O°125 O'oIs 0°018+0°005 
0°327 0°022 0°029 + 0'028 

If the measurer of a spectrum were always presented with intensity profiles 
of the observed lines, and if he set his cross-wire on the centroid of a line, then the 
asymmetry of the apparatus function would only Jead to a change of origin in 
wave-length. This would have no effect on the determination of sight-line 
velocities, or on the measurement of any other displacement where true, sym- 
metrical profiles are bodily displaced. Measurement, however, is never so simple 
a matter. Quite apart from the transformation of the observed intensity profile by 
the photographic emulsion and of its appearance by Fechner’s law, it is unknown 
where in this transformed profile the experienced measurer sets his cross-wire. 
The art of measurement is learned by finding a method of setting which leads to 
the greatest consistency, and the setting which gives this result may well be 
different for lines of different equivalent width and indeed for different measurers. 
It is therefore almost certain that the setting on the asymmetrical profiles produced 
by the Oxford spectroscope will be displaced by an amount which is a function of 
the equivalent width. ‘This would for example be the case if the setting were 
always made on the centre of intensity of the observed profile. 

The only way to overcome this difficulty is to make differential measures, 
that is, the position of a particular line must always be measured from an origin 
given by the same line photographed under standard conditions where its wave- 
lengthisknown. ‘Thus no matter how the setting is made, the measured displace- 
ment is correct and is independent of the equivalent width of the line. Standard 
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conditions are provided by a spectrum of the centre of the Sun, and the differential 
measurement is most conveniently made by using a Hartmann spectrocomparator 
or Evershed’s positive-on-negative method. 

3.3. Discrimination of spectra.—Fine details on the solar surface, more or less 
imperfectly imaged on the slit (as a result of seeing and the telescope) give spectra 
still further obliterated by the vertical apparatus function. A convenient test 
object of this phenomenon is provided by solar granulation. According to 
Keenan (9), granules have a diameter of approximately 1-1 seconds of arc, a 
brightness ratio to intergranular regions of 1-15, and occupy about 26 per cent 
of the area of the solar surface. Assuming the granules to be evenly spaced, a slit 
0-41 seconds of arc in width (corresponding to a normal Oxford slit-width of 
0:04 mm) will cut centrally through a granule once every 8-9 seconds of arc, and 
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Fic. 7.—‘* True’ (dashed) and observed (continuous) profiles of Fraunhofer lines with 
0°030 (left) and =0-069 (right). Circles are values computed from eqn. (4). 





in the same interval will cut off-centre through a number of others. The illu- 
mination, observable just behind the slit, produced by a somewhat simplified 
form of granulation of this kind and imaged in perfect seeing by a perfect 
telescope is shown in Fig. 8. The variation in the heights (proportional to the 
flux per unit length through the slit) of the granules is due in part to the slit not 
cutting the granules centrally and in part to an assumed variation of brightness 
due to the finite life of the granules. 

The spectra of these idealized granules tends to be obliterated by the vertical 
apparatus function. The resulting intensity distribution, computed from the 
vertical apparatus function and an equation similar to (4), is shown as a smooth 
curve in Fig. 8; this is simply the intensity distribution of the streaks produced 
in the spectrum by the granules. ‘The amount of obliteration is considerable, 
but it seems reasonable to suppose that spectroscopic details of granules, affected 
in a known way by the intergranular regions, could be obtained for the three or 
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four brightest granules shown in the 4-mm height of slit in Fig. 8. Observations 
of granular spectra during the last four years show indeed much the same intensity 
distribution as has been theoretically calculated. Solar detail, very much finer 
than granulation, even if it were adequately imaged on the slit, would however 
be so obliterated by the vertical apparatus function as to be nearly invisible. 
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0 0-4 0-8 1-2 1-6 2:0 2:4 2:8 3-2 36 mm 
Fic. 8.—Ideal granulation (rectangular profiles) as reproduced by spectroscope (continuous 
curve). The circles are intensities computed from integral equation. Ordinates are intensities 
and abscissae are millimetres in focal plane perpendicular to direction of dispersi 


Help has been received at all stages in the design and construction of the 
spectroscope. I wish to mention specifically that received from the late Mr R. M. 


Abraham of Casella’s, to whom is largely due the success of the mechanical 
parts, from Mr J. W. Perry of Hilger’s, who designed the objective and suggested 
the interpretation of its focal properties and finally from Mr G. H. Neale, who 
supervised the construction of the tunnel, did all the joinery and assisted in the 
erection of the instrument. 


University Observatory, 
Oxford : 
1952 January 3. 
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ON SPHERICALLY SYMMETRICAL ACCRETION 


H. Bondi 


(Received 1951 October 3) 


Summary 


The special accretion problem is investigated in which the motion is 
steady and spherically symmetrical, the gas being at rest at infinity. ‘The 
pressure is taken to be proportional to a power of the density. It is found 
that the accretion rate is proportional to the square of the mass of the star 
and to the density of the gas at infinity, and varies inversely with the cube 
of the velocity of sound in the gas at infinity. The factor of proportionality 
is not determined by the steady-state equations, though it is confined within 
certain limits. Arguments are given suggesting that the case physically 
most likely to occur is that with the maximum rate of accretion. 





1. ‘he importance of the accretion of interstellar gas by stars has been 
recognized since the work of Hoyle and Lyttleton (1, 2, 3). In their work the 
problem, later investigated in detail by Bondi and Hoyle (4), was that in which 
the rate of accretion was limited principally by the relative motion of the star 
and the gas cloud, the effects of pressure being considered negligible in comparison 
with the dynamical effects. The result derived in these papers was that 


dM/dt =27«a(GM)*V-*p,,, (1) 


where M is the mass of the star, dM/dt is the rate of accretion, p,, the density 
of the gas cloud far from the star, V is the relative velocity of the star and the 
distant (undisturbed) parts of the cloud, G is the constant of gravitation, and « 
is a numerical constant which was first estimated to be equal to 2. Later work 
(4) showed that the steady-state equations did not determine «, although it seemed 
likely that it should always be between 1 and 2. It was also shown that if the 
star entered a cloud of uniform density with a plane boundary, « settled down to 
a value near 1°25. 

In all this work pressure effects were neglected, the argument being that any 
heat generated would be radiated away rapidly, so that the temperature of the gas 
was always very low. Considerable mathematical simplification is introduced 
by this assumption, and it was shown that it was likely to be satisfied in most 
cases of astrophysical interest (3). ‘The mathematical difficulties of the more 
general problem, in which both dynamical and pressure effects are considered, 
seem insuperable at present. However, the extreme case of negligible dynamical 
effects is again far simpler, and will be discussed in this paper. It may reasonably 
be expected that the case discussed here together with the case discussed previously 
bracket the complete problem. 

2. The problem to be discussed may be defined as follows: 

A star of mass M is at rest in an infinite cloud of gas, which at infinity is also 
at rest and of uniform density p,, and pressure p,. ‘The motion of the gas is 
spherically symmetrical and steady, the increase in mass of the star being ignored 
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so that the field of force is unchanging. ‘The pressure p and density p are related 
everywhere by 


Pipa =(p Po)” (2) 


where y is a constant satisfying I <y < §. 

With a suitable choice of y, equation (2) is equivalent to the physical condition 
that no heat is radiated or conducted away. Hence the solution should provide 
the most complete contrast possible with the problem previously investigated. 
The equations governing the problem are easily set up. If we take r to be the 
radial coordinate and wv the inward velocity of the gas, the equation of continuity is 


4nr*ov = constant = A (say), (3) 


where A is the accretion rate. 
Bernoulli’s equation is 


2 Pp 7M 
pi +{ 2 .. = =constant (=0). (4) 


oa 
The constant is readily seen to vanish by virtue of the boundary conditions at 
infinity. Combining (2) and (4) we have 


2 " yl GM 
+o bel (2) -1]-—. (5) 
2 Y¥—-I Pao Pow r 


Equations (3) and (5) are two equations for the two variables v and p in terms of 7, 
the distance from the centre of the star. 

The equations may be made non-dimensional by the appropriate use of the 
velocity of sound in the gas at infinity, which as usual we denote by c. By the 
well-known formula 

C=yp. Po (6) 
Let us introduce non-dimensional variables, x, y, z, to replace r, v, p, respectively, 
as follows: 
r=xGM/c*, 
v=ye, (7) 
P= Pa 
Then (3) and (5) take the non-dimensional form 
x*yz =A, (8) 
”1— 1)(y—1)=1/x, (9) 
where A is given by 
A= 4rX( G M)Pc 3p, . (10) 
Accordingly A is the non-dimensional parameter determining the accretion rate. 
It plays the same role as « in equation (1). It will also be observed that the 
relative velocity V of equation (1) has been replaced by c¢ in (10). 

3. The explicit solution of equations (8) and (9) for general y is possible not 

in terms of the variables y and z but only if an auxiliary variable depending only 


on y?/z’1 is introduced. It is particularly interesting that mathematical 
requirements lead to the introduction of this variable, since 


u=yz {y¥—1)/2 (11) 
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has immediate physical significance as the ratio of the local bulk velocity v of 
the gas to the local velocity of sound (yp/p)"*. Substituting (11) into (8) and 
solving for y and z we have 


y =u2-1(Q/x2)-WO+0), (12) 
2 =(A/x?u)?-0, (13) 


Then (9) becomes 
A \20-Dv+D I A \2-Div+) I 
a+ { & RY fat Ss 
iu (5) * y-I (=) x y-1I° (14) 
Rearranging the terms and multiplying by (x?/A)*”-%°'+» we find that (14) 
takes the form 
fu) = )\-*- Dr 4 Yo(x), (15) 


where 


? “vs 


= / 


Alu) = fut +D + oe y- 2 +) = yA +) ‘ 4 ‘ (16) 
3 y—1 


( ) = xD +1) E . = abinctes (5—3Y)(Y +1) (1 ) 
gx) = . a | = “ape +X ° 7 

A study of the functions f and g serves to determine u as a function of A and x. 
The variables y and z are then readily found by (12) and (13). 

4. We shall first assume that 1 <y<§. The two limiting cases y=I, 
y= will be examined later. With this assumption both f and g are each the 
sum of a positive and negative power of the respective variables and hence each 
of them has a minimum. The minimum of f(u) occurs for u=u,,=1 and is of 
value 

f= Uy + )/(y—1) =fm (say). 
The minimum of g(x) occurs for x =x,, = }(5 —3y) and is of value 

g= - > [i(5—3y)} om =g,, (say). 

4y—I1 
In our problem x varies between infinity and the value corresponding to the surface 
of the star. This last value is very small indeed. As an example, if the star is 
taken to be the Sun, and c= 1 km/s corresponding to a gas temperature at infinity 
of nearly 3 000 deg. K, the surface value of x is only 5 x 10-®. Even for a red 
giant the surface value of x would be less than 10~? unless the temperature of the 
gas at infinity were quite improbably high (more than, say, 5 x 10*deg. K). 
Accordingly x will attain the value x,, in the physically significant interval. 
Hence, somewhere in that interval, the right-hand side of (15) will reach a value 
as low as A-*”-DIO +g, But the lowest value f can reach is f,,. Hence A cannot 
exceed A, where 


oe (¥+)D/Av-)) 1 \(+bD/20/-1)) 5 sy (5-3yy&y-1) : 
A. ={— — ‘ (18) 
Fea \e 4 


Our first result is therefore that the accretion rate A cannot exceed the value 
47A,(GM)c-p,.. (19) 
Table I gives the value of A, for a few values of y. 
TABLE I 


y I 1°2 14={ 1°5 8 
re pe 1-12 .. }(0°7)-*5==0°872... 0°625 0*500 0°250 
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5. In order to obtain a more detailed picture it is necessary to discuss f(u) 
and g(x) more fully, and to take account of the boundary conditions of the problem. 
Since, at infinity, v vanishes but p and p tend to finite limits, it follows that u tends 
to zero there. 

It is now easy to consider the problem graphically. Fig. 1 shows f and g 
drawn for the typical case y=4. The resulting variation of u as a function of x 
is shown in Fig. 2 for the cases (i) A= }A,, (ii) A=A, and (iii) A= 4A,. It will be seen 
that the boundary condition at infinity implies that u is very small near infinity 
(points beyond A on the graphs). For A<A,, as x diminishes, u rises gradually 
to a maximum (B) and diminishes to zero (C) as x tends to zero. ‘The closer 
A approaches A, from below, the sharper the maximum B. No part of the curve 
A’B’C’ (on which uw is very large both at infinity and near x=0) is of physical 
significance, since u given by this curve does not satisfy the boundary condition 
at infinity. No jump from the curve ABC to A’B’C’ is possible, since this would 
imply an infinite acceleration. Along ABC the variable u is always less than 
unity, so that the motion is subsonic. Along A’ B’C’ the value of u exceeds unity. 

The case A=A, is quite different. For in this the curves have contact at 
B=B’. Coming from A, the physically significant curve can continue either 
toCortoC’. Inthe first alternative, the curve is the limiting form of the curves 
for A<A,. The curve has a discontinuous tangent at B and there is hence a 
finite jump in the acceleration. Although this is perhaps physically not very 
plausible, there does not seem to be any argument disallowing this motion altogether. 
It may be significant that at B the value of u is unity, so that the bulk velocity 
equals the velocity of sound. 

‘The curve ABC’ is perfectly smooth and monotonic. For x>x, the motion 
is subsonic, while for x<x, the motion is supersonic. ‘The system is in a state 
quite different from any state possible for A<A,. 

If A exceeds A, then the pattern of the curves changes as indicated, and no 
solution is possible. 

We see hence that there are two quite different types of motion. Type I 
exists for A<A,. ‘The motion is everywhere subsonic (except at x= x,, if A=A,), 
and u has a single maximum which is less than or equal to unity. The bulk 
velocity v has a maximum if y<# but not if y>#. (This follows from a simple 
consideration of y’ in terms of u.) For y<# the velocity v tends to zero as r 
tends to zero, for y = } it tends to a limit (equal to }cA), and for y># it tends to 
infinity. The density is always a monotonic function of the radius. Type II 
exists only if A=A,. In this case u, v, p are all monotonic functions of the radius. 

‘The special case A=0 may be briefly referred to here. In this case the gas 
is at rest, forming a tenuous continuation of the star. Since y =0 we have by (9) 


2 1=1+(y—1)/x. (20) 


Figs. 3, 4 and 5 show y and z as functions of x for y=1, 7 and § respectively. 
Three values of A are taken in each case, namely A=A,, A= 4A, and the case A=0 
just referred to. It will be seen from the figures that z does not depend very 
critically on A, varying only slightly between the extreme cases A=0 and A=A,, 
especially in the case of the higher y values. 
6. It remains to discuss the two limiting cases y= 1 and y= % respectively. 
If y=1 equation (g) becomes 


by? +Inz=1/x. (9’) 
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Fic. 1.—f(u) and g(x) for y=}. 






































Fic. 2.—u as function of x for y=}. 
(i) A=; 
(ii) A=; 
(iii) A=4A¢. 
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It is easy to eliminate z between (8) and (9’) since, by (11), y=w in the present 
case. Accordingly 


sy?—Iny = —InA+(1/x+21nx). (14’) 
The minimum of the left-hand side occurs for y=1 and equals 4, while the 
minimum of the bracket on the right-hand side occurs for x=4 and equals 
2—2In2. Accordingly 
A,=}fe?? =1-120..., InA,=0-1138.... (18’) 
This is also the limit of expression (18) as y—>1, so that A, is continuous at y =I. 
Fig. 4 represents this case. 





(iv)and (v) 
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Fic. 3.—y and 2 as functions of x for y=}. 
(i) z (A=o); 
(ii) z (A=, Type 1); 
(iii) z (A=, Type ID; 
(iv) y (A=, Type ID); 
(v) y (A=, Type 1); 
(vi) y (A= 42¢). 





The case y =§ is an even simpler extension of the ordinary case. The only 
real change is that x,, (which decreases monotonically as y increases) now equals 
zero. Accordingly there is now no difference between Type I and Type II 
motions for A=A, =} (Fig. 5). 
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It can be seen that u, y and 2, considered as functions of x and y, are continuous 
in I <y<§. 

7. The final question that must be considered is what determines the value 
of Ain any actual case. It has been seen that the steady-state equations possess a 
solution whenever 0<A<A,. The particular value of A actually occurring must 
therefore be determined by other considerations. ‘This is analogous to the 











4 
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Fic. 4.—y and 2z as functions of x for y=1. 
(i) z(A=0); 
(ii) z (A= Ae, Type 1); 
(iii) 2 (A=A¢, Type I); 
(iv) y (A=, Type ID; 
(v) y (A=A-, Type I); 
(vi) y (A= }A¢). 


velocity-limited case of accretion where « is not determined by the steady-state 
equations (4). ‘The method used there to determine a specific value of « was to 
consider the case in which the star entered a cloud of gas with a plane boundary. 
The use of such a model would lead to very great mathematical difficulties in 
the present case, since the conditions of time-independence and spherical 
symmetry would have to be dropped simultaneously. 
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The boundary conditions at the surface of the star do not seem to be of help 
in the problem, since the star will swallow up any material falling into it without 
imposing any real conditions on the velocity, density, or pressure of the incoming 
materials. 
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Fic. 5.—y and z as functions of x for y=%. 
(i) z (A=o); 
(ii) 2 (A=); 
(iii) y (A=e); 
(iv) y (A=). 

There remains the possibility of investigating the stability of the system with 
respect to small disturbances. Even if only spherically symmetrical perturbations 
are admitted, a partial differential equation of considerable complexity results. 
It is easily seen from it that disturbances are in part propagated with the velocity 
of sound relatively to the material, but the nature of the part that remains behind 
is not easily found. However, this may be a possible method of approach to 
the problem. 

There is yet another possibility of investigating the stability of the system, 
and that is by comparing the energy of the system in its various states. ‘The state 
with the lowest energy would then be expected to be the only stable state. Owing 
to the fact that our system is not isolated (the star itself not being considered part 
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of the system) the validity of this approach is not quite assured, and similarly 
doubts may arise owing to the infinite extent of the system. Nevertheless, since 
the comparison of the energy of the system in every spherical shell leads tothe same 
result, it seems very likely that the method gives the correct answer. 

The energy of the gas per unit mass is constant by virtue of our assumptions 
(cf. equation (4)). Accordingly a comparison of the densities is all that is involved. 
Consider now z as function of A for fixed x. ‘Then it may be seen from equations 
(13), (15) and (16) that for uw less than unity z decreases as A increases. Accord- 
ingly the energy of the system in every spherical shell is lower for the Type I 
state with A=A, than for any other Type I state. Comparing now, for A=A,, 
the Type I and the ‘Type II states, it is immediately seen that, for x >x,,, the 
densities are the same, but that for x<.x,, the density (and hence the energy) is 
lower in the Type II state. 

Accordingly the system has, in the sense described, the lowest energy in the 
Type II state, and we may expect to find a natural system in this state with A=A,. 
If y=, the difference between the Type II state and the Type I state with 
A=A, disappears, and we would expect to find the system in this joint state. 
The result that the Type II state is the one most likely to be realized is very 
satisfactory, since the behaviour of all the functions is most uniform and smooth 
in this state. ‘The result is also in agreement with the intuitive idea that, since 
there is nothing to stop the process of accretion, it takes place at the greatest 
possible rate, i.e. with A=A,. 

8. The two cases of accretion that have been examined so far may be called 
velocity-limited and temperature-limited respectively. ‘The intermediate range 
of cases presents far greater difficulties. However, it may be possible to con- 
jecture what the result is in the following way. 

In the velocity-limited case of accretion the accretion rate A is given (4) by 


A=2:57(GM)*V-p,,, (21) 


3 


while for, say, y = in the temperature-limited case the result is 


A =2n(GM)?c~p,,. (22) 
If we therefore write down the formula 
A==2n(GM)V2+2)-32p,,, 


it seems likely that it represents the order of magnitude of the accretion rate in 
the intermediate case, in which a star of mass M moves with relative velocity V 
in a uniform cloud of gas, in which the undisturbed density and the velocity of 
sound have the values p,, and c respectively. ‘This formula, in agreement with 
intuitive ideas, suggests that if c exceeds V, temperature (pressure) imposes the 
chief limitation on the rate of accretion; whereas if V exceeds c, dynamical 
limitations are of greater importance. 

The limitations due to pressure have probably been somewhat overestimated 
in this work. For if the cloud is able to radiate away some of the heat of com- 
pression then the adiabatic law will not apply, the pressure near the star will 
be diminished, and the accretion rate somewhat increased. How large this 
effect will be depends on the composition of the cloud. If there is a high pro- 
portion of constituents (such as hydrogen molecules (2)) that easily radiate at 
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moderate temperatures, then the effect will be appreciable. In this case the 
effective value of y will be closer to unity than to the standard value for the gas in 
question. If y equals unity the process is isothermal. 

The work of the present paper, together with previous work, is likely to give 
a fair estimate of the order of magnitude of accretion in all cases of physical 
interest. Further progress in this field will probably require the consideration 
of non-steady states. 


Trinity College, 
Cambridge : 
1951 October 2. 
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Summary 

It is shown that the retarding force upon a star exerted by interstellar 
matter through which it moves is effectively the same as that previously 
evaluated by Bondi and Hoyle, whether or not accretion takes place in the 
way they consider, so long as the interstellar matter has the general properties 
they assume. It is, however, suggested that the actual accretion may occur 
in the way considered without the operation of the rest of the mechanism 
envisaged by Bondi and Hoyle. If there is a density-gradient transverse 
to the path of the star, there is a transverse force, and this is evaluated. In 
this case any accreted or captured material must possess angular momentum 
about the centre of the star, and this also is evaluated. Some of the calculations 
illustrate the power of the vector methods in dynamics that were expounded 
by the late Professor E. A. Milne. 





1. Introduction.—The theory of accretion by stars according to the mechanism 
originally investigated by Hoyle and Lyttleton* has been elaborated by Bondi 


and Hoyle.t ‘The possible importance of this process in regard to problems of 
stellar evolution is now widely appreciated. Actually there are two fairly distinct 
fields in which its application is important, though naturally the two effects of the 
mechanism are or may be produced simultaneously. ‘lhe first application is 
to the evolution of an individual star: here all that is relevant about the mechanism 
is the possibility it affords of supplying the star with fresh material. ‘The second 
application is to stellar dynamics and so to the evolution of stellar systems, the 
relevant factor here being the retarding force produced by the mechanism upon 
the star (with its possibly increasing mass). 

The latter application was briefly discussed by Bondi and Hoyle in the paper 
quoted, where they evaluate the retarding force and solve the resulting equation 
of motion. 

In the present paper we show that the force is effectively the same as that 
evaluated by Bondi and Hoyle whether or not actual accretion takes place in the 
way they consider. We evaluate the force when the density of the interstellar 
material is non-uniform, in which case the force has a component transverse to 
the path of the star. If accretion occurs in this case, we show that there is also 
an effect upon the angular momentum of the star, though its astrophysical signi- 
ficance requires further investigation. 

2. Particle-orbits—We require certain properties of a hyperbolic orbit of 
a particle P moving under an inverse square law of force acting towards a centre O. 
Some of them could be quoted from other sources. However, it probably takes 
no more space to derive these, together with the other formulae required, by the 

* F, Hoyle and R. A. Lyttleton, Proc. Camb. Phil. Soc., 35, 495, 1939; 36, 424, 1940. 
+ H. Bondi and F. Hoyle, M.N., 104, 273, 1944. 
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methods expounded by the late Professor E. A. Milne.* It is an opportunity 
to illustrate the way in which these vector methods lead directly from the equation 
of motion to results expressed in terms of the desired parameters. They avoid 
any necessity for introducing auxiliary parameters, such as the eccentricity, which 
are not needed in the final results. 


--> 

We write OP=r and r=|r|, with a similar use of the italic letter for the 
modulus in the case of other vectors. The force of attraction is taken to be 
p/r* per unit mass. Then we have as first integrals of the equation of motion 
(Milne, op. cit., Section 269 (2), (3), (4)) 

rxfr=h, (2.1) 
7? —p/r= W, (2.2) 
hxr= —pr/r+B, (2.3) 
where h, B, W are constants of integration. 

Particle projected from infinity.—We require to consider the case in which P 
is projected “from infinity’’ with velocity U, the perpendicular from O to the 
line of projection being s. We write 

U=Ui, S=sp, q=ixp, (2.4) 
where i, p, q are fixed orthogonal unit vectors. 





Fic. 1. 
As initial conditions we have therefore 
r=U, r/r=—i, j : (2.5) 
Using these in (2.1)-(2.3) we obtain 
h=sxU= —sUq, (2.6) 
W=}U?, (2.7) 
B= —pi-—sU"p. (2.8) 
Asymptotic velocity.—Now let F.,.=u,,i+v.p be the final velocity when P 
is receding ‘“‘to infinity”. So |f,|=U and hence, when f=F,,, we have 
r/r=f,/U. Then (2.3) gives 
h xf. = —pr,,/U+B. 
This yields 


hx(hxr,)=—phxer,/U+hxB 


re ~~ Pp B+hxB. 





* E. A. Milne, Vectorial Mechanics, Sections 265-274, London, 1948. 
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Noting that the left-hand member reduces to —h?r,,, (2.9) gives 
2 
te £hd oo Oe 
(: + ir) te 78 hx B. 
Using (2.6), (2.8) this leads to 
2U 
I+s* U4 . 
2sU% 
p(t + sus pe”) 3 
Point C.—We require to find the characteristics of the point C of the orbit 
> 
such that OC is parallel toU. At C let 


r=xi, r=F,=ui+ vp. (2.12) 


Un = (2.10) 


Vp = 


(2.11) 


Then from (2.1), (2.3) 
xixf,=h, (2.13) 
hxf,= —pi+B. (2.14) 
If we multiply (2.14) by i. and note, using (2.13), 
hxr,.i=h.r, x i= —h*/x, 
we find 
x=s*U?/2u. (2.15) 
From (2.14), on multiplying by ix and px, respectively, we find further that 
u=U, (2.16) 
v= —2p/sU. (2.17) 
Particle projected from general point—asymptotic velocity.—We shall require 
the asymptotic velocity rf, of a particle projected with initial velocity fy from a 
general point ry. Let unit orthogonal vectors i, j, k be taken such that 


and write Fr, =u,i1+2,j, |r,|[=a. (2.19) 
Then from (2.1), (2.3), (2.18) 
h = xvok, (2.20) 
B = (pu — xv9”)i + xUgUoj. (2.21) 
Since r/r=f,/w for the asymptotic motion, (2.3) gives 


h xf, = —pr,/w+B. 
This yields 
h x (hx #,) = —(p/w)(h x F,)+hxB 
= —(y/w)(—pr,/w+B)+hxB 
(h? + ?/w?)r, =(./w)B—h x B. (2.22) 
Using (2.20), (2.21) we then find 
Uy =("? — Uxv9? + x*ugry?>w)w/ (pL? + x°v_"2*), (2.23) 
y= (yxugry — wxvyw + x*vy80)e0/(u? + x2 40?). (2.24) 
Also, from (2.2), w® = ug? + vp? — 2u/x. (2.25) 
3. Calculation of the retarding force.—F ollowing Bondi and Hoyle, we consider 
a star of mass M moving with velocity —U through a distribution of interstellar 
matter. This matter, which we shall call ‘“‘the cloud”’, is taken to have the 
15* 
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general properties postulated by these authors. In particular, the thermal 
motions of its particles are assumed to be negligible. We take an origin O fixed 
at the centre of the star and consequently regard the particles as being projected 
at a large distance from O with velocity U relative to O. Here we contemplate 
only the conditions after a steady state is assumed to have been established. 

We wish to calculate the force by which the cloud retards the motion of the 
star. ‘The force we are about to determine is actually equal and opposite to 
that required to hold the star in a Newtonian frame of reference when the cloud 
streams steadily past, having velocity U at a large distance ““up-stream”’. So 
long as the force so determined is sufficiently small it may be taken to be the 
force we wish to calculate. 

If OX is the axis in the i-direction, consider any plane II normal to OX and 
meeting it at a point N a large distance up-stream from O. Let s be the position 
vector relative to N of any point of II. 

Let p(s) be the density of the cloud at the point s. It may be any function 
ofs. But p may not be supposed to vary in the direction of OX, if the conditions 
are taken to be strictly steady. For such a variation would result in a variation 
with time of the amount of material flowing past the star. 

Case of no collisions.—If we entirely exclude the possibility of collisions 
between the particles of the cloud, even in the vicinity of the line OX, then every 
particle will describe a hyperbolic orbit about O of the sort studied in Section 2 
with x. =yM, where y is the constant of gravitation. We first evaluate the force 
in this case. As usual in such work, we neglect the effect of particles colliding 
directly with the surface of the star itself. 

Let dS be an element of II at position s. ‘Then the mass of the particles 
crossing dS per unit time is 

p(s)U dS. (3-1) 


The change of momentum per unit mass of these particles is U—r,. Hence 
the force F upon the star exerted by particles which cross II within a circle of 
centre N and radius s is 


F = | p(s\(U—F,,)UdS, 


where the integral is over the area of this circle and U—r,, is given by (2.10), 
(2.11). 

If F has components X, P along and perpendicular to OX, when we use the 
expressions (2.10), (2.11) in (3.2), we find 


r p(s) dS 
1+s*U/4 pe’ 

p— 2U* [ p(s)sdS 
pm J 1+ 8U4Yp? 


‘This gives the general result required. 

Case of linear variation of density.—For any applications of the results it is 
almost certainly adequate to use a linear approximation to the variation of the 
density. ‘This case will certainly serve to illustrate the nature of the results. 
But we may also notice that the force P which arises from the variation of p(s) is 
in fact unaffected by quadratic terms in this variation. 


X=2U?| (3.3) 


(3-4) 
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Let j be the direction in which the variation of p(s) occurs. This is the 
direction of grad p, which must be perpendicular to i in the case of a strictly 
steady state, as already remarked. 

We take then 


p(S) =(1 + bs cos $)pp, (3-5) 
where pp, 6 are constants and ¢ is the azimuthal angle in II measured from the 
j-direction. 

Substituting in (3.3), (3.4) we then find, since dS =s d¢ ds and 
s=scos¢j+ssin dk, 


2 2774 
X = 2p rf In (: + — ; (3.6) 


y — Pus we #Ut 
Y = mbpop E i In (: + 2 )] ’ (3-7) 
where P = Yj. 


We shall now write 
X= p/ uF, t= EXo, (3.8) 
and when s, x are related according to (2.15) we then have 


s* = 2x,%. (3-9) 

So (3.6), (3.7) may be written in terms of € instead of s as 
X = 27po(p?/U*) In (1 + 2€), (3.10) 
Y = wbp,(u3/U*)[2é — In (1 + 2€)}. (3-11) 

We shall be interested in the case of é large compared with unity. In that 
case we may note that, approximately, 

z = bx, a (é large). (3.12) 

Recalling the derivation of (2.15) we have that X, Y in (3.10), (3.11) give 
the force on the star exerted by particles of the cloud whose orbits in the gravitational 
field of the star meet the axis OX at all distances up to x = %9€. 

Before discussing these results we have to obtain the corresponding ones in 
the next case. 

[It should be noted that (3.10), (3.11) do not imply that the force tends to 
infinity when Uo. For in any application of (3.6), (3.7) the value of s is cer- 
tainly bounded above by the presence of other stars. ‘The value of X, Y must 
not exceed the value given by (3.6), (3.7) when s is replaced by its upper bound. 
But it is then seen that the latter values tend to zero when Uo. 

This remark is intended mainly to remove any misapprehension of a mathe- 
matical sort. In most applications it seems that U is likely to be sufficiently 
large for the upper bound of the distance €xy to be more significant than the 
upper bound of s. But if the formulae are to be applied to a range of values of U, 
including very small values, then it must be remembered that the approximations 
appropriate to large € will not be valid throughout the range. ] 

Case of collisions on the axis.—-Under the conditions of the preceding case 
the orbit of every particle intersects the axis OX. We shall now calculate the 
consequences of supposing that all the material simultaneously arriving along 
such orbits at any element (x,x-+dx) of the axis coalesces there before continuing 
its motion in the gravitational field of the star. ‘The purpose of considering this 
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case will be discussed below (Section 5), where the way in which it differs physically 
from the case considered by Bondi and Hoyle will be noted. 
We shall continue to use the density law (3.5). 
From (2.15) we have 
sds =x, dx, (3-13) 
and so from (3.1) the mass of the particles reaching dx per unit time is m dx, where 
m=poxoU | (1+bscosd) dp =2mpyxyU. 3-14) 
“0 


Using (2.17) the total momentum of these particles in the negative j-direction is 


ZpoX ott dx | (1 + bs cos $) cos 6 db/s = 2rbpyxopu dx. 
J0 


Consequently the velocity in this direction immediately after coalescence is 
vy where, from (3.14), (3.15), 
Up = by / U. (3.16) 


Also it follows at once from (2.16) that the velocity up in the i-direction is 


Up = U. (3-17) 
It should be noted that m, uo, vy are all independent of x. 

The motion after coalescence is then given by (2.18)—(2.25) with uo, vp having 
the values (3.16), (3.17). 

We notice in the first place that (2.25) gives 


w* = U? + by?/U? — 2y/x. (3-18) 
It follows that all the material reaching OX in 0 <x<axg, where from (3.18) 
a=2/(1+6*x,*), (3.19) 


must remain gravitationally bound to the star. We shall return to this feature 
in Section 5, noting in particular the way in which it differs from the corresponding 
result in the conditions envisaged by Bondi and Hoyle. Meantime we calculate 
the force on the star exerted by material reaching OX in axy <x < &x9. 
Evaluating the rate of change of momentum of this material, the required 
force has components 
rbd, §Xy 
X=m| (U—u,) dx, Y=m v, dx, (3.20) 
/ at, / ax, 
where u,, v, are given by substituting the present values of up, vo, w in (2.23), 
(2.24). 
The resulting expressions are highly complicated. However, for our purpose 
we need write down only the dominant terms when we take bx,<1 and €>1. 
These are easily found to give 


X ~2mp F tag Y ~2nb, we (3.21) 
9U2 ’ Po [ja ° J 


The result to be noticed is that these are precisely the asymptotic forms of 
(3.10), (3.11) (remembering that In(1 + 2) ~In2é =In€ +1n2 ~Iné for large é). 

{In this case, as in the case of no collisions, further considerations have to 
be taken into account when U is small. ‘The results will be given elsewhere 
in connection with certain applications of the work: they are not needed in the 
present paper. ] 
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Case of accretion column.—The treatment in the preceding section will not 
apply to the particular case where the density of the undisturbed cloud is uniform, 
i.e. where b=0. In that case the coalesced material will proceed to move along 
the axis OX and will encounter other incoming material as it does so. It is in 
fact the case considered in detail by Bondi and Hoyle. It should be remarked 
that the term “coalesce”’ is not intended to have any special physical connota- 
tions. It is convenient merely as a general description of the state of affairs 
treated by the mathematics. 

Bondi and Hoyle show in particular that there is a value x=ax, such that 
all the cloud-material reaching the axis in 0 <x<ax, proceeds to move towards 
O, so that it becomes “‘ accreted’”’ by the star, while all the material reaching the 
axis along x >axy proceeds to stream away from the star. ‘They point out that 
the parameter a cannot in this case be evaluated by steady-state considerations 
alone. But they give reasons for inferring that the value of a must be between 
I and 2. 

The authors proceed to evaluate the force upon the star exerted by the material 
reaching the axis between ax, and x, when >a and obtain, in our notation, 


X = 2mpo(u*/U*)In€g. (3-22) 


This is the correct asymptotic form resulting from their equations, though their 
method actually yields Iné+ A in place of In €, as a closer approximation, where 
A is a constant which the method does not serve to determine. 

Bondi and Hoyle give reasons for using a “ cut-off value”’ of € in applications 
of (3.22). Here it need only be recalled that the existence of such a value depends 
upon the presence of other stars in the vicinity. ‘Ihe same cut-off effect must be 
considered to operate in all the cases we have considered. 

4. Discussion.—For present purposes we shall call the force denoted by X 
the resistance and that denoted by P the deflecting force. 

The resistance (3.22) as calculated by Bondi and Hoyle appears to depend 
upon the existence of their outward-flowing stream along the accretion axis. 
This depends in the first place upon the density of the cloud ahead of the star 
being sufficiently uniform. For, if the density varies appreciably in a direction 
transverse to the motion of the star, i.e. if the coefficient b is appreciable, the 
stream will not form. In this case, even if collisions between cloud-particles 
along the axis have their maximum possible effectiveness, equation (3.16) shows 
explicitly that the material will drift to one side of the axis instead of accumulating 
along the axis. 

In the second place, even when the density is uniform, the existence of the 
stream depends upon the attainment of a steady state. ‘I’his in turn depends 
upon the particle-orbits passing through the accretion axis. But perturbations 
by other stars, though insufficient to change the general character of the orbits 
below the cut-off value of €, may nevertheless cause them to miss the axis by a 
small amount and so prevent the formation of the stream. Such effects will 
apply more to particles associated with the larger values of s concerned (which 
are those contributing the main part of the force) than to particles associated with 
the smaller values of s. ‘Thus the latter particles may still produce an inward- 
flowing stream on the axis near the star in the manner of the Bondi—Hoyle— 
Lyttleton theory, even where the outward-moving stream is not established, or 
only partially established. 
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The main object of the present paper is to point out that the value of the 
resistance does not depend upon the realization of the special conditions required to 
produce streaming along the accretion axis. ‘This statement applies only to, the 
dominant term in the force X. Since, however, the problem is in any event 
somewhat idealized, only this term need be considered. ‘The conclusion then 
rests upon the agreement between (3.10) and (3.22) in respect of this term. 

The conclusion is what would probably be expected on general physical 
grounds. But in view of the complicated motion of the material on the axis 
in the case of Bondi and Hoyle, we would not feel justified in presenting a physical 
argument without the backing of the mathematical calculations. 

The chief point to notice is that the resistance is, to the first order, the same 
whether or not collisions between the cloud-particles occur along the accretion 
axis. For in this regard the respective derivations of (3.10) and (3.22) may be 
considered to deal with the extreme cases. ‘This being so, any intermediate 
case is expected to lead also to the same value. Indeed, (3.21) may be regarded 
as dealing with one such intermediate case. 

Fundamentally the conclusion is that the force is the same whether the 
particles act cooperatively over part of their motion, in the special manner 
considered by Bondi and Hoyle, or independently as in the derivation of (3.10). 
This is important because the derivation of (3.10) does not depend physically 
upon the assumption of a steady state. ‘The result will continue to hold for a 
non-steady state if pp is taken to be a suitable mean value of the density ahead of 
the star. It will also continue to hold if the path of the star is not rectilinear, 
owing to the deflecting force either of the cloud or of another star. 

Turning now to the deflecting force, (3.12) may be interpreted as meaning 
that this may be significant for the motion of the star if the density of the cloud 
varies appreciably over distances comparable with interstellar distances. For 
the cut-off value of € is such as to make xpé in (3.12) a distance of this order. 

The conditions which would make P significant are ones that would preclude 
the formation of an axial stream. Here the extreme conditions are those leading 
to (3.11), (3.21), respectively, and these results agree in their dominant terms 
as already noted. 

The deflecting force may, in general, not be important itself. It is, however, 
not clear how large the transverse density-gradient must be in order to upset the 
mechanism of accretion considered by Bondi and Hoyle. ‘This must remain 
a subject for future study, but it may be that the effect of a transverse gradient 
on the mechanism of accretion is important even when it is too small to make 
the transverse force important. 

5. Accretion._-From the non-convergence of the expressions for the forces 
X, Y it follows that these forces are contributed principally by those cloud- 
particles which are most remote from the star concerned but whose orbits are 
yet determined primarily by the gravitational field of that star and not by the 
field of neighbouring stars. Accretion, on the other hand, if it does occur, is 
contributed by cloud-particles whose orbits lie much closer to the star. It 
follows that various perturbation effects may modify the motion of the particles 
mainly responsible for the forces, while the mechanism of accretion still operates 
in the manner envisaged by Bondi and Hoyle. That is why it was important to 
see that the value of the forces does not in fact depend upon the Bondi—Hoyle- 
Lyttleton mechanism. 
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There is, however, some interest in seeing how a variation in cloud-density 
may tend to affect the accretion process. Some consideration of the problem 
has recently been given also by Gething.* 

In the first place we have to recall that when an accretion column has been 
established in accordance with the Bondi-Hoyle—Lyttleton theory, the incoming 
cloud-particles collide with material already in the column and not directly 
with each other. Now the consequence of a transverse density-gradient in the 
cloud is to give the incoming material a resultant transverse momentum. In 
that case there cannot therefore be an accretion column lying symmetrically along 
the accretion axis. If the resultant transverse momentum is sufficiently large, 
there can in fact be no accumulation of material along the axis for the incoming 
particles to collide with; the only collisions to be considered would be collisions 
along the axis between the incoming particles themselves. If also the cloud- 
density were sufficiently small, such collisions would be rare and there would 
be no accretion. 

However, for any given cloud-density pg, there must be a range of values of 
the density-gradient as measured by 6 which leads to a situation intermediate 
between that of the Bondi-Hoyle—Lyttleton theory (6=0) and the extreme case 
of no accretion. Such an intermediate case must have the general character of 
that treated in (3.13)—(3.19). 

In this case all the material arriving at the axis OX in 0 <x<axy, where 
a is given by (3.19), is captured by the star. Also, since we are now dealing with 
the case where d is small, this equation gives a=2. ‘Thus the amount of material 
‘captured’ is approximately the same as that ‘“‘ accreted”’ in the Bondi-Hoyle— 
Lyttleton theory. In the present case, however, the captured material may not 
all fall on to the star but some may go to form a distribution revolving in the 
vicinity of the star. ‘The latter possibility can be shown, however, to require 
a somewhat large value of b. 

We notice that, whereas the parameter a is indeterminate on steady-state 
considerations in the case h=0, it becomes determinate in the case of 640 here 
considered. 

It is important to appreciate how this case differs from that considered by 
Bondi and Hoyle. In the present case the incoming particles are assumed to 
‘““soalesce’’ on the accretion axis, in the sense considered in Section 3. ‘This 
means merely that, owing to collisions in the vicinity of the axis, the particles 
arriving there at any instant acquire a common velocity. Owing to the assump- 
tion that 6 is not zero (though it may be very small) this velocity at any point of 
the axis has a non-zero component perpendicular to the axis. It is part of the 
assumption that on account of this lateral drift the material concerned is not 
involved in further collisions with subsequently arriving incoming particles. 
In the case considered by Bondi and Hoyle, on the other hand, the material 
proceeds to move along the axis where it does encounter further incoming 
particles. ‘Thus, although the present case is intermediate between that of 
Bondi and Hoyle and that of no collisions, it does not yield either of these cases 
as limits. In particular, it is not possible to infer that a must be determinate in 
Bondi and Hoyle’s case by attempting to treat this as the limit of our case. We 
do not question that in Bondi and Hoyle’s case a is not determined by the steady- 
state conditions but depends upon the way in which the state is established. 

*P. J. D. Gething, M.N., 111, 468, 1951. 
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In fact one of us (K.N.D.) hopes shortly to publish some new calculations of 
this effect. 

The more interesting point is that the angular momentum about O per unit 
mass of the captured material is, from (3.14), (3.16), 

I (%b P 
ri A 7 xdx=b ie (5.1) 
taking a=2 as before. If we take yu to correspond to the mass of the Sun and 
take U =5 km/sec as a typical value (being that also used by Bondi and Hoyle), 
this gives about 1-4 x 10°° bcm? sec"!. 

We may compare this with the angular momentum per unit mass of the Sun 
itself. According to values given by Chapman’, this is about 5 x 10!4 cm? sec !. 
‘Thus the value given by (5.1) is greater than or equal to this if b-? is less than 
about 100 parsecs. Now d-! is the distance over which the cloud-density changes 
by 100 per cent according to (3.5), and 100 parsecs would be a large value for the 
radius of a cloud.t 

It is then tempting to conclude that any significant accretion of mass is likely 
to be accompanied by a significant change in the star’s angular momentum. 
Nevertheless it has to be noted that, although the density-gradient required to 
produce this result is so small, its component transverse to the star’s path would 
have to be in the same direction over a great length of this path. Without much 
further investigation it is therefore impossible to say whether the effect has any 
general significance in regard to the phenomenon of stellar rotation. It may 
serve merely to endow separate portions of accreted material with varying amounts 
of angular momentum about the star in the manner suggested by Gething in 
relation to Lyttleton’s theory of comets. 


Royal Holloway College : 
1951 October 22. 


*S. Chapman, .W.N., 108, 236, 1948. (See p. 251.) 
t W. H. McCrea, The Observatory, 70, 100, 1950. 





THE ALIGNMENT OF GALACTIC DUST 
T. Gold 
(Received 1951 December 5) 


Summary 

It is shown that the dynamical interaction between galactic gas and those 
dust particles which are elongated will, in certain circumstances, lead to a 
partial alignment of those particles. ‘The direction which is preferred by 
the long dimension of the particles is that of the relative velocity between gas 
and dust. 

The possible relation between this effect and the observed polarization 
of starlight is discussed. 





1. Statement of the problem.—We aim to find the probability distribution of 
the orientation of elongated solid particles which results from their collisional 
interaction with a tenuous gas. ‘The conditions we shall assume are those 
which are presumed to exist within the Galaxy. ‘The essential simplifications 
are then as follows. 

(a) Impacts between the solid particles and the gas molecules are discrete 
events; the problem is one of dynamics of collisions rather than fluid dynamics. 


(6) Impacts of the solid particles with each other are too rare to be taken into 
consideration. 

(c) Photon interactions result only ina negligible transfer of angular momentum 
to the particles compared with that resulting from molecular collisions. (‘This 
may be untrue in the close vicinity of bright stars.) 

(d) Both gas and dust are sufficiently tenuous for interpenetration to occur. 
(For a dust particle as small as 5 x 10 °cm and density 3, the distance travelled 
in gas of density 10-* until the relative velocity is halved, is of the order of 10” cm.) 

(e) Thermal velocities of the gas molecules are small compared with the 
relative velocity of the bulk motions. 

Owing to the different penetrating power of gas and dust, we expect the two 
media to pursue substantially different courses within the Galaxy; processes of 
collision and interpenetration of dust clouds with gas clouds must then be 
common. We shall consider the effect of such interpenetration on the orientation 
of elongated particles. 

In the following discussion we are unable to use precise data for the shapes 
of dust particles ; it is also uncertain how small the thermal velocities are compared 
with the bulk velocities. We shall therefore discuss at each stage the ideal 
condition, in which we consider needle-shaped particles of negligible thickness 
and ignore thermal agitation, and we shall then mention the effects of departures 
from these ideal conditions. 

For the description of the orientation of a particle we shall use spherical polar 
coordinates (0,4), the direction 9=0 being that of the relative bulk velocity 
between gas and dust. 
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2. Dynamics of an elongated dust particle. (a) The direction of the angular 
momentum vector.—Collisions between dust and gas are assumed to be inelastic. * 
The angular momentum given to the dust particle at such a collision is then 
P’ =r = mv, where F is the position vector from the centre of mass of the particle 
to the point of impact, m the mass of the molecule, and v the relative velocity. 
(If w is the angular velocity of the particle, we ignore w x r in comparison with v.) 

In the idealized case v is entirely in the direction @=0; all contributions 
like P’ are thus confined to the plane @ = 7/2, and thus the total angular momentum 
P remains always in that plane. 

In the actual case the total angular momentum after n collisions is 


P=P,+<P,’ or P,+ =r, m,v,, 
n n 
where P, is the initial angular momentum. ‘The departure of P from the plane 
( =7/2 now arises from P,, which we take to be of thermal origin, and from the 
thermal variability of the direction of the v,. In the case when the original 
temperature of the dust equals that of the gas, the variability of the v, does not 
add to the spread of direction of P, which is already fully accounted for by the 
term P,. P is thus brought close to the plane 6=7/2 when <r, x my is large 


n 


compared with a term like Py, which may be taken to represent the initial thermal 
angular momentum of the particle, or that appropriate to the temperature of 
the gas, whichever is the greater. This condition can be satisfied only if the 
velocity of interpenetration is supersonic. 

(6) The relation between the angular momentum vector and the orientation of 
the particle.—We denote the principal moments of inertia of the particles by 
I, 1,, 1,, the x direction being taken to coincide with the long dimension. In the 
idealized case the vector r must lie entirely in the x direction, since the particle is 
assumed to have negligible y and z dimensions. All contributions to the angular 
momentum P,' are in the y-z plane; and in the simple case when P is confined 
to the plane @ = 7/2, the x direction of a particle remains in a plane of constant ¢ ; 
each particle has, therefore, to move through the polar direction 6=0, whilst 
there is no such preferred direction in the equatorial plane. 

In areal case where the particle has certain dimensions in the y and z directions, 
two new considerations arise. Firstly r is not entirely confined to the x direction, 
and the angular momentum imparted at an impact is thus not entirely confined to 
the y—z plane. Secondly J, is different from zero, and as the particle performs 
free nutations between collisions, its x-axis is then not completely constrained 
to move ina plane. We now aim to show that these two effects do not completely 
destroy the previous argument, but that 4 restriction on the y and z dimensions 
of the particle in comparison with the x dimension will imply a restriction of the 
angle by which the motion may carry the x direction out of a plane of constant ¢. 


* This will in fact be a close approximation for the conditions which are envisaged. ‘The 
temperature 7 of the dust particle is low, being of the order of the radiation temperature of the 
space. Gas molecules with kinetic energies much in excess of $kT will lose most of their energy 
to crystal lattice vibrations, and at higher energies also to lattice displacements and ionizations. 
This energy is subsequently lost by the dust particle through radiation, with a negligible momentum 
transfer. ‘The gas molecule may then be embedded or stick, or move off with a small velocity; 
the momentum transferred to the dust is very nearly that of the incoming molecule. The value of 
the minimum velocity for collisions to be inelastic depends on the substance of the dust, but is 
unlikely to be more than one kilometre per second. 
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Most position vectors r, (and all the large ones) will be still dominantly in 
the x direction, and most contributions P,’ (which are all in the fixed plane 
§ =7/2) will therefore be close to the y—z plane of the particle; they will thus be 
associated with the two major axes of inertia J, and J, rather than with /,. 

To consider the free motion of a particle it is convenient to use Poinsot’s 
theorem. ‘The condition that the angular momenta P,’ are principally in 
the y—z plane then implies that the “ ellipsoid of inertia”’ is rolling with its centre 
close to the “invariable plane’. As the distance of the centre from the invariable 
plane cannot change for a free motion, the long axis of the ellipsoid (if it is much 
longer than the other two, i.e. if 7, is much smaller than J, and J,) cannot depart 
far from parallelism to the invariable plane. If, as before, the invariable plane 
is constrained to be a plane of constant ¢, the x (long) axis of the ellipsoid, and 
hence of the particle, is constrained to remain nearly parallel to it. 

3. Probability distribution of orientation.—We define an ‘‘ alignment factor’ 
for an assembly of particles by projecting each particle into one particular plane 
of constant ¢ (say ¢ = 0); in that plane we then sum all components in the direction 
?=o and in the direction 0=7/2, and we then take the ratio of the two sums. 
(This definition may be of value in considering the optical polarization effect 
on light travelling in or near the @=7/2 plane.) 

In the idealized case all planes of constant ¢ are equally populated by spinning 
particles, and within each such plane all angles of @ are equally probable. The 
alignment factor is then simply given (for arbitrarily large numbers of particles) by 


on[2 -n/[2 n/2 on/2 
| dé| cosé ao | | dd| sin 0 cos ¢ dé = 7, 2. 
0 ~ 0 “0 “0 


We have not calculated the alignment factor for any other case; but we 
estimate, for example, that it will differ from 7/2 by only a few per cent if the 
particles are more than ten times as long as they are thick. 

4. The polarization of starlight.—It now seems probable that the polarization 
of starlight has to be explained by selective absorption in interstellar space. 
Hence, one enquires into agencies which will result in alignment of elongated 
dust particles, as no other suitable absorbing medium is known. 

Previous explanations (1, 2) have used the hypothesis that sufficiently strong 
interstellar magnetic fields exist, and that a sufficient fraction of elongated particles 
have the required ferromagnetic or hysteresis properties. ‘There has been no 
independent evidence for such strong magnetic fields, and certain theoretical 
considerations show the required strength to be improbable; and there has been 
no evidence for the otherwise implausible hypothesis that a major fraction of 
the dust possesses those particular physical properties. 

The present effect must be 2xpected to operate on all those particles which 
are elongated enough to have the desired optical property. ‘The aligning process 
of double streaming is itself a necessary consequence of conditions which are 
generally believed to exist in the Galaxy. 

Although there will be a diminution of the relative velocity in the course of 
an encounter between a gas and a dust cloud, such relative velocities cannot 
generally die out so long as there are encounters between clouds at all. For if 
we consider an encounter between two mixed clouds in which the dust possesses 
the same bulk motion as the surrounding gas, a separation would again occur, 
the motion of the gases being much more severely modified by the encounter 
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than that of the dust. A similar conclusion is reached if one considers any kind 
of turbulent motion of a mixture of gas and dust. 

The finite thickness of the Galaxy implies that all the material above and 
below the central plane must at some stage penetrate through it; on dynamical 
grounds we must expect the maximum interpenetration of material to occur in 
the central plane, with the relative velocity directed normal to it. Double stream- 
ing within the central plane in other directions may also exist, but will be subject 
to much more rapid damping. 

The observed polarization of starlight shows in most cases that the electric 
vector normal to the galactic plane is attenuated more than the vector in the plane 
(3). This sense of the effect, and its coherence over wide regions of the Galaxy, 
agree with the effect which must be expected from the process of double streaming. 
The maximum magnitude of the observed effect amounts to a variation of only 
5 per cent in the extinction coefficient of the dust, depending on the orientation 
of the electric vector (2). 


Trinity College, 
Cambridge : 
1951 December 3. 
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Summary 


This paper describes radio echo observations of the major night-time 
streams for the period 1946-1951. Radiant positions and velocities have 
been measured for the Perseids, Geminids and Quadrantids, and orbits 
have been computed for each stream. ‘The characteristics of these streams 
are discussed, and evidence is presented for the existence of a core of larger 
meteors in the region of maximum density in the streams. A progressive 
movement of the radiant point has been detected for the Perseid and Geminid 
showers, which is in agreement with previous observations. ‘The velocity 
distributions indicate that there may be a spread in heliocentric velocity 
within the streams of the order of 1-2 km/s. Radiant positions and hourly 
rates are also given for the Orionids, 8 Aquarids, Lyrids and the night-time 
‘Taurids and Arietids. 





1. Introduction 


In the Northern Hemisphere the richest night-time meteor displays recurring 
annually are the Perseids, Geminids and Quadrantids, which, together with less 
active displays such as the Orionids and Lyrids, are classed as the major night-time 
meteor streams. Quetelet (1x) first discovered the recurrent nature of the Perseid 
shower in 1836, and the regular displays of the Geminids (6) and Quadrantids 
were commented on by Herrick in 1838-39. ‘These streams have since been 
exhaustively studied by visual methods (1-13) as far as moonlight, clouds and 
other adverse factors will allow. Early research was impeded by the inaccurate 
determinations of meteor velocities, but after Elkin (14) had initiated the technique 
of interrupted photography reliable results became possible, and Whipple and 
others have recently determined exact orbits for several bright Geminid (15) 
and Perseid (16) meteors. ‘The development of radio echo methods has enabled 
measurements of velocity and radiant position of the.fainter and more numerous 
meteors of a stream to be made independently of visual observing conditions. 
Studies by these radio techniques of the Geminid, Perseid, Quadrantid and other 
less active streams are described in this paper. 


2. Methods of observation 


(i) Radiant co-ordinates.—The apparatus used for radiant determination has 
been fully described elsewhere (17). By recording the range and time of 
appearance of radio echoes from meteor trails in two azimuthal directions it is 


* Received in original form 1951 November 3. 
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possible to measure the radiant co-ordinates of a meteor stream with an accuracy 
ofabout +1°. ‘This accuracy represents the ideal limit obtained with a favourably 
placed point radiant supplying about 50 echoes an hour. For radiants of weaker 
activity the errors are greater, but seldom exceed 5°. It is also possible to detect 
a spread of meteor directions about the mean from the scatter of points on the 
range-time plots. ‘This scatter of points gives a measure of the mean diameter of 
a radiant, although the precise shape of a radiating area cannot be found. In cases 
where a finite radiant diameter exists, the accuracy of position determination 
is again reduced and it is difficult to deduce the exact centre of maximum activity 
within the radiating area. 

(ii) Meteor velocities—'l'wo equipments operating on wave-lengths of 4m 
and 8 m were used for measuring velocities, as in previous experiments (18, 19), 
and it has been shown that these equipments are capable of an experimental 
accuracy of +5 percent. ‘lhe use of a narrow beam aerial reduces the number 
of sporadic meteors included in the records. 

(iii) Hourly rate of echoes—The most reliable measure of the activity of a 
stream was given by the radiant survey apparatus (17), where the transmitted 
power, receiver gain and other parameters were checked regularly in an attempt 
to obtain constant operating conditions. A further check on the stability of the 
operating conditions was provided by the observations of the rate of sporadic 
meteors. On the assumption that the mean rate of sporadic meteors remains 
constant at any particular position in the Earth’s orbit, any variations in the 
observed radio echo rate can be used to calibrate the equipment from year to year. 
‘The sensitivity of the apparatus was set at a value which gave an hourly rate of 
echoes comparable to that obtained by a single observer under good visual 
conditions. 

Additional hourly rate measurements have been obtained from the velocity 
equipments, and also from the records made with the apparatus used in the early 
experiments (20, 21). 

(iv) Long-duration echoes.—In early experiments, durations were estimated 
by visual observation of a range—amplitude cathode ray tube display. In the 
more recent work, durations have been recorded by the continuously moving film 
of the radiant survey apparatus. As an example of the method, Fig. 1 (Plate 3) 
shows a record of echoes displayed on a range—intensity cathode ray tube, photo- 
graphed on a film moving at 72cm/h across the range axis. A short-duration 
echo appears as a bright dot*, whereas a long-duration echo is extended along 
the time axis so that its duration can be measured to within +18, Additional 
information on long-duration echoes has been obtained from the photographic 
recorders (18) used in conjunction with other equipments. ‘These recorders 
were triggered by a radio echo and photographed the variations in echo amplitude 
for a period of 0*-1 at intervals of 18. By this means the lifetime of long-duration 
echoes may be measured to the nearest second. 


3. Observational data 

(i) Radiant positions.—Observations have been made of the radiant co-ordinates 
of the night-time streams over the period 1949 September to 1951 August, and 
the results for the Perseids, Geminids and Quadrantids are given in ‘Table I. 


*A single echo appears double on the display because twin transmitter pulses were used, 
separated by 50 km in range. ‘This transmitter code facilitated the identification of echoes. 
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These measurements were made at the time of transit of the radiant and the 
change of date was made at 0? U.'T. From the characteristics of the range-time 
plots (17) the radiants have been divided into three classes, A, B and C, weighted 
in order of importance as 2:1:0. A weighted mean radiant position may then 
be deduced as in Table II. For the less active streams the observational errors 
are somewhat increased, so that day-to-day positions have not been given. In 
these cases the weighted mean position (W) has been given in ‘Table I and the error 
in this value is about + 3°. 

(ii) Meteor velocities.—Velocity observations have been made on the night-time 
streams during the period 1947 December to 1951 August. No significant 
variations were found in the velocities of each stream from day to day and the 
mean velocity for each year is given in Table I. Similarly the distributions of 
velocities obtained in each year did not justify separate treatment, and all velocities 
obtained for the Perseids, Geminids and Quadrantids have been combined into 
the histograms of Fig. 2. In the case of the Orionids, Leonids and other less 
active streams it was not possible to obtain enough velocity measurements to 
form a reliable distribution. 


TABLE I 
Radiant 
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Quadrantids 
1947 Jan. 
1948 Jan. 
1949 Jan. 
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Lyrids 

1947 Apr. 
1948 Apr. 
1949 Apr. 
1950 Apr. 
1951 Apr. 
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* Apparatus unstable and figures less reliable, so that the period 1950 December 28 to 1951 
January 7 was taken as the reference epoch for sensitivity corrections. 

t In 1951 a subsidiary burst of meteor activity was noticed near the Quadrantid radiant and 
has been included in the Quadrantid sequence. However, no velocity measurements were made, 


so that it is impossible to ascertain the exact relationship of this later radiant with the main 
Quadrantid stream. 
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Radiant 





‘ Echo H.R. 
No. of —_ -~ Mean 
R.A. Dec. Class Diam. echoes © Obs. at |Red.to vel.in No. of 
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& Aquarids 

1949 July 29 ... ve an igi a 
1950 July 339 115 
1951 July 27 336 87 


Perseids 
1946 Aug. 
1947 Aug. 
1948 Aug. 
1949 Aug. 
1950 July 
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Night-time 
Arietids 
1950 Oct. 14 37 ! W 130 6199°5 19 





* On 1950 August 12 and August 13, 50 per cent of the echoes corresponded to a sharp sub-centre, 
the remainder were scattered in a circle of diameter 10°. 
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TABLE I (contd.) 
Radiant 
A —~ Echo H.R. Mean 
No. of r A . vel.in No. of 
Date .A. Dec. Class Diam. echoes ©) Obs. at| Red. to atmos. vels. 
used transit} 1950 km/s  measd. 








Orionids 

1946 Oct. 

1947 Oct. 

1948 Oct. = ‘om pos 
1949 Oct. +13 W 
1950 Oct. +g W 


Leonids 

1946 Nov. 
1947 Nov. 
1948 Nov. 
1949 Nov. 
1950 Nov. 


Night-time 
Taurids 

1946 Nov. 
1947 Nov. 
1950 Nov. 


Geminids 
1946 Dec. = ie uae =e we is “ee = 
1947 Dec. ne sed = ree be oe 34°4+1'2 10 
1948 Dec. : i iis ove 
1949 Dec. 
Dec. 
Dec. 
Dec. 
Dec. 
Dec. 
Dec. 
Dec. 
Dec. 
Dec. 
Dec. 
Dec. 
Dec. 
Dec. 
Dec. 
Dec. 
Dec. 


A 


A 


I Ph WOREAWOWYUAWHWWWWYW! 


A 


\ 
/\ 


QPF rrr rDNA SPrewaeTre: 





Ursids 

1947 Dec. 
1948 Dec. 
1949 Dec. 
1950 Dec. 
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(iii) Hourly rate of echoes.—The hourly rates of the various streams are 
given in Table I. For the three richest streams day-by-day rates are given as 
recorded by the radiant survey apparatus. ‘This rate is a mean taken over the 
range-time envelopes of each aerial. ‘The mean sporadic rates over the period 
of each shower were as follows :— 

Quadrantids Perseids Geminids 

1949 stl sa 2°8 

1950 3°4 g°2 a7 

1951 5:0 II‘l 55 
Approximately 400 hours of observation have been used in evaluating each figure, 
and the statistical variation expected from purely random fluctuations is +0-I 
meteors an hour. ‘These figures have been assumed to give a measure of the 
relative sensitivity of the apparatus, and the observed hourly rates in each case 
have been reduced to the sensitivity of the corresponding period in 1950. 
Corrected rates for three streams are plotted against the Sun’s longitude in 
Figs. 3(a), 4(a) and 5(a). ‘The shape of these curves is confirmed by the mean 
hourly rates (Figs. 3(6), 4(6), 5(6)) deduced from the equipments referred to 
in Section 2(iii), between 1946-1951, but the standard deviation introduced 
by these inhomogeneous data is necessarily large. 

The remaining streams, the Lyrids, 5 Aquarids, night-time Arietids, Orionids 
and night-time Taurids are on the limits of observation by present radio echo 
techniques. However, the measured radiant positions, maximum hourly rate 
and date of maximum activity are given in Table I. The figures of maximum 
rate necessarily include a number of sporadic meteors and hence these figures 
must be taken as an upper limit to the activity of these streams. 


4. Determination of orbital elements 


The orbital elements given in Table II have been computed from the radio 
observations of the Perseids, Geminids and Quadrantids. ‘The limits of error 
in the orbital elements have been derived from the estimated error of radiant 
determination, and the probable error of the velocity distributions. 

Corrections for zenithal attraction and diurnal aberration have been applied 
to the mean radiant positions, and the mean velocities have been corrected for 
the attraction of the Earth and for atmospheric deceleration. An approximate 
value for the deceleration has been obtained from the equations proposed by 
Whipple and Opik (22), which govern the dissolution of rapidly moving bodies 
in the upper atmosphere. For a meteor of mass m, moving with a velocity 7, 
at a height / and incident at a zenith angle x, the laws of conservation of momentum 
and energy give 

m3 dy = yApyve~" sec x dh, (1) 
Ldm = }AyAm**pyv*e-"4 sec x dh, (2) 


where H is the scale height and py the air density at h=o. ‘The constant 4 
depends on the shape of the meteor particle, y is the fraction of air molecules 
trapped in its surface, A is the efficiency of the heating process and / is the energy 
required to evaporate unit mass of the particle. 

Division and integration of these two equations gives 


m =m. EXP [( ~ 5) (0.7 — “) | , 
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where m,, and v,, are the mass and velocity of the particle outside the Earth’s 
atmosphere, 


Substituting in (1) for m and integrating again yields 


2yAHm;,"*p, sec x =e"* {> x 'e-° dx, (4) 
z 
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Fic. 2.—Histograms of Perseid, Geminid and Quadrantid velocities. 


Ordinates: Number. 
Abscissae t Observed velocity in atmosphere (km/s). 


where x =Av*/12/ and p, is the air density at the point where the deceleration is 
required. 
Since the deceleration is small we may write 
ayAHmz 8p, sec x = }(x'z'e-*@ +37 4e-*)(%_—2), (5) 




















a 
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from which it is possible to deduce (v.,—v), the change in velocity experienced 
by a meteor during its passage through the atmosphere. 

The mass of the meteor outside the Earth’s atmosphere, m,,, has been derived 
from temperature measurements of the trails by assuming that the major part 
of the kinetic energy is converted into heat (23). For the meteors studied in 


rS ‘ 
° ---- Denning 


Hawkins x 2 














Fig. 3c 
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Fic. 3.—(a) Hourly rate of Perseid stream, corrected for variations in sensitivity, plotted 
against Sun’s longitude. (b) Mean hourly rate observed by all equipments plotted against Sun's 
longitude. The vertical lines indicate the probable error in the mean values. (c) Percentage of 
echoes with durations greater than one second plotted against the Sun’s longitude. The vertical 
lines indicate the probable errors in the values. 


this work the masses were found to be of the order of 0-005g. The constants 
A(=0'5), 1(=6 x 10! erg/g) and A(=0-5) have been adopted from the data of 
Whipple (22) and the upper-air data have been taken from Jacchia (25). ‘The 
mean height, H, derived from radio echo data (23, 24) was taken as 95 km for 
Perseid meteors and 85 km for Geminid and Quadrantid meteors. 
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Substitution of these data in (5) gives the following corrections to be applied 
to the observed velocities : 

Perseids 0-3 km/s, Geminids 0-4km/s, Quadrantids 0-4 km/s. 

The orbital elements corrected for this deceleration are given in Table II. 


5. Discussion 


(i) The Perseid stream.—The mean radiant position of this stream, and the 
large diameter of the radiating area, are in agreement with visual observations 
(34,5). The curve of corrected echo rate against Sun’s longitude in Fig. 3 (a) 
gives the density of matter across the stream in the cross-section traversed by 
the Earth. Early visual observations on stream density obtained by Denning 
(9), together with more recent observations, are plotted for comparison in 
Fig. 3(a). It can be seen that the position of maximum density of the stream 
has remained fixed over the period of observation. There is evidence of a 
change in mass distribution across the stream as shown in Fig. 3(c), where the 
percentage of echoes of duration greater than I* is plotted against Sun’s longitude. 
Recent work (26) on the scattering of radio waves from meteor trails has shown 
that the duration of the radio echo is a measure of the ionization density and 
hence of the mass of the meteor. The results shown in Fig. 3(c) therefore 
indicate that the greatest proportion of large meteors is concentrated in regions 
of greatest stream density. A further characteristic of the core is disclosed by 
the contraction of the radiant area towards maximum, a feature of both radio 
and visual observations (4). This contraction indicates that meteors in the centre 
of the stream are less disturbed than meteors on the fringes. 

The association of the Perseids with Comet 1862 III was demonstrated by 
Schiaparelli (2) in 1866. From Table II it can be seen that the radio results 
confirm this association. It has been suggested (27, 28) that the Perseid meteors 
could be produced by an ejection of material from the cometary nucleus at 
perihelion with subsequent secular perturbations by the planets, which, in time, 
would lead to a diffuse distribution of meteor orbits, with a spread in heliocentric 
velocity. The radiant area confirms the spread in orbits, and the evidence for 
a spread in heliocentric velocity is considered later (Section 5 (v)). 

(ii) The Geminid stream.—Radio echo, visual (7, 8, 31) and photographic (15) 
determinations of the radiant position and of the velocity of the stream are in 
agreement, but there is disagreement in the measurements of radiant diameter. 
The radio observations show that most meteors appear from a radiant area of 
diameter not greater than 4°, whereas the visual paths define a radiant diameter 
of about 12°. There is, however, a visual concentration conforming to the radio 
position, 50 per cent of the meteors delineating an area 4° x 3° around the mean 
visual position (29). A comparison of the measurements of stream density is 
shown in Fig. 4(a@) and it can be seen that both techniques show an asymmetrical 
distribution of material about the mean orbit, the Earth encountering less débris 
after maximum than before. The percentage of long-duration echoes shown in 
Fig. 4(c) also follows this asymmetry. Thus the Geminids possess a core of 
heavy meteors displaced to one side of the stream. In contrast to the Perseids, 
however, the radiant areas do not show any marked tendency for improved align- 
ment of meteor paths as the core is traversed. 

(iti) The Quadrantid stream.—Visual determinations (4, 10-13) of the radiant 
of the Quadrantids show a scatter of +7° about a mean value of 232°, +52°. 
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‘I'he present radio determinations fall within these limits. The radiant area is 
found by both techniques to be large, and visually (29) there is evidence of perhaps 
three sub-centres of radiation. In 1951 the radio data suggested the possibility 
of two centres of equal intensity at 239°, +49° and 234°, +54°, but the radiants 
were on the limits of resolution and a mean position has been taken. The 
contraction of the radiant area at maximum is again pronounced, but the evidence 
for a core of large meteors is not conclusive. The similarity of the radio and 
visual data (10) regarding the duration of the stream is shown in Fig. 5§(a). Both 
results indicate that the stream is of short duration. 
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Fic. 4.—(a) Hourly rate of Geminid stream, corrected for variations in sensitivity, plotted 
against Sun’s longitude. (b) Mean hourly rate observed by all equipments plotted against Sun’s 
longitude. (c) Percentage of echoes with duration greater than one second plotted against Sun's 


longitude. 
The orbit of this meteor stream is unique in that it exhibits both a short 
period and high inclination to the plane of the ecliptic, and is therefore markedly 
dissimilar from the orbit of any comet, minor planet or meteor stream at present 
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known. This exceptional orbit is confirmed by the radiant motion (13), and by 
photographic observations by Jacchia (30) of one Quadrantid meteor in 1951, 
the orbital elements of which are given in Table I]. The mean radio echo orbit 
is drawn in Fig. 6. 
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Fic. 5.—(a) Hourly rate of Quadrantid stream, corrected for variations in sensitivity, plotted 
against Sun’s longitude. (b) Mean hourly rate cbserved by all equipments plotted against Sun’s 
longitude. (c) Percentage of echoes with duration greater than one second plotted against 
Sun’s longitude. 


(iv) Hourly rates.—The hourly rates of the Geminid and Quadrantid streams 
when reduced to the sensitivity of 1950 are constant within the limits to be 
expected from statistical fluctuations. For these short-period streams this 
implies that the material is uniformly distributed around the length of the orbit. 
For the Perseids, the observations do not cover a long enough period for the 
presence or absence of a periodicity to be established. ‘The absence of a marked 
peak in 1951 is conspicuous and outside the limits of statistical fluctuations. This 
may indicate local irregularities in the density of the Perseid stream. 
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Although the radio sensitivity was set to produce a rate comparable with the 
visual rate, there is some doubt as to whether the ratio of the rates remains constant 
throughout the various showers. So far it appears that the visual activity (29) 
of the Perseids, Geminids and Quadrantids are in the ratio 50:60:40, whereas 
the corresponding radio echo ratios are 60:80:90. Further observations are 
necessary to establish if any significance can be attached to these differences. 


& 





Fic. 6.—Mean orbit of the Quadrantid stream drawn in the plane of the stream. 


(v) Distribution of stream velocities.—In dealing with the origin and subsequent 
perturbations of a stream, any spread in the heliocentric velocity of each stream 
is significant. ‘The observed spread in geocentric velocity is caused in part by 
experimental errors in measurement, and in part by external factors, such as the 
inclusion of sporadic meteors in the sample. ‘Thus the standard deviation op 
of the observed geocentric distribution is given by 


Ge" =0,"+ CR” +a," +o," +o,’, 


where o, =standard deviation caused by random errors in the experimental 

technique, 

o, = standard deviation which would result from the variation of apparent 
elongation of radiant points, 

og=standard deviation introduced by sporadic meteors, 

o,=standard deviation caused by differential atmospheric deceleration, 
and 

oy = Standard deviation introduced by a real spread of heliocentric velocities 
in the stream. 


Now og may be found from the errors of individual velocities (19), on from the 
observed radiant positions and diameters, og from the measured background 
rate of sporadic meteors and a, may be estimated by using equation (5), assuming 
that two-thirds of the meteors observed were contained in a layer of the atmosphere 
of depth 15km. It is therefore possible to estimate og and hence the spread in 
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heliocentric velocity: this has been done in Table III. This indicates that 
there may be a velocity spread of I or 2km/s in each stream, though the errors 
in oy are of the same order as og and reduce the significance of the result. The 
spread in velocity found from the radio echo observations is considerably greater 
than the value found from photographic observations. In the case of the 
Geminids, Whipple (15) has published seven measurements of the velocity of 
individual meteors, and these velocities have a true spread of 0:23km/s. This 
discrepancy may be due to the different range of mass groups observed with the 
two methods. 


Tasce III 


Perseids Geminids Quadrantids 
Og 1°7 km/s 1:2 km/s 1-0 km/s 
Og 2°9 3°5 I‘! 
Op 2°! 2° 1°7 
4, 073 073 073 
V oy? + 05"+ 0+ 0,4? 4°0 4°0 
Og 4°6 4°6 





Oy 2°3 2°3 


(vi) Motion of the radiant point.—Applying the Gaussian method of weighted 
least squares it is possible to derive the mean daily motion of the radiant position, 
and this is given in Table IV for two of the radiants. ‘The Geminid and Perseid 
radiants are found to have a diurnal change in R.A. of 1°-3 and 0°-6 respectively, 


which confirms the eastward drift found by visual and photographic methods 
(3, 7, 8,16). No significant variations were found for the change in declination 
of these radiants. 

TABLE IV 


The Daily Motion of the Perseid and Geminid Radiants 





Perseids Geminids 


R.A. 





Jodrell Bank : Jodrell Bank 
(Radio Echo) to:rto: (Radio Echo) 
Whipple (16) Whipple (16) 
(Photographic) Lo +or (Photographic) + 1°1+0°1 
Denning (3) King (8) 

(Visual) + 1° “4 (Visual) ‘2 
Maltzev (7) 

(Visual)* +1'0 


° > 


r°4 7 0°% 








* The value given by Maltzev is based on a reinterpretation of King’s original data. 


A daily eastward motion of 1° parallel to the ecliptic maintains the radiant at 
a constant elongation from the apex of the Earth’s way, and these results indicate 
that the meteors in the Perseid and Geminid stream move in orbits which show 
a progressive rotation in the line of apsides. In the case of the Quadrantids, 


more results are required before it will be possible to give a reliable value for the 
radiant motion. 
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ON THE INTERNAL STRUCTURES OF JUPITER AND SATURN 


B. Miles and W. H. Ramsey 
(Received 1952 January 24) 


Summary 


One of the authors (/.N., 111, 427, 1951) has shown that the proportion 
of hydrogen in Jupiter is 76 per cent by mass and in Saturn 62 per cent by 
mass, assuming that the heavier elements are distributed uniformly throughout 
each planet. In the present calculations the heavier elements are assumed 
to be in part concentrated in a central core and in part distributed uniformly. 
This relaxation of the condition of chemical homogeneity has little influence 
on the estimated hydrogen contents, the extreme estimates for Jupiter being 
76 and 84 per cent by mass and for Saturn 62 and 69 per cent by mass. The 
ellipticities of the planets due to their axial rotations suggest that each planet 
has a central concentration of heavy elements of the order of ten times the 
mass of the Earth; but the possibility of chemical homogeneity is not completely 
excluded by the present empirical data. 





1. One of the authors has shown in a recent paper (1) that hydrogen is the 
predominant constituent of both Jupiter and Saturn. The ordinary molecular 
modification of solid hydrogen is stable only at the low pressures in the outer 
layers of the planets; in the deep interiors hydrogen will be metallic. Metallic 
hydrogen is the simplest of all solids from the point of view of atomic physics, 
and its pressure—density relationship can be computed accurately. ‘The computed 
densities of solid molecular hydrogen, though less reliable than those of metallic 
hydrogen, are sufficiently accurate for calculations on the bulk characteristics of 
the planets. The theoretical pressure—density relationship refers to the absolute 
zero of temperature, but the temperatures in the interiors of the planets are too 
low to alter the densities appreciably ; likewise, the presence of an atmosphere is of 
little importance, since its contributions to the mass and volume of the planet 
are negligible. Both Jupiter and Saturn are denser than hydrogen planets of 
the same mass, but by less than a factor of two. The relatively small proportions 
of heavier elements necessary to account for the observed mean densities were 
assumed in these calculations (1) to be distributed uniformly throughout each 
planet; in other words, the planets were assumed to be chemically homogeneous. 
The proportion of hydrogen in Jupiter was estimated to be 70 per cent by mass, 
and in Saturn 62 per cent by mass. ‘The moments of inertia computed from the 
theoretical density distributions were close to, but definitely larger than, the 
values derived by means of the Radau—Darwin formula from the polar flattenings 
of the planets due to their axial rotations. 

The principal weakness in the earlier calculations (1) is the assumption that 
the planets are chemically homogeneous. Each planet could contain a central 
concentration of heavy elements, as has been suggested by Jeffreys (2) and by 
Wildt (3). A central core of heavy elements is to be expected if the planets were 
formed by an aggregation process; only the heavier and less volatile materials 
would have been captured when the mass of the planet was small. Also, even if 
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a planet were initially chemically homogeneous, the heavier elements may have 
gravitated towards the centre. The discrepancies between the empirical and 
the computed moments of inertia indeed suggest central concentrations of heavy 
materials; but it will be shown later in the paper that the disagreement is partly 
due to the inadequacy of the Radau-Darwintreatment. Inthe present calculations 
the heavier elements will be assumed to be in part concentrated in a central core 
and in part distributed uniformly. It will be shown that this relaxation of the 
condition of chemical homogeneity has little influence on the estimated hydrogen 
contents of the planets. The central concentrations of mass influence primarily 
the computed moments of inertia of the planets and the ellipticities arising from 
their axial rotations. The extents of the central condensations will be estimated 
from the empirical data by a purely numerical method which is more accurate 
than the Radau—Darwin approximation. 

2. If the density p of the material of a planet is known as a function of the 
pressure p, the bulk properties of the planet can be computed on the assumption 
of hydrostatic equilibrium. In the absence of rotation, the equation to be solved 


1S 


dp/dr = — pg(r), (1) 


where r is the distance from the centre of the planet and where g(r) is the local 
acceleration due to gravity: 


g(r) =(47Gir9) pl a)x* de, (2) 


G being the constant of gravitation. Equation (1) can be used for the Earth, 
since the centrifugal force is never more than 0-3 per cent of gravity. But on 
the surface of Jupiter the centrifugal force at the equator is 9 per cent of gravity, 
and on Saturn 16 per cent of gravity; it must therefore be taken into account 
in the hydrostatic equation (1). This can be done by the method adopted in 
the earlier paper (1). Any pressure p defines in the planet an oblate surface 
which is approximately a spheroid. Let r be the radius of the sphere which 
contains exactly the same volume as this surface. The pressure p is a function 
only of the parameter 7, which will therefore be adopted as the independent 
variable. It can easily be shown that the equation for hydrostatic equilibrium is 


dp|dr = —p{g(r) — ur}, (3) 
w being the angular velocity of axial rotation; as before, g(r) is defined by 
equation (2). This equation is to be solved as though the planet had spherical 
symmetry; the oblateness is allowed for in the definition of the parameter r. 
Equation (3) is not exact, fourth and higher powers of the angular velocity w 
having been ignored ; but it has been shown in reference (1) that it is adequate for 
the present purposes. This equation will be used throughout the paper; the 
angular velocity assumed for Jupiter corresponds to a rotational period of 9" 55™, 
and for Saturn to a period of 10% 38™, 

Equation (3) will be integrated numerically on the assumption that the planet, 
apart from a central concentration of heavy elements, consists of a chemically 
homogeneous material whose pressure—density relationship is known. It is 
most convenient, as Bullen (4) found in analogous calculations on the terrestrial 
planets, to start the integration at the surface of the planet and to integrate inward 
to the core of heavy material. A single numerical integration is valid irrespective 
of the dimension of the central core. ‘The mass of the central core depends on its 
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mean density, but it will be shown that the dependence is so slight that the density 
may be assumed to be infinite in these calculations. The only empirical data 
required in the integrations are the masses and radii of the planets. ‘The mass 
and mean radius of Jupiter are 

M = 1-902 x 10°g, R=69 900 km, (4) 
and of Saturn 

M =0:569 x 10%" g, R=57550km. (5) 
The masses of the planets are known accurately, but the radii are uncertain by 
about I per cent. The mean density of Jupiter is 1-33+0-04g/cm*, and of 
Saturn 0°71 + 0:02 g/cm’. 

3. Planetary models have been constructed on the assumption that the planets 
consist entirely of hydrogen, apart from central concentrations of heavy materials ; 
a numerical pressure—density relationship for hydrogen is contained in the earlier 
paper (x). The results of the calculations are presented in Tables I and II. 
In these tables r is the parameter introduced in equation (2), p and p are the 
corresponding pressure and density, and M(r) is the mass contained within the 
surface defined by r: 


M(r)= 47 p(s)s? dx. (6) 


The hydrogen content of Jupiter is estimated to be about 80 per cent by mass and 
of Saturn about 65 per cent by mass, but the precise proportions depend on the 
densities of the central cores of heavy material. If the density of the heavy 
material is assumed to be 14 g/cm? at the pressures at the centres of the planets, 
the estimated hydrogen contents of Jupiter and Saturn are 76 and 63 per cent 
by mass respectively. The computed proportions of hydrogen are larger the 
denser the heavy material is assumed to be, and in the limit of infinite density 
the estimated hydrogen contents are 84 and 69 per cent by mass. At these 
enormous pressures the heavy material is unlikely to be less dense than 14 g/cm’, 
which is only three times the corresponding density of hydrogen. ‘Thus, irres- 
pective of the dimensions of the central cores, the hydrogen contents compatible 
with the present planetary models are restricted to narrow limits, namely, 76 to 
84 per cent by mass for Jupiter and 63 to 69 per cent by mass for Saturn. These 
small uncertainties in the compositions of the planets are not in themselves of 
especial interest in this paper. Moreover, the uncertainties have little influence 
on the computed moments of inertia and polar flattenings of the planets; the 
dimensioris of the central cores are so small that the moments of inertia can be 
altered by only about 2 per cent. ‘The dimensions of the central cores are of 
even less importance in the planetary models considered later in the paper. 
It will therefore be assumed that the central concentrations of heavy elements 
are so dense that they may be treated as mass points. 

4. Planetary models have also been constructed on the assumption that each 
planet consists of a homogeneous mixture of hydrogen and helium, together 
with a dense central core of heavier elements. It has been shown in the earlier 
paper (x) that the density of a mixture containing m atoms of hydrogen to each 
atom of helium is 

p={(n+4)/(n+4o0)}on, (7) 
where py is the density of hydrogen at the pressure concerned and where « is 
the ratio of the densities of helium and hydrogen. ‘The value of « is 4 in the 
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non-metallic outer layers of the planets, and it is 2-5 in the deep interiors where 
hydrogen is metallic. Equation (7) gives the pressure—density relationship of 
any mixture of hydrogen and helium in terms of the data for hydrogen, and this 
relationship has been used in constructing the planetary models. It has been 
TaBLe | 
Estimated Density Distribution in Fupiter 
(Model J,—heavy elements assumed concentrated at centre) 
r p p M(r) 
(10° cm) (10!? dyne/cm*) (g/cm*) (1077 g) 
° 009 1902 
0°03 o'17 1895 
o'l4 0°24 1871 
0°28 0°28 1842 
0°45 o°31 1812 
0°64 0°33 1781 
080 0°35 1757 
0-80 0°77 1757 
2°10 0o'99 1595 
4°08 1°26 1401 
6°71 1°54 1203 
10°10 1°84 1012 
14°4 2°17 835 
19°9 2°53 680 
27°2 2°93 551 
37°7 3°43 45! 
54°8 4°13 379 
93 5°8 331 
Rea ne 303 


TABLE II 
Estimated Density Distribution in Saturn 
(Model S,—heavy elements assumed concentrated at centre) 

r p p M(r) 
(10° cm) (10!? dyne/cm?) (g/cm*) (107? g) 

5°755 ° 0°09 569 

5°6 0°02 O'rs 562 

0°06 0°20 548 

ovll 0°23 534 

0°17 0°25 518 

0°36 0°30 479 

0°62 0°33 443 

080 0°35 425 

o*80 0°77 425 

ris 0°83 397 

2°19 1°O1 336 

3°69 1'21 284 

6-o1 1°48 241 

10°! 1°85 209 

20°0 2°53 188 
62:2 4°39 178 

om ne 176 

shown in the previous paper (1) that the mean densities of the planets can be 
explained without central concentrations of heavy material if n is assumed to be 
13°0 for Jupiter and 6-4 for Saturn. Smaller helium contents—that is, larger 
values of n—have been adopted in the present models and, as before, the masses 


17 
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of the central cores of heavy material have been determined by integrating the 
hydrostatic equation (3) inward from the surface of the planet. The values of 
n adopted for Jupiter are 45 and 22, and for Saturn 18 and g. The results of 
the calculations are presented in Tables III to VI. The masses of the central 
cores of heavy material are smaller the smaller the value of n, and the various 


TABLE III 
Estimated Density Distribution in Fupiter 
(Model J,; n=45) 


r p p M(r) 
(10° cm) (10!” dyne/cm?*) (g/cm*) (107? g) 

° 0°09 1902 
0°03 o'1g 1895 
Orls 0°26 1868 
0°30 0°30 1837 
0°48 0°33 1804 
0°69 0°36 1770 
0-80 0°37 1754 
080 o'81 1754 
2°33 1°08 1558 
4°41 1°36 1347 
7°11 1°66 1133 
10°47 1°97 928 
14°6 2°29 740 
19°6 2°63 577 
25°8 3°00 444 
34°1 3°43 341 
46°7 4°00 271 
7°°7 4°93 227 
eee eee 201 


TABLE IV 
Estimated Density Distribution in Jupiter 
(Model J,; n=22) 


r p p M(r) 
(10° cm) (10? dyne/cm*) (g/cm*) (107? g) 
6°99 ° o'10 1902 
0°03 0°20 1894 
0°16 0°28 1865 
0°32 0°33 1832 
0°52 0°36 1797 
o'74 0°39 1761 
0°80 0°40 1752 
080 0°85 1752 
2°56 z°R7 1519 
4°74 1°48 1291 
7°50 1°79 1060 
10°81 2°10 839 
14°7 2°42 640 
19°1 2°73 469 
24° 3°05 332 
29°9 3°39 229 
37°4 3°78 161 
49°7 4°35 120 
eve dies 99 
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planetary models lead to nearly the same hydrogen contents. The estimated 
hydrogen contents are smallest if the planets are assumed to be chemically homo- 
geneous and largest if all the heavier elements are assumed to be concentrated at 
the centre. The extreme estimates for Jupiter are 76 and 84 per cent by mass, 
and for Saturn 62 and 69 per cent by mass. Thus the spatial distributions of the 
heavier elements within the planets have little influence on the estimated chemical 
compositions. ‘The location of the heavier materials affects primarily the computed 
moments of inertia of the planets and the polar flattenings due to their rotations. 
The moment of inertia J is greatest if the heavier elements are distributed 
uniformly and least if all the heavier materials are concentrated at the centre. 
The extreme values of the ratio // MR? are 0-223 and 0-282 for Jupiter and 0-189 
and 0-285 for Saturn, cf. Tables VII and VIII. The values of J/ MR? derived 
from the polar flattenings of Jupiter and Saturn using the Radau—Darwin 
approximation are 0-25 and 0-22 respectively, but the approximation may be in 
error by I0 per cent. 


‘TABLE V 


Estimated Density Distribution in Saturn 
(Model S,; n=18) 


r p p M(r) 
(10° cm) (10% dyne/cm?) (g/cm*) (107? g) 
"755 ° o°10 569 
6 0°02 o'19g 560 

0°07 0°24 544 
o'13 0°28 526 
0°20 0°30 507 
- 0°43 0°35 460 
o'72 0"40 418 
080 0°40 409 
0°80 0°86 409 
1°43 0°99 349 
2°47 1°17 277 
3°82 1°38 217 
5°69 1°61 169 
8°63 1°93 134 
14°7 2°46 113 
39°1 3°90 102 
ane nn 96 


UONGMNEONDWMOOUNdGNSE 


5 
5 
5 
5 
5 
4 
4 
3 
3 
3 
3 
2 
2 
I 
I 
° 
° 


5. The oblateness of a rotating planet can be calculated on the assumption | 
that the interior is in hydrostatic equilibrium. ‘l'his assumption is certainly 
valid in massive planets such as Jupiter and Saturn. Deviations from hydro- 
static equilibrium cannot exceed the strength of the material, which is of the 
order of magnitude of 10° dyne/cm? for most solids. ‘This pressure is negligible 
by planetary standards and could not introduce an appreciable error. Even for 
the Earth the ellipticity computed on the assumption of hydrostatic equilibrium 
is not in error by more than 0-5 per cent (§). A complete theoretical treatment of 
the figure of a planet requires a knowledge of the internal density distribution 
which was formerly not available for Jupiter and Saturn. ‘The moments of 
inertia of the planets had therefore to be derived from the empirical ellipticities 
by means of the approximate Radau-Darwin formula. ‘This approximation is 
reliable if the ratio // MR? is not much smaller than 0-4, but it can be seriously in 
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error when J/ MR? is as small as for Jupiter and Saturn; cf. reference (6). Also, 
second and higher powers of the ellipticity are ignored in the Radau—Darwin 
formula, which is not justified for Jupiter and Saturn. An accurate numerical 
treatment will therefore be used in preference to the Radau—Darwin approxi- 
mation. The ellipticities of the planets will be computed on the basis of the 
various planetary models in the hope that some of these models can be eliminated 
by the empirical data. ‘The empirical data for Jupiter are 
€,*=15'4t0°2, &,/m=0°73 + 0°02, (8) 
and for Saturn 
e, '=9°5+0°2, e,/m = 0°65 + 0°02, (9) 
é, being the surface ellipticity and m the ratio of the centrifugal and gravitational 
forces on the equator. 
Tas.e VI 
Estimated Density Distribution in Saturn 
(Model S;; n=9) 
r p p M(r) 
(10° cm) (10!? dyne/cm?) (g/cm") (10?? g) 
5°755 fo) or12 569 
5°6 0°03 0°22 559 
5°4 009 0°28 540 
5°2 O'ls 0°32 519 
5"0 0°23 0°35 497 
4°5 0°48 o'41 443 
4°00 0°80 0°45 394 
4°00 0°80 0°94 394 
3°5 1°61 I°I2 303 
3°0 2°57 1°30 219 
2°5 3°66 1°48 153 
2°0 4°93 1°66 103 
1°5 6°54 1°87 69 
1'o 9°08 2°16 49 
o's 16°6 2°84 40 
° oe pes 39 
6. Darwin (7) and de Sitter (8) have given a comprehensive treatment of the 
figure of a planet which is more than adequate for the present purposes; the 
hydrostatic equations in the interior are treated accurately and only third and 
higher powers of the ellipticity are ignored. Bullard (§) has reduced this treatment 
to an entirely numerical procedure which is applicable if the internal density 
distribution is known; this method has been adopted in the present calculations. 
The equation of a surface of constant pressure is of the form 


B =b{1 —esin® ¢ + ($e? — x) sin? 24}, (10) 
where £ is the distance from the centre of the planet and ¢ the angle of latitude ; 
e is the ellipticity of the surface and the quantity « allows for the departure from an 
ellipsoid. As in equation (3), it is convenient to introduce the radius r of the 
sphere which contains exactly the same volume as the surface (10). ‘The theory 
is most simply formulated in terms of the quantity 

e' =e — 5e*/42 + 4x/7 (11) 
and its logarithmic derivative 
» = d(in e’)/d(inr). (12) 
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De Sitter (8) has shown that 7 satisfies the equation 


(1 +n)!2=(5/pr°)| F(n)pr*ar, (13) 
~0 


where p is the mean density within the surface defined by r and where 

F(y) =(1+ 49 - jou? + A)/(t+ 0), (14) 
A being a small quantity proportional to the ellipticity. ‘lhe importance of A 
in equation (14) has been examined and its omission alters the calculated 
ellipticities by less than 0-3 per cent; it has therefore been ignored in these 
calculations. For values of » less than unity the function F(») Goes not deviate 
from unity by more than 1 per cent. In the Radau-Darwin approximation 
F(n) is replaced by unity and equation (13) reduces to 

(1+)'* =${1— 3(/, MR’)}, (15) 

I being the planet’s moment of inertia. This approximation is not accurate 
for Jupiter and Saturn as the values of 7 are large, F(y) being up to 7 per cent 
less than unity. Equation (13) has been integrated numerically using the internal 
density distributions derived above; the boundary condition is that » must 
vanish at the centre of the planet. ‘The surface values of » (denoted by »,) for 
the various planetary models are given in Tables VII and VIII. ‘The magnitude 
of the parameter », is all that need be known about the internal density distri- 
bution so far as the figure of the planet is concerned. For a given value of », 
and a given angular velocity w the ellipticity e, of the planet can be found from the 
boundary condition at the surface 


€, 1, = 5/2 — 2, + 1Op?/2I + 4e,?/7 — Ope, /7, (16) 

where e,’ is the value of e’ at the surface and where 
y= (w®R5/ MG), (17) 
M and R being the mass and mean radius of the planet. Equation (16) involves 


the quantity « indirectly through the parameter e,’, cf. equation (II). De Sitter 
(8) has set an upper limit to the magnitude of « at the surface 

|x|<$e,u — te,?, (18) 
and even this extreme value leads to an uncertainty of only I per cent in the 
computed ellipticities of the planets. ‘The departure from an ellipsoid is therefore 
too small an effect to be included in the present calculations and « has accordingly 
been put equal to zero. The results of the calculations are presented in ‘Tables 
VII and VIII. 

7. The agreement between the computed and empirical ellipticities of the 
planets cannot be expected to be better than 4 per cent. Not only are the 
empirical ellipticities uncertain by about 2 per cent, but the ratio » defined by 
equation (17) may be in error by 2 or 3 per cent because of inaccuracies in the 
adopted values of the radii. The ratio /; MR? in best agreement with the 
empirical ellipticity is 0-267 for Jupiter and 0-262 for Saturn; these values are 
respectively 8 and 18 per cent larger than those indicated by the Radau—Darwin 
formula. The empirical data can determine the ratios /;MR?® only to about 
5 percent. ‘The probable masses of the central concentrations of heavy materials 
are 11+ 10 Mg, for Jupiter and 5 +4 M, for Saturn; the mass My, of the Earth is 
5°976 x 10°7g. The models J, and S, give the best representations of the planets, 
but the chemically homogeneous models J, and S, cannot be completely excluded. 
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8. The present calculations were undertaken partly in an attempt to establish 
or refute the condensation theories of the origin of the solar system. This mode 
of formation of the planets would require the central concentrations of heavy 
materials in Jupiter and Saturn to be at least as massive as Neptune, that is, at 
least 17 My; cf. reference(1). ‘The above calculations show that this requirement 
is compatible with the empirical data for Jupiter, but the central concentration 
in Saturn seems to be definitely less massive than Neptune. The discrepancy 
cannot be attributed to the assumption that the planets are chemically homo- 
geneous, apart from the central concentrations of heavy materials; the estimated 
masses of the central cores are decreased if the mean atomic weight is assumed 
to increase gradually with depth. But the condensation theories of the origin 
of the solar system cannot be dismissed merely on account of this discrepancy 
for Saturn. The axial rotation of Saturn, and of Jupiter to a lesser degree, 
exhibit features which are not understood (g). The period of rotation adopted 
for Saturn is 10"38™, but the apparent period is 10%14™ at the equator and 
about 11 in high latitudes. The variation of period with latitude is presumably 
due to circulation currents in the atmosphere but, in contrast to the Earth, the 
apparent period of rotation increases with latitude. ‘The estimated mass of the 
central core in Saturn would be close to the mass of Neptune if the period 
indicated by equatorial spots were adopted. The observational data cannot be 
interpreted unambiguously until the nature of the atmospheric circulations 1s 
known. Should more reliable estimates of the parameter . defined by equation 
(17) become available, the ellipticities corresponding to the various planetary 
models can easily be found from equation (16) using the values of y, given in 
Tables VII and VIII. Second-order perturbations of the planets on their 


closest satellites may also provide a check on the planetary models (10). 
Taste VII 


Computed Ellipticity of Jupiter * 
Ns i" I/MR* 
‘67 17°4 0'223 
"50 16°6 O'241 
"32 15°7 0'260 
"10 14°7 0282 


Tasie VIII 
Computed Ellipticity of Saturn * 


Model "Ns e,? I/MR* 
1.97 11°5 0189 

1°67 10°6 0'220 

S, 1°34 9°68 0'256 
S, 1°04 8-79 0°285 

* In models Jg and Sq the planets are assumed to be chemically homogeneous; the estimated 
internal density distributions are contained in the earlier paper (1). 

The aggregation theories of the origin of the solar system relate the chemical 
compositions of Jupiter and Saturn in an interesting way. On these theories 
one would expect the central mass /, of heavy materials to be the same for the 
two planets. An upper limit to the mass M, is set by the following consideration. 
The value of n for Jupiter cannot be less than that for Saturn; that is, the material 
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captured by the larger planet must be at least as rich in hydrogen as that captured 
by the smaller planet. The value of n indicated by the mean density is 23 for 
both planets if M, is assumed to be 19 My: 

M,=19 Mg, n=23. (19) 
Larger values of M, would lead to larger values of n for Saturn than for Jupiter. 
On the other hand, M, must be at least 17 My, the mass of Neptune. The mass 
of the central concentrations of heavy materials is therefore restricted to narrow 
limits and the value of is nearly the same for the two planets. In other words, 
the difference in the hydrogen contents of Jupiter and Saturn is due mainly to 
the fact that the central core accounts for a larger fraction of the mass of the 
smaller planet. The value of m in relations (19) corresponds to 85 per cent 
hydrogen by mass and this is presumably the hydrogen content of the material 
from which the planets were formed. ‘This relationship between the chemical 
compositions of Jupiter and Saturn is a welcome simplification ; but it cannot be 
regarded as a proof that the planets were formed by an aggregation process. 
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